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The  focus  of  this  dissertation  is  on  the  design  and 
performance  analysis  of  multiuser  receivers  for  asynchronous 
code  division  multiple  access  communication  channels.  In 
many  applications,  such  as  mobile  radio  networks,  the  power 
levels  of  the  received  signals  at  the  central  base  station 
receiver  are  unequal,  resulting  in  stronger  users 
interfering  with  weaker  users.  Historically,  adaptive  power 
control  has  been  suggested  as  a solution;  more  recently, 
however,  multiuser  detectors  have  been  demonstrated  to 
operate  independently  of  the  received  power  levels  and  offer 
better  system  performance  than  could  be  obtained  through  the 
use  of  power  control  alone.  One  particular  multiuser 
detection  method,  the  decorrelating  strategy,  has  the 
advantage  over  other  methods  of  not  requiring  the  side 
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information  of  the  received  power  levels.  The  decorrelating 
strategy  suffers,  however,  from  long  decoding  delays,  huge 
memory  requirements,  and  the  burden  of  computing  the  inverse 
of  an  enormous  matrix,  all  of  which  render  it  impractical. 

In  this  work,  a method  of  implementing  a decorrelating 
receiver  which  does  not  suffer  from  long  decoding  delays, 
huge  memory  requirements,  or  the  inversion  of  an  enormous 
matrix  is  proposed  for  both  additive  white  Gaussian  noise 
channels  and  f requency-nonselective  Rayleigh  fading 
channels . 

A new  performance  measure,  the  near- far  spread  of  the 
receiver,  is  introduced.  The  performance  of  the  proposed 
receiver  is  analyzed  for  both  the  additive  white  Gaussian 
noise  and  f requency-nonselective  Rayleigh  fading  channels. 
Numerical  examples  successfully  demonstrate  the  receiver's 
immunity  to  unequal  received  energies.  The  sensitivity  of 
the  proposed  receiver  to  errors  in  the  timing  information  it 
requires  is  examined,  as  well  as  the  computational 
complexity  associated  with  its  operation. 
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CHAPTER  1 
INTRODUCTION 

Wireless  Multiple  Access  Communication  Systems 

Wireless  multiple  access  communication  systems  are 
among  the  fastest  growing  technologies  in  the  field  of 
communications.  Cellular  telephones,  cordless  telephones, 
and  pagers  have  become  commonplace  items,  seen  almost 
everywhere.  Indoor  wireless  communication  networks  and 
personal  communication  networks  are  becoming  more  popular  in 
the  business  and  public  arenas.  A typical  wireless  multiple 
access  communication  system  is  illustrated  in  Figure  1.1. 

The  system  is  characterized  by  multiple,  possibly  mobile, 
users  sharing  a common  channel.  Each  user  has  a wireless 
link  to  a common,  centralized  base  station.  The  base 
station  is  typically  connected  to  the  existing  telephone 
infrastructure . 

One  of  the  first  decisions  which  must  be  made  when 
designing  any  multiple  access  communication  system  is  how 
the  users  will  be  allowed  to  share  the  channel.  Currently 
there  are  three  competing  technologies  for  this  task: 
frequency  division  multiple  access  (FDMA) , time  division 
multiple  access  (TDMA) , and  code  division  multiple  access 
(CDMA) . In  FDMA,  the  available  bandwidth  is  divided  into 
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narrow,  nonoverlapping  channels.  Multiplexing  is  achieved  by 
assigning  each  user  to  a different  channel.  In  TDMA,  time 
is  divided  into  nonoverlapping  slots,  and  multiplexing  is 
achieved  by  assigning  each  user  to  a different  time  slot. 

For  the  duration  of  each  time  slot,  the  corresponding  user 
has  access  to  the  entire  available  bandwidth.  In  CDMA,  all 
users  share  the  entire  available  spectrum  at  all  times.  The 
focus  of  this  work  will  be  on  direct  sequence  code  division 
multiple  access  (DS-CDMA) , in  which  multiplexing  is  achieved 
by  assigning  a unique  signature  waveform  to  each  user. 
Knowledge  of  a particular  user's  signature  waveform  is 
normally  required  in  order  to  decode  that  user's  data  from 
the  resulting  multitude  of  signals. 
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Recently,  CDMA  has  received  much  attention  as  an 
attractive  alternative  to  FDMA  and  TDMA  for  urban  digital 
mobile  radio  networks  [1]— [16] . Several  papers  suggest  that 
a CDMA  system  could  support  more  simultaneous  users  than  an 
FDMA  or  a TDMA  system  [1],  [2],  [7],  [8].  One  reason  for 

the  increased  capacity  in  voice  networks  is  that  users  are 
speaking  only  roughly  35%  of  the  time;  the  rest  of  the  time 
they  are  listening  [1] . This  voice  activity  cycle  thus 
serves  to  reduce  the  overall  mutual  interference.  It  has 
also  been  suggested  that  a CDMA  system  could  overlay  across 
the  same  frequency  bands  currently  being  used  for  other 
purposes  without  the  performance  of  either  system  being 
affected,  thus  saving  valuable  spectrum  [3],  [17]— [24] . An 

important  parameter  in  cellular  networks  is  the  frequency 
reuse  factor,  which  dictates  how  often  the  same  carrier 
frequency  can  be  used  across  adjacent  cells.  Typically, 

FDMA  and  TDMA  cellular  systems  repeat  the  same  frequency  in 
only  one  of  every  seven  cells  [2].  For  cellular 
applications,  CDMA  allows  for  a frequency  reuse  factor  of 
one  since  every  cell  can  use  the  same  frequencies.  In 
multipath  fading  channels,  multiple  reflections  of  the 
transmitted  signal  can  arrive  at  the  receiver,  resulting  in 
destructive  interference.  Since  DS-CDMA  signals  are 
wideband,  they  have  the  capability  of  resolving  multipath 
signal  components  which  can  then  be  recombined  to  provide 
improved  performance  in  multipath  fading  channels  [25] -[28]. 


4 


Review  of  DS-CDMA  Systems 

Before  proceeding,  it  will  be  helpful  to  review  the 
architecture  of  a direct  sequence  code  division  multiple 
access  communication  system.  Figure  1.2  illustrates  the 
system  model.  There  are  K users  sharing  an  additive  white 


s R(  t - i T ) 


Figure  1.2  - Synchronous  DS-CDMA  System  Model 


Gaussian  noise  (AWGN)  channel.  Here  the  users  are  assumed 
to  employ  binary  phase  shift  keying  (BPSK)  modulation.  Each 
user  transmits  on  the  same  carrier  frequency  denoted  by  coc  ; 
multiplexing  is  achieved  by  assigning  a unique  signature 
waveform  to  each  user.  The  signature  waveform  of  the  k 
user  can  be  written  as 


5 


Sv(t)=* 


J|’Xa^PTc(t-(3-l)Tc)cos(tOct  + 0k),  0 < t < T 

V l j=i  (1.1] 


0, 


otherwise 


where 

ake{±l},  j = l,2,-,N  (1.2) 

th 

is  the  signature  sequence  of  the  k user,  pT  is  the  pulse 

Ac 

shape  of  the  signature  chip  waveform  whose  duration  is 
Tc  = T/N  seconds  where  T is  the  common  bit  duration,  and  0k 
is  the  carrier  phase  of  the  kth  user.  Note  that  in  this 
model  the  signature  waveforms  are  normalized  to  unit  energy 
and  are  nonzero  only  during  the  time  interval  [0,T)  . Each 
user  is  assumed  to  transmit  a message  of  length  M bits,  and 
for  simplicity  the  system  is  assumed  to  be  synchronous;  that 
is,  the  ith  bits  of  all  users  arrive  simultaneously  at  the 
base  station.  The  more  general  asynchronous  system,  in 
which  each  user  may  experience  a different  propagation 
delay,  will  be  considered  in  the  remaining  chapters.  The 
received  signal  at  the  base  station  can  be  written  as 

M K 

r(t)  = X Z wk(i)bk(i)sk(t  - iT  - t)  + n(t)  (1.3) 

i=l  k=l 

where  wk(i)  is  the  received  energy  of  the  ith  bit  of  the  kth 
user,  bk(i)e{±l}  is  the  ith  transmitted  bit  polarity  of  the 
kth  user,  T is  the  common  propagation  delay  from  each  user 
to  the  base  station,  and  n(t)  is  additive  white  Gaussian 
noise  with  power  spectral  density  a2  . 

The  job  of  the  base  station  receiver  is  to 
simultaneously  demodulate  the  data  for  each  of  the  K users. 
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It  was  shown  by  Verdu  [29]  that  a sufficient  statistic  for 
determining  the  maximum  likelihood  transmitted  sequences 
bk(i)  for  k = l,2,--,K  and  i = l,2,*-*,M  can  be  obtained  by 
passing  the  received  signal  r(t)  through  a bank  of  filters 
matched  to  the  allowed  signature  waveforms.  The  matched 
filter  bank  is  shown  in  detail  in  Figure  1.3.  The  received 
signal  is  correlated  with  appropriately  delayed  replicas  of 


SR  (t-iT-T) 

Figure  1.3  - Matched  Filter  Bank  for  Synchronous  System 

each  of  the  allowed  signature  waveforms  and  the  results  are 
then  sampled  at  the  end  of  each  bit  interval.  Define  the  K 
by  K element  system  cross-correlation  matrix  R such  that  its 
(m,n)th  element  is 

= Jsm(t)sn(t)dt 


R 


(1.4) 
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The  output  of  the  kth  matched  filter  sampled  at  the  end  of 


the  ith  bit  interval  is  then 


(i+l)T+T 


Yk(i)  = | r(t)sk(t  — iT  — x)dt 


(1.5) 


j=i 

i*k 


where 


nk  (i)  = J n(  t)  sk  (t  - iT  - t)  dt 


(1.6) 


is  the  noise  component  out  of  the  matched  filter.  Finally, 
the  vector  of  matched  filter  outputs  sampled  at  the  end  of 


To  understand  the  nature  of  the  near- far  problem,  it  is 
necessary  to  first  understand  the  conventional  method  of 
demodulating  the  received  signal.  The  conventional  DS-CDMA 
receiver  consists  of  the  matched  filter  bank  followed  by  a 
symbol-rate  sampler  and  threshold  device  as  shown  in  Figure 
1.4.  The  ith  input  to  the  kth  user's  threshold  device  will 
thus  be  given  by  equation  (1.5).  The  performance  of  the 
conventional  DS-CDMA  receiver  has  been  extensively  analyzed 
and  many  different  bounds  and  approximations  for  its 
probability  of  error  have  been  found  [30] -[33].  Note  that 
the  decision  statistic  yk(i)  is  the  sum  of  three  distinct 
components:  a desired  signal  component,  a multiple-access- 


the  ith  bit  interval  may  be  written  as 


Y = RWb  + n 


(1.7) 


where  Y = [y^i)  ,y2(i)  ,--,yk(i)]T  , b = [bx(i)  ,b2(i)  ,---,bk(i)]T  , 
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r ( t ) 


Figure  1.4  - The  Conventional  DS-CDMA  Detector 


interference  (MAI)  component,  and  a noise  component.  It  is 
the  presence  of  the  MAI  component  which  accounts  for  the  so- 
called  near-far  problem.  Since  in  practice  the  chosen 
signature  waveforms  will  not  have  zero  cross-correlations 
(particularly  in  the  more  general  asynchronous  case) , the 
magnitude  of  the  MAI  component  will  increase  as  the  received 
energies  of  the  K-l  other  users  increase.  As  a result,  even 
for  low  signature  waveform  cross-correlations,  the  MAI 
component  can  dominate  the  desired  user's  decision  statistic 
when  the  received  energy  of  any  other  user  becomes  much 
greater  than  that  of  the  desired  user.  Thus,  weaker  users 
will  not  be  able  to  communicate  in  the  presence  of  stronger 
users.  Since  unequal  received  energies  could  occur  at  the 
base  station  if  some  users  were  farther  away  than  others, 
this  problem  has  been  named  the  near- far  problem. 
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Until  recently,  the  most  popular  strategy  for  combating 
the  near-far  problem  has  been  the  user  of  transmitter  power 
control,  in  which  the  base  station  continually  monitors  the 
received  power  levels  of  all  users  and  instructs  them  to 
either  increase  or  decrease  their  transmitter  powers  so  that 
all  received  signals  are  of  approximately  equal  strength  at 
the  base  station  [1],  [34] -[36].  There  are  two  major 

disadvantages  to  using  a power  control  scheme,  however. 
First,  the  rate  at  which  power  control  updates  could  be 
performed  is  an  open  issue.  While  it  may  be  feasible  to 
counteract  differences  in  received  energies  due  to  radial 
distances  from  the  base  station,  a much  greater  update  rate 
would  be  required  to  combat  the  effects  of  fading. 

Secondly,  and  more  importantly,  the  use  of  power  control 
will  not  eliminate  the  near- far  problem  but  will  only  reduce 
it.  That  is,  even  if  perfect  power  control  was  possible, 
there  would  still  be  an  MAI  component  present  in  each  user’s 
decision  statistic.  The  presence  of  the  MAI  component 
results  in  an  irreducible  error  floor;  consequently,  the 
conventional  detector  can  never  obtain  the  same  level  of 
performance  that  a single  user  would  have  in  the  same 
channel . 


Previous  Work 


As  an  alternative  to  using  power  control  schemes, 
recent  work  has  focused  on  using  information  known  about  the 
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signature  waveforms,  time  delays,  received  energies, 
received  carrier  phases,  and  channel  parameters  to 
completely  eliminate  the  MAI  component  from  each  user's 
decision  statistic.  The  resulting  structures  are  called 
multiuser  detectors,  and  they  are  immune  to  the  near- far 
problem.  All  multiuser  detectors  have  a front-end 
consisting  of  the  matched  filter  bank;  the  difference 
between  the  various  proposed  multiuser  receivers  is  the 
algorithm  chosen  to  follow  the  filter  bank. 

The  optimum  multiuser  detector  was  derived  by  Verdu 
[29];  it  consists  of  the  matched  filter  bank  followed  by  a 
dynamic  programming  algorithm  which  searches  for  the  set  of 
maximum  likelihood  transmitted  sequences,  for  which  the 
Viterbi  algorithm  may  be  used.  Although  this  receiver  is 
able  to  find  the  maximum- likelihood  transmitted  sequences, 
there  are  several  practical  limitations  associated  with  its 
implementation.  First,  the  number  of  computations  required 
per  demodulated  data  bit  is  exponential  in  the  number  of 
active  users.  For  instance,  for  BPSK  transmission  the 
complexity  is  0(2K),  where  K is  the  number  of  active  users. 
Furthermore,  the  optimum  multiuser  detector  requires  the 
perfect  knowledge  of  the  received  energies  and  carrier 
phases  of  all  bits.  Perfect  knowledge  of  each  user's 
propagation  delay  is  also  required.  The  acquisition  of  this 
side  information  in  the  hostile  urban  digital  mobile  radio 
environment  may  be  extremely  difficult. 
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Motivated  by  the  need  to  reduce  the  complexity  of  the 
optimum  multiuser  detector,  Lupas  and  Verdu  considered 
restricting  the  receiver  to  have  a complexity  which  is 
linear  in  the  number  of  users  and  arrived  at  the  so-called 
decorrelating  detector  [37],  [38].  A linear  detector  will 

form  its  decisions  as 


where  L is  a linear  transformation.  It  is  interesting  to 
note  that  the  conventional  detector  is  a linear  detector 
which  simply  chooses  L = I (the  identity  matrix)  as  its 
transform.  The  decorrelator  chooses  L = R_1  as  its 
transformation  and  therefore  forms  its  decisions  as 


= sign[wb  + R xnj 

It  is  clear  from  equation  (1.9)  that  the  decorrelator 
would  make  no  errors  in  the  absence  of  noise  since 


MAI  component  present  in  equation  (1.9);  consequently,  the 
decorrelator's  decision  is  independent  of  the  received 
energies  of  the  other  users  and  hence  the  decorrelator  is 
immune  to  the  near- far  effect.  It  should  be  mentioned  that 
the  inverse  of  the  matrix  R is  guaranteed  to  exist  provided 
that  the  signature  waveforms  are  linearly  independent  of  one 
another,  a rather  mild  restriction  compared  to  that  of  the 
conventional  DS-CDMA  detector  which  would  require  the 
signature  waveforms  to  be  orthogonal  in  order  to  achieve 


b = signjLY] 


(1.8) 


(1.9) 


Furthermore,  note  that  there  is  no  longer  an 
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roughly  the  same  level  of  performance.  One  advantage 
offered  by  the  decorrelator  is  that  it  does  not  require 
knowledge  of  the  received  energies;  however,  for  an 
asynchronous  system  it  does  require  knowledge  of  the 
propagation  delays  and  received  carrier  phases  of  all 
received  bits.  Furthermore,  in  an  asynchronous  system  the 
decorrelator  must  observe  the  entire  transmitted  sequences 
of  all  users  before  it  can  begin  decoding.  Thus,  the 
decorrelator  will  suffer  from  long  decoding  delays  and  large 
memory  requirements.  Also,  the  inverse  of  an  MK  by  MK 
matrix  must  be  computed,  where  M is  the  message  length  in 
symbols  and  K is  the  number  of  active  users.  Thus,  the 
decorrelator  is  rendered  impractical  for  most  asynchronous 
systems.  As  an  alternative  implementation,  a linear  time- 
invariant  transfer  function  was  also  found  for  the 
decorrelator  [38],  but  this  approach  cannot  accommodate 
changes  in  the  timing  configuration  or  the  number  of  active 
users.  Finally,  it  is  interesting  to  note  that  when 
demodulation  of  only  one  user's  data  is  required,  as  in  the 
mobile  receivers,  for  example,  then  the  matched  filter  bank 
may  be  replaced  with  only  one  matched  filter.  Instead  of 
correlating  the  received  signal  with  a replica  of  the 
desired  user's  signature  waveform,  the  received  signal  is 
correlated  with  the  projection  of  the  desired  user's 
signature  waveform  onto  the  subspace  which  is  orthogonal  to 
the  space  spanned  by  the  signature  waveforms  of  the  other 
active  users. 
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Nonlinear  multiuser  detection  techniques  have  also  been 
proposed.  These  receivers  generate  tentative  decisions  for 
each  user  and  later  use  them  for  cancellation  of  the  MAI 
[39]  — [41]  . The  multistage  multiuser  detector  [39]  utilizes 
the  conventional  DS-CDMA  detector  to  form  the  initial 
polarity  decisions  in  the  first  stage.  These  decisions  are 
then  used  to  reconstruct  the  MAI  contributed  from  each  user. 
The  MAI  contributions  are  summed  together  and  subtracted 
from  the  desired  user's  initial  decision  statistic  in  the 
second  stage.  At  this  point,  final  decisions  could  be  made, 
or  additional  stages  of  MAI  re-creation  and  subtraction 
could  be  performed.  The  decorrelator  has  also  been 
suggested  as  a replacement  for  the  conventional  detector  in 
the  first  stage.  The  main  disadvantage  associated  with 
these  nonlinear  detectors  is  that  in  order  to  correctly 
reconstruct  the  MAI,  they  require  knowledge  of  all  the 
received  energies,  carrier  phases,  and  propagation  delays  of 
all  bits  from  all  users. 

Another  nonlinear  technique  proposed  for  multiuser 
decoding  is  to  first  decode  the  strongest  user,  then  use 
that  decision  to  help  decode  the  next  strongest  user,  and  so 
on  [42]  — [44]  . The  idea  is  that  the  strongest  user  will  have 
the  least  amount  of  MAI  corruption  present  in  the  decision 
statistic,  and  so  a conventional  decision  may  be  made.  Then 
the  MAI  from  the  strongest  user  may  be  re-created  and 
subtracted  from  the  initial  decision  statistic  of  the  next 
strongest  user.  This  process  then  repeats  itself  until  the 
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weakest  user  has  been  demodulated.  Such  a strategy  also 
requires  knowledge  of  the  received  energies,  carrier  phases, 
and  propagation  delays.  Interestingly,  the  receiver 
presented  by  Dent  [44]  performs  the  MAI  cancellation  by 
nulling  a component  in  a spectral  domain  and  hence  avoids 
the  need  to  reconstruct  the  MAI  in  the  time  domain. 

Others  have  suggested  an  adaptive  filter  based  receiver 
structure  for  multiuser  detection  [45] -[47].  Both  a minimum 
mean  squared  error  criterion  and  a weighted  least  squares 
criterion  have  been  investigated  [45] . Although  not  a 
multiuser  detector,  the  receiver  structure  presented  by 
Miller  [46]  is  able  to  obtain  timing  information  on  a single 
user  and  thus  can  be  used  for  joint  estimation/detection. 
Interestingly,  this  receiver  does  not  require  knowledge  of 
the  desired  user's  signature  waveform. 

The  extension  of  multiuser  receivers  to  noncoherent 
communication  systems  was  presented  by  Varanasi  [48],  [49]. 

Unlike  their  coherent  counterparts,  these  receivers  do  not 
require  knowledge  of  received  carrier  phases. 

The  extension  of  multiuser  detectors  to  fading  channels 
includes  [50]— [55] . The  linear  time- invariant  transfer 
function  version  of  the  decorrelator  was  presented  and 
analyzed  for  both  flat  Rayleigh  fading  and  frequency 
selective  fading  channels  by  Zvonar  [55];  as  with  their  AWGN 
counterpart,  these  receivers  cannot  accommodate  changes  in 
the  number  of  active  users  or  their  relative  delays. 
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Outline  of  the  Dissertation 

Motivated  by  the  promising  features  offered  by  the 
decorrelator,  this  work  proposes  a method  by  which  a 
decorrelating  detector  may  be  implemented  in  an  asynchronous 
communication  system  while  avoiding  the  inversion  of  the  MK 
by  MK  matrix  and  the  need  to  observe  and  store  the  entire 
transmitted  sequences  of  all  users.  One  strategy  to 
accomplish  this  goal  is  to  apply  the  decorrelating  procedure 
to  a very  small  portion  of  the  received  sequences  at  a time, 
thus  in  effect  windowing  the  received  sequences.  Such  a 
strategy  has  been  suggested  [45] , and  a sliding  window 
version  of  the  decorrelator  was  proposed  [56] . Neither  work 
presented  analytical  results  on  the  performance  of  such  a 
receiver,  however.  Furthermore,  the  sliding  window 
decorrelator  [56]  is  dependent  upon  the  system  employing  a 
rate  1/2  convolutional  code  for  error  correction  purposes, 
thereby  restricting  the  overall  applicability  of  such  a 
receiver.  The  main  contributions  of  this  work  are  to 
provide  an  analytical  treatment  of  system  performance  when  a 
windowed  decorrelating  strategy  is  used  and  to  propose  a 
method  for  implementing  such  a receiver  independently  of  any 
error  correction  schemes,  for  both  AWGN  and  frequency 
nonselective  Rayleigh  fading  channels. 

Chapter  two  presents  a windowed  decorrelating  receiver 
suitable  for  use  in  additive  white  Gaussian  noise  (AWGN) 
channels.  The  performance  of  the  receiver  is  analyzed  and 
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compared  to  that  of  the  conventional  detector  with  power 
control  and  the  full  length  decorrelator. 

Chapter  three  presents  a windowed  decorrelating 
receiver  suitable  for  use  in  frequency  nonselective  Rayleigh 
fading  channels.  The  performance  of  the  receiver  is 
analyzed  and  compared  to  that  of  the  conventional  detector 
with  power  control  and  the  full  length  decorrelator. 

Chapter  four  discusses  some  of  the  practical  issues 
associated  with  the  implementation  of  the  windowed 
decorrelator.  Computational  complexity  and  sensitivity  to 
errors  in  the  timing  information  are  investigated. 

Chapter  five  summarizes  the  windowed  decorrelator  in 
both  AWGN  and  flat  fading  channels.  Areas  for  further  work 
related  to  the  windowed  decorrelator  are  suggested. 


CHAPTER  2 

WINDOWED  DECORRELATING  MULTIUSER  RECEIVER  FOR  ADDITIVE  WHITE 

GAUSSIAN  NOISE  CHANNELS 

System  Description 

The  link  from  the  users  to  the  base  station  is  assumed 
to  be  asynchronous;  that  is,  no  effort  is  made  by  the  users 
to  coordinate  their  transmitted  sequences  in  time.  The  link 
from  the  base  station  to  the  users  is  assumed  to  be 
synchronous;  that  is,  the  base  station  would  transmit  the 
ifck  bit  of  each  user  simultaneously.  The  asynchronous 
transmission  makes  the  design  for  the  base  station  receiver 
much  more  challenging  than  the  receiver  for  each  of  the 
mobile  units.  To  indicate  why  this  is  so,  Figure  2.1 
displays  a three  bit  wide  window  in  both  a synchronous 
system  and  an  asynchronous  system.  The  window  spans  bits  i 
through  i+2 . In  the  synchronous  system,  there  is  no 
interference  from  bits  outside  the  window;  consequently,  the 
standard  decorrelating  strategy  may  be  applied  to  each 
consecutive  window  independently.  In  the  asynchronous 
system,  however,  bits  outside  the  window  overlap  with  bits 
within  the  window  thus  introducing  MAI  from  outside  the 
window.  For  instance,  in  Figure  2.1  the  (i  + 3)r<^  bits  of 
users  1 through  K-l  will  interfere  with  the  (i+2)nd  bit  of 
user  K.  Furthermore,  since  the  receiver  will  be  attempting 
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Figure  2.1  - Comparison  of  a Window  in  Synchronous  and 
Asynchronous  Systems 


to  determine  the  polarity  of  the  bits  contained  within  the 
current  window,  the  polarities  of  the  (i+3)rd  bits  of  users 
1 through  K-l  are  unknown.  Consequently,  the  standard 
decor relating  procedure  cannot  be  applied  independently  to 
each  window  and  a more  complex  demodulation  strategy  must  be 
chosen.  The  focus  of  this  work  will  therefore  be  to  propose 
a receiver  for  the  base  station  and  to  analyze  its 
performance . 

The  communication  system  of  interest  is  one  in  which  K 
users  employ  BPSK  transmission  through  a common  AWGN  channel 
through  use  of  direct  sequence  code  division  multiple 
access,  as  illustrated  in  Figure  2.2.  For  purposes  of 
comparison,  the  notation  used  by  Lupas  [38]  will  be  adopted. 
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s R(  t - i T ) 


Figure  2.2  - Asynchronous  DS-CDMA  System  Model 

Each  user  transmits  a sequence  of  length  M bits,  and  without 
loss  of  generality  the  users  are  ordered  according  to  their 
relative  time  delays  0 < x1  < x2  <•■•<  xK  < T , where  T is  the 

common  bit  duration.  The  received  signal  at  the  base 
station  is  written  as 

r(  t ) = s(t,  b)  + n(  t)  (2.1) 

where  n(t)  is  white  Gaussian  noise  with  power  spectral 
density  a2  and 

s(t,  b)=  (2.2) 

i=l  k=l 

where  bk(i)  € {±1}  is  the  kth  user's  ith  data  bit,  wk(i)  is 
the  kth  user's  received  power  in  the  ith  bit  interval,  and 
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sk(t)  is  the  user's  signature  waveform  which  is  zero 
outside  the  time  interval  [0,T]  and  normalized  to  unit 
energy.  The  K by  K signature  waveform  cross-correlation 
matrices  are  defined  by  their  (j/k)^  component 

oo 

Rjk(i)=  Jij(t-X.)sk(t  + iT-Tk)dt  (2.3) 

— oo 

From  equation  (2.3)  it  can  be  seen  that  R(i)  = 0 for  |i|  > 1 and 
that  R(— 1)  = Rt(1)  . 

It  is  important  to  note  that  the  signature  waveforms 
are  RF  waveforms  and  thus  each  user's  RF  carrier  waveform  is 
contained  within  each  user's  signature  waveform.  The  base 
station  is  assumed  to  have  knowledge  of  each  user's 
signature  waveform;  thus,  it  is  implied  in  this  model  that 
the  received  carrier  phase  of  each  user  is  known  by  the  base 
station.  This  assumption  will  be  justified  for  BPSK 
modulation  since  knowledge  of  the  received  carrier  phases  is 
required  for  coherent  demodulation. 

Another  important  observation  from  equation  (2.3)  is 
that  baseband  transmissions  will  maximize  the  signature 
waveform  cross-correlations.  To  see  why  this  is  so, 
equation  (2.3)  can  be  re-written  with  baseband  signature 
waveforms  modulating  RF  carrier  waveforms  as 

oo 

Rj,k<i)=  Js.(t-l.)cos(wct  + e.)sk(t  + iT-ik)cos(wct  + ek)dt  (2.4) 
where  sk(t)  now  represents  the  baseband  signature  sequence 

_ , . th 

of  the  k user.  Using  an  appropriate  trigonometric 
identity,  equation  (2.4)  becomes 
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Rjik(i)  = J | I.ft-T^t  + iT-Tjcos  (0j-0k)dt 

+ \ | Sjft  - Xj  )sk(t  + iT  - xk)cos^2wct  + ©j  + 0k)dt 


(2.5) 


For  wct  » 1 , the  second  integral  in  equation  (2.5)  is  zero 
and  the  signature  waveform  cross-correlations  become 


Clearly,  equation  (2.6)  will  be  maximized  for  baseband 
transmissions.  Since  the  performance  of  an  asynchronous  DS- 
CDMA  communication  system  will  in  general  become  better  as 
the  signature  waveform  cross-correlations  become  smaller  in 
magnitude,  baseband  transmissions  will  therefore  represent  a 
worst  case  performance  with  respect  to  carrier  phases.  For 
the  remainder  of  this  chapter,  baseband  polar  signaling  will 
be  the  assumed  modulation  format. 

It  was  shown  by  Verdu  [29]  that  a sufficient  statistic 
for  determining  the  maximum- likelihood  transmitted  sequences 
through  an  AWGN  channel  can  be  obtained  by  passing  the 
entire  received  signal  r(t)  through  a bank  of  filters 
matched  to  the  K signature  waveforms.  Therefore,  the 
proposed  multiuser  receiver  will  consist  of  a matched  filter 
bank  followed  by  an  appropriate  decision  algorithm,  as  shown 
in  Figure  2.3.  The  matched  filter  bank  is  shown  in  more 
detail  in  Figure  2.4,  from  which  it  is  apparent  that  both 
bit  synchronization  and  code  synchronization  are  required  in 
order  to  operate  the  matched  filter  bank  properly.  Bit 
synchronization  is  necessary  in  order  to  sample  the 


(2.6) 
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Figure  2.3  - Structure  of  Multiuser  Detector 


SK  (t-iT-V 

Figure  2.4  - Matched  Filter  Bank  for  Asynchronous  System 

integrators  at  the  proper  intervals,  and  code 
synchronization  is  necessary  in  order  to  align  the  replica 
of  each  user's  signature  waveform  with  the  corresponding 
arriving  waveform.  Acquisition  of  each  user's  time  delay  is 
thus  necessary  to  achieve  code  synchronization.  The  output 
of  the  k*-h  user's  matched  filter  sampled  at  the  end  of  the 
ith  interval  is 
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iT+T+Tk 

yk(i)  = J r(t)sk  (t  - iT  - xk  )dt 

iT+Tk 

K K 

= X RW  (1)VW: (i  “ 1)b3 (i  " 1}  + X R*.j  (0  (i)  (2.7) 

j=k+l  j=l 

k-1 

+ XRk.j  (-1)7wj(i  + 1)bj(i  + 1)  + nk(i) 

j=i 

Defining  y(i)  = [yx(i)  ,y2(i)  ,•••  ,yK(i)]T  to  be  the  vector  of 
matched  filter  outputs  sampled  at  the  end  of  the  it*1  bit 
interval,  b(i)  = [bx(i)  , b2(i)  ,••• , bK(i)]T  to  be  the  vector  of  it*1 

transmitted  bit  polarities,  w(i)  = diag[A/w1(i),  Vw2<i>,-,VwK<i>] 

to  be  the  matrix  of  received  energies  during  the  i^h  bit 
interval,  and  n(i)  = [n1(i)  ,n2(i)  ,--nK(i)]T  to  be  the  vector  of 

matched  filter  noise  outputs  sampled  at  the  end  of  the  it*1 
bit  interval,  equation  (2.7)  may  be  written  more  compactly 
as 

y(i)  = R(— l)W(i  + l)b(i  + 1)  + R(0)w(i)b(i)  + R(l)w(i  - l)b(i  - 1)  + n(i)  (2.8) 

For  a multiuser  receiver  operating  on  a data  window  spanning 
bits  i through  i+L-1,  the  entire  vector  of  matched  filter 
sampled  outputs  for  the  window  can  be  written  as 

yL  = RLKWLbL  + + PWebe  (2.9) 

where  yL  =[y(i)T,  y(i  + l)T,y(i  + 2)T,  ”-,y(i  + L - 1)T]T  , 
bL  = [b(i)T,b(i  + l)T,-”,b(i  + L-l)T]  , WL  = diag[w(i),W(i  + l),--- ,W(i  + L-l)]  , 
nL  =[n(i)T,H(i  + l)T,---,E(i  + L-l)T]T,  w„  =diag[W(i-l)  ,W(i  + L)]  , 
be  = [b(i-l)T,b(i  + L)T]  , and  the  LK  by  LK  and  LK  by  2K  matrices 
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It  is  important  to  note  that  if  the  last  term  in 
equation  (2.9)  was  missing,  then  equation  (2.9)  would  be  the 
same  equation  as  for  the  vector  of  matched  filter  outputs 
for  the  full  length  decorrelator  [38] . The  additional  term 
PWebe  is  needed  in  equation  (2.9)  to  take  into  account  the 
edge  effects  of  the  windowing  operation.  Referring  back  to 
Figure  2.1,  the  matrix  We  contains  the  received  powers  of 
the  bits  immediately  preceding  and  immediately  following  the 
current  window;  likewise,  the  vector  be  contains  the 
polarities  of  these  bits.  The  matrix  P accounts  for  the 
cross-correlations  between  the  signature  waveforms  modulated 
by  these  out-of-window  interfering  bits  and  those  modulated 
by  the  bits  within  the  window.  For  the  full  length 
decorrelator,  the  decorrelating  operation  is  performed  on 
all  M bits  of  all  users  at  once,  and  thus  the  full  length 
decorrelator  may  be  thought  of  as  a windowed  decorrelator 
which  uses  only  one  window  of  length  M bits.  Since  the  K 
users  are  not  transmitting  prior  to  the  first  bit  nor  after 
the  Mth  bit,  there  is  no  source  of  out-of-window  MAI  and 
hence  the  term  PWebe  is  not  needed. 

Performance  Measures 

The  performance  measures  of  interest  are  the 
probability  of  error  and  the  multiuser  asymptotic  efficiency 
of  each  user,  originally  defined  by  Verdu  [29].  For  the  k^ 
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user  operating  at  a probability  of  error  Pe  with  a received 

energy  wk  , the  asymptotic  efficiency  is  defined  as 

woff 

T|k(i)  = lim — *11.  (2.10) 

a->o  Wk 

where  wef£  is  the  received  energy  required  by  a single  user 
to  achieve  the  same  probability  of  error. 

From  equation  (2.10)  it  is  evident  that  an  asymptotic 
efficiency  equal  to  one  would  allow  the  k^h  user  to 
communicate  with  the  same  probability  of  error  as  a single 
user  while  expending  no  additional  transmitter  power.  For 
non-zero  asymptotic  efficiencies,  the  k^h  user  can  attain 
the  same  probability  of  error  as  a single  user,  but  will 
require  additional  transmitter  power.  Of  course,  an 
asymptotic  efficiency  of  zero  would  mean  that  the  k^h  user 
could  never  achieve  the  same  probability  of  error  as  a 
single  user,  regardless  of  how  much  additional  transmitter 
power  was  available.  These  observations  lead  to  the 
conclusion  that  if  the  k^h  user  is  operating  at  a non-zero 
asymptotic  efficiency,  then  in  the  absence  of  noise  the 
probability  of  error  would  be  zero.  In  other  words,  the  k^*1 
user's  performance  is  not  interference  limited.  On  the 
other  hand,  if  the  kfc^  user  is  operating  at  an  asymptotic 
efficiency  of  zero,  then  in  the  absence  of  background  noise 
there  will  still  be  an  irreducible  error  floor  due  to  the 
presence  of  the  other  users.  Thus,  the  kth  user's 
performance  will  be  interference  limited  and  hence  the 
receiver  will  not  be  immune  to  the  near- far  effect. 
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Finally,  just  as  probability  of  error  is  in  general  a 
function  of  signal-to-noise  ratio,  asymptotic  efficiency  is 
in  general  a function  of  the  differences  between  received 
energies.  A useful  performance  measure  therefore  is  the  kt^1 
user's  near-far  resistance,  defined  as  the  minimum  value  of 
the  k^h  user's  asymptotic  efficiency  over  all  possible 
energies  of  the  other  bits,  both  interfering  and 
noninterfering . 


The  standard  decorrelating  strategy  of  multiplying  the 
vector  of  matched  filter  outputs  by  the  inverse  of  the 
signature  sequence  cross-correlation  matrix  will  no  longer 
yield  a decision  statistic  vector  entirely  free  from  MAI  due 
to  the  presence  of  the  out-of -window  MAI,  namely  the  term 

PWebe  in  equation  (2.9).  The  strategy  chosen  here  is  to 
first  subtract  off  an  estimate  of  the  out-of -window  MAI, 
call  it  PWebe , from  yL  before  multiplying  by  R“*  . Since 
the  matrix  P is  known,  this  strategy  amounts  to  obtaining 
reliable  estimates  of  the  received  energies  contained  in  We 
and  the  bit  polarities  contained  in  be . The  final  decision 
statistic  vector  for  the  current  window  will  then  become 


Proposed  Multiuser  Receiver  Structure 


(2.11) 


If  the  receiver  was  able  to  produce  perfect  estimates  or  had 
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knowledge  of  the  true  received  energies  and  bit  polarities, 
then  the  last  term  in  equation  (2.11)  would  equal  zero  and 
equation  (2.11)  would  be  the  same  as  given  by  Lupas  [38]  for 
a full  length  decorrelator  for  which  all  users  transmit  L 
bits.  In  practice,  of  course,  the  receiver  will  not  be  able 
to  produce  perfect  estimates;  however,  from  equation  (2.11) 
it  is  reasonable  to  expect  that  the  windowed  decorrelator’s 

performance  will  improve  as  the  estimates  We  and  be  become 
closer  to  the  true  quantities  We  and  be . A block  diagram 
depicting  the  general  structure  of  the  windowed  decorrelator 
is  presented  in  Figure  2.5.  The  delay  block  inserted  before 
the  matched  filter  bank  is  included  since  knowledge  of  the 
timing  information  will  be  necessary  in  order  to  operate  the 
matched  filter  bank.  The  problem  of  extracting  this  timing 
information  from  the  received  signal  is  currently  an  area  of 
active  research  [57] -[59].  In  this  work  it  will  be  assumed 
that  perfect  timing  information  is  available  to  the  receiver 
and  that  the  matched  filter  bank  is  ideal  and  operating 
according  to  equation  (2.9). 

The  requirement  of  producing  an  estimate  of  the  out-of- 
window MAI  affecting  the  current  window  can  be  divided  into 
two  separate  tasks:  estimating  the  received  energies  of  the 

interfering  bits  and  estimating  the  polarities  of  the 
interfering  bits.  For  an  RF  system,  the  received  carrier 
phase  of  the  interfering  bits  will  be  required  as  well; 
however,  since  it  is  assumed  that  the  multiuser  receiver  is 
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Figure  2.5  - Generalized  Block  Diagram  of  Windowed 

Decorrelator 

tracking  the  carrier  phases  for  the  purpose  of  coherent 
demodulation,  these  quantities  would  not  need  to  be 
estimated.  An  interesting  point  to  note  is  that  if  a 
noncoherent  modulation/demodulation  technique  were  employed 
by  the  system,  such  as  DPSK,  the  carrier  phases  of  the  out- 
of -window  MAI  bits  would  need  to  be  known  and  hence  would 
require  estimating.  Thus  the  windowed  decorrelator  would  in 
a sense  be  self-defeating,  since  the  advantage  of  using  a 
noncoherent  scheme  is  that  knowledge  of  carrier  phases  is 
not  required.  Nonetheless,  for  situations  in  which  tracking 
the  received  carrier  phases  is  either  not  possible  or  not 
desirable,  a noncoherent  technique  could  be  used  with  a 
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windowed  decorrelator  provided  that  estimates  of  the  carrier 
phases  before  and  after  each  window  could  be  obtained. 

To  eliminate  the  need  for  estimating  the  polarities  of 
the  out-of -window  MAI  bits,  pilot-symbol-assisted-modulation 
(PSAM)  can  be  used.  In  PSAM,  known  symbols  (pilot  symbols) 
are  periodically  inserted  into  the  transmitted  sequence. 

PSAM  is  commonly  used  to  obtain  information  on  the  channel 
characteristics  in  a dynamically  evolving  environment  [28] . 
For  the  multiuser  scenario  considered  here,  the  (L+l)th  bit 
of  each  user  would  be  a pilot  bit  as  illustrated  in  Figure 
2.6.  The  pilot  bits  thus  serve  to  separate  the  data  windows 
from  one  another  and  the  multiuser  receiver  will  simply 
discard  them  at  the  output  of  the  matched  filter  bank. 


L data  bits 

pilot 

bit 

L data  bits 

pilot 

bit 

Figure  2.6  - Data  Block  Structure  of  Each  User 

Incorporating  PSAM  into  the  system  offers  two  important 
advantages  over  other  polarity  estimation  techniques. 

First,  since  the  polarities  of  the  pilot  bits  are  known 
apriori  by  the  multiuser  receiver,  PSAM  will  result  in 
computational  savings  compared  to  any  method  which  actually 
computes  out-of-window  MAI  polarity  estimates.  The  second 
advantage  derived  from  using  PSAM  is  that  be  = be  for  all 
windows,  whereas  any  algorithm  which  attempts  to  compute  be 
will  always  be  subject  to  errors  introduced  by  a combination 
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of  the  background  noise  and  MAI.  That  is,  just  as  in 
decision-directed  feedback,  it  is  possible  that  a polarity 
estimate  would  be  in  error  and  could  therefore  cause  errors 
in  the  final  polarity  decisions  of  the  data  bits  within  the 
current  window.  If  these  data  bits  are  then  used  by  the 
algorithm  to  estimate  the  out-of-window  MAI  polarities 
affecting  the  next  window,  the  end  result  could  be  a 
"snowball  effect"  in  which  one  estimation  error  leads  to 
many  final  bit  decision  errors.  Of  course,  using  PSAM 
eliminates  this  potential  threat. 

There  are,  however,  two  disadvantages  in  using  PSAM. 
First,  each  user  will  have  a overhead  due  to  the 

presence  of  the  pilot  bits.  The  overhead  results  in  a 
reduction  of  each  user's  data  rate  by  1^%  and  an  increase 
in  total  required  transmitter  power  by  The  reduction 

in  data  rate  and  increase  in  transmitter  power  occur  since 
each  pilot  bit  has  no  message  content;  the  same  information 
previously  transmitted  in  L bits  will  now  require  L+l  bits. 
Since  the  mobile  units  will  most  likely  be  battery  powered, 
any  increase  in  the  total  amount  of  transmitter  power 
required  can  be  costly.  For  this  reason,  care  should  be 
used  when  incorporating  PSAM  to  keep  the  overhead  as  low  as 
possible.  The  second  disadvantage  in  using  PSAM  is  that  the 
multiuser  receiver  must  now  acquire  and  maintain  pilot 
synchronization.  In  other  words,  the  link  from  the  mobile 
units  to  the  base  station  is  no  longer  asynchronous  in  the 
strictest  sense  since  the  receiver  must  align  the  m^h  pilot 
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bit  of  each  user  at  the  output  of  the  matched  filter  bank  to 
within  T seconds.  From  the  viewpoint  of  each  user,  though, 
the  system  is  still  completely  asynchronous  since  the  K 
users  do  not  need  to  coordinate  their  transmitted  sequences 
in  time;  this  task  is  relegated  to  the  base  station.  Recall 
that  both  code  synchronization  and  bit  synchronization  are 
required  to  properly  operate  the  matched  filter  bank  (Figure 
2.4).  Assuming  that  the  receiver  has  already  obtained  these 
levels  of  synchronization,  obtaining  pilot  synchronization 
would  then  seem  quite  feasible.  For  the  remainder  of  this 
analysis  it  will  be  assumed  that  pilot  synchronization  has 
been  accomplished. 

All  that  remains  is  for  the  multiuser  receiver  to  form 
the  matrix  of  required  energy  estimates,  We  . For  a full 
length  decorrelator,  the  magnitudes  obtained  from  the 
decorrelating  solution  are  maximum- likelihood  estimates  of 
the  corresponding  received  energies  [29].  That  is,  the 
decorrelator  bases  its  final  polarity  decisions  upon  the 
sign  of  the  appropriate  element  in  the  decision  statistic 
vector,  but  the  magnitude  of  that  element  yields  a maximum- 
likelihood  received  energy  estimate.  Of  course,  a full 
length  decorrelating  receiver  does  not  require  any  knowledge 
of  received  energies;  however,  this  available  side 
information  can  be  used  to  provide  the  required  energy 
estimates  for  the  windowed  decorrelator.  For  slowly  varying 
received  energies,  the  multiuser  receiver  can  use  the  energy 
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estimates  derived  from  the  previous  window's  decorrelating 
solution  as  its  out-of -window  MAI  energy  estimates  for  the 
current  window.  The  energy  estimates  used  for  the  very 
first  window  will  be  zero,  since  no  past  decorrelating 
solutions  are  available.  Due  to  the  presence  of  the  out-of- 
window MAI,  this  estimation  scheme  will  no  longer  be 
maximum- likelihood  but  its  savings  in  computational 
complexity  make  it  a very  attractive  choice. 

The  proposed  receiver  structure  is  now  complete.  A 
detailed  block  diagram  of  the  receiver  is  shown  in  Figure 
2.7.  It  is  important  to  note  that  the  cross-correlation 
matrix  RLK  is  a function  of  only  the  active  signature 
sequences  and  the  relative  time  delays  between  the  active 
users.  As  these  quantities  may  remain  static  for  many 
consecutive  windows,  the  computation  of  R~*  will  only  need 
to  be  performed  when  the  receiver  senses  a change  in  either 
the  number  of  active  users  or  their  timing  configuration. 

Performance  Analysis 

The  probability  that  the  ith  bit  of  the  kth  user  is  in 
error  will  be 

Pk(i)  = Pr{sign[zL]“  = 1 1 bk(i)  = -l}  (2.12) 

where  the  notation  [x]“  represents  the  vt*1  element  in  the 
uth  window  of  the  vector  x . The  window  index  u is  given  by 

u = int[^]  + l (2.13) 

and  the  vector  element  v is  given  by 
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Figure  2.7  - Detailed  Block  Diagram  of  Windowed  Decorrelator 

v = (i  - 1)K  + k - intp^jLK  (2.14) 

where  the  notation  int[x]  represents  the  integer  part  of  x. 
Now, 

Pr{M  >0|bk(i)  = -l}=  Pr{  [wLbLJv+[R^HL]; 


R«p(webe-WebeJ  > 0 1 bk(i)  = — 1 } 


(2.15) 


which  reduces  to 


Pk(i)  Pr|[RLKnL]v  > -\/wk(i)  j^RLKP^Webe  WebeJ 


(2.16) 


Let  cm,n  denote  the  element  in  the  m^h  row  and  n^^1  column  of 


the  matrix  R * . Then 
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Rlk  = C = 


Cl,l  Cl,2 


'1,LK 


C2,l  C2,2  C2,LK 


CLK,1  CLK,2  ‘*’CLK,LK 


(2.17) 


The  noise  component  in  equation  (2.16)  will  be  Gaussian  with 
zero  mean  and  a variance  which  is  the  diagonal  element 

of 

e{[r'^l  ]u  ([r~A  ]u  )' } = (a  2Rlk  )r 


= °2r^ 


(2.18) 


Thus 


var 


{[RLKnL];}=o2cv,v 


(2.19) 


From  equation  (2.16),  the  probability  of  error  will  be 


Pk(i)  = Q 


V^i)-[R-p(webe-weLe) 


(2.20) 


where  the  Q-function  is  defined  as 

Q(x)  = ^Je^dt 


(2.21) 


Note  that  if  PSAM  was  not  used,  then  the  probability  of 
error  given  by  equation  (2.20)  would  need  to  be  averaged 

2k  

over  the  2 possible  MAI  polarity  vectors  be  and  would  thus 
be  a sum  of  Q- functions. 

Before  proceeding  further,  it  is  of  interest  to 
consider  the  ideal  case  in  which  the  receiver  has  perfect 
knowledge  (or  perfect  estimates)  of  the  out-of-window  MAI. 
This  special  case  will  indicate  the  best  performance 
possible  from  a windowed  decorrelating  receiver.  In  this 
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scenario  Webe=Webe  and  equation  (2.20)  reduces  to 


f 

Pk(i) = Q 


Vw^iy" 

^VCv-v  J 


(2.22) 


which  is  a function  of  only  the  k^h  user's  received  energy 
and  not  of  the  other  users'  energies.  Since  the  probability 
of  error  for  a single  user  employing  BPSK  in  an  AWGN  channel 
is 


( 


Pk(i)  = Q 


■i 


w 


eff 


(2.23) 


V ° 

then  clearly  weff  = wk(i)/cv v and  from  equation  (2.10)  the 

asymptotic  efficiency  for  the  ith  bit  of  the  kth  user  will 
be 


Tlk(i) 


(2.24) 


Since  the  asymptotic  efficiency  is  a nonzero  constant,  the 
near-far  resistance  will  be  equal  to  this  constant  as  well 
and  the  receiver  is  indeed  immune  to  the  near- far  effect. 

The  actual  value  of  the  near- far  resistance  will  depend  upon 
the  cross-correlations  between  the  signature  sequences 
through  the  term  cvv . It  should  be  noted  that  equation 

(2.24)  is  the  same  result  which  would  be  obtained  for  a full 
length  decorrelator  operating  on  a transmission  of  L bits. 
Another  interesting  observation  is  that  both  the  probability 
of  error  and  the  asymptotic  efficiency  of  the  windowed 
decorrelator  are  a function  of  the  transmitted  bit  index; 
more  specifically,  they  are  a function  of  the  bit  position 
within  each  window  through  the  term  c 


Thus  the  kth 
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user's  probability  of  error  and  asymptotic  efficiency  will 
differ  from  bit  to  bit  within  each  window,  a property  which 
is  different  from  that  of  an  infinite  length  decorrelator 
whose  performance  measures  are  independent  of  the  bit  index. 

To  evaluate  the  probability  of  error  given  by  equation 
(2.20)  further,  the  last  term  in  the  numerator  must  be 
investigated.  For  convenience  it  will  be  assumed  that  all 
pilot  bits  have  a polarity  of  +1,  although  in  practice  this 
is  not  a necessary  condition.  The  uth  window  of  the  vector 


[w,b  — W b 1 

lee  e e 


will 


[webe-webe] 


then  become 

Vw  i ((u  - 1)(L  + 1))  - 7wx  ((u  - 1)(L  + 1)) 

Vw  2 ((u  “ !)(L  + !))  “ V^2  ((u  ~ !)(L  + !)) 

Vwk((u  - 1)(L  + 1))  - ^/wK  ((u  - 1)(L  + 1)) 
Vw,  (u(L  + 1))  - A/w1  (u(L  + 1)) 

Vw2  (u(L  + !))  - Vw2  (u(L  + 1)) 


Vwk  (U(L  + !))  - V^k  (U(L  + !)) 


(2.25) 


Define  the  error  in  the  energy  estimate  for  the  jt*1  bit  of 
the  kth  user  as 

Awk(j)  = Vwk(3)  ~ (2.26) 


Actually,  equation  (2.26)  represents  an  error  in  an 
amplitude  estimate;  indeed,  the  receiver  actually  needs  to 
form  amplitude  estimates  rather  than  energy  estimates.  For 
consistency,  however,  equation  (2.26)  will  be  referred  to  as 
an  energy  estimation  error.  Expanding  out  the  term 

[R-^p(webe-Webe)]"  yields 
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[R^p(webe-webe)]u  = X X°  vdRU  (DAw,  ((u  - D(L  + D) 

v 1=2  j=l 

K-l  LK 


(2.27) 


+ X Xcv,ri,-(l-dk  (“DAW!  (u(L  + 1)) 

1=1  j=(L-l)K+l+l 

Recall  that  the  subscript  v denotes  the  vt^1  element  in  the 
corresponding  vector,  and  thus  in  equation  (2.27)  v will 
range  from  1 to  LK.  Substituting  equation  (2.27)  into  the 
equation  for  probability  of  error  given  by  equation  (2.20) 
produces 


P*  (i)  = Q - X X c vjRu  d^i  ((u  - !)(L  + !)) 

I 1=2  j=l 

K-l  LK 

“ X Xcv.jRj-(L-DK,i  (-DAwj.  (u(L  + 1))  j- 

1=1  j=(L— 1JK+1+1  <7 

Note  that  Awk(i)  may  be  either  positive  or  negative, 

depending  upon  whether  the  corresponding  received  energy  was 
underestimated  or  overestimated,  respectively.  It  is 
therefore  necessary  to  consider  the  worst-case  probability 
of  error  which  corresponds  to  the  following  least  favorable 
event 

s ignjAw,  ((u  - 1)(L  + 1))}  = s ign{c  vlR (1)},  1 = 2,3 , • • ■ , K 

and  (2.29) 

s ignfAw,  (u(L  + 1))}  = s ign{c  (-1)},  1 = 1,2,  • • • , K - 1 

The  worse-case  probability  of  error  will  therefore  be 

f K 1-1 


(2.28) 


-Jcv.v 


P,(iLx  = Q W^-XXFvjRj. iWlKCCu-iXL  + Dj 

l 1=2  j=l 

K-l  LK  - , 

“X  X |c  v,j  R j-(L-i)K,l  (-1)||  Aw  1 (u(L  + l))| }— 7=  i 

1=1  j=(L-l)K+l+l  <7VCv,v  J 

The  actual  value  of  Pk(i)max  will  depend  upon  how  large  in 


(2.30) 


magnitude  the  out-of -window  MAI  energy  estimation  errors 
are.  Of  course,  for  perfect  energy  estimates  equation 
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(2.30)  reduces  to  equation  (2.22).  The  worst-case 
asymptotic  efficiency  can  be  found  from  equation  (2.30); 
applying  equation  (2.10)  gives 


r\^iKin=maX'4 


r r JLfcJ  |Aw,  ((u  - 1)(L  + 1))| 

{ 0.  1-IE|cv,1R11(1)|I  * - 

1 1=2  j=l  VW 

AWj(u(L  + l))|  , i i 
Vw~(i)  ■* 


K— 1 LK 


(2.31) 


-I  £ Kjr 

1=1  j=(L-l)K+l+l 


j-(L-l)K 


The  form  of  equation  (2.31)  provides  some  valuable  insight 
into  the  behavior  of  the  windowed  decorrelator's  asymptotic 
efficiency.  First,  as  the  magnitudes  of  the  out-of-window 
MAI  energy  estimation  errors  approach  zero,  the  asymptotic 
efficiency  approaches  the  constant  value  found  in  equation 
(2.22).  As  the  energy  estimation  errors  become  larger  in 
magnitude,  the  asymptotic  efficiency  will  decrease  until  the 
two  summation  terms  yield  a combined  value  of  one,  after 
which  the  asymptotic  efficiency  will  remain  zero.  Thus  the 
near-far  resistance  of  the  windowed  decorrelator  is  zero. 

It  is  important  to  observe  from  equation  (2.31),  however, 
that  the  errors  in  the  estimates  of  the  interfering 
amplitudes  are  divided  by  the  amplitude  of  the  desired  user 
in  determining  the  desired  user's  asymptotic  efficiency. 
Thus,  the  estimation  errors  are  in  effect  normalized  by  the 
desired  user's  amplitude.  The  stronger  the  desired  user's 
received  signal  is,  the  less  effect  the  energy  estimation 
errors  will  have  upon  the  desired  user's  performance.  Thus, 
even  though  the  near-far  resistance  is  zero,  the  differences 
between  received  energies  at  which  the  asymptotic  efficiency 
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becomes  zero  may  be  sufficiently  large  to  cover  any 
reasonable  near- far  scenario  which  might  be  encountered  in 
practice . 

All  that  remains  to  complete  the  performance  analysis 
of  the  windowed  decorrelator  is  to  substitute  appropriate 


expressions  for  Awd ((u - 1)(L  + 1))|  and  AWj (u(L  + 1))|  into  the 

expression  for  worst-case  asymptotic  efficiency  given  by 
equation  (2.31).  Decorrelating  the  first  window  will  yield 
energy  estimates  of  the  K users  for  the  bit  intervals  1 
through  L.  The  second  window  will  use  the  estimates 
corresponding  to  the  most  recent  bit  interval,  the  Lth 
interval,  to  modify  its  matched  filter  output  vector. 
Decorrelating  the  second  window  will  produce  energy 
estimates  of  the  K users  for  the  bit  intervals  L+2  through 
2L+2 . The  third  window  will  then  use  the  estimates 
corresponding  to  the  (2L+2)nd  interval  to  modify  its  matched 
filter  output  vector,  and  so  on. 

Consider  the  very  first  window  (u=l)  spanning  bits  1 
through  L.  No  past  estimates  are  available  and  so 

w.(L  + l)  = 0,  j = 1,2,"-, K (2.32) 


The  receiver  forms  the  decision  statistic  vector 

= WLbL+yL+R^PWebe 


(2.33) 


where  the  noise  component  vector  yL  has  been  defined  as 

Yi,=R«nL  (2.34) 

The  last  term  in  equation  (2.33)  represents  the  out-of- 
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window  MAI  corrupting  the  first  window's  decision  statistic 
vector.  Expanding  out  this  term  and  denoting  it  by  1(1) 
yields  its  v*-*1  element  to  be 

K— 1 LK 

[I(1)]v  = I Scv,jRj— , l-i,k.i(-1)VW3(L  + 1)  (2.35) 

1=1  j={L— DK+l+l 

The  estimate  of  the  kth  user's  received  energy  will  be 
formed  as 

>/wk(L)  =|A/wk(L)bk(L)  + Yk(L)  + Ilt(l)|  (2.36) 

Since  the  expression  for  asymptotic  efficiency  given  by 
equation  (2.31)  was  derived  in  the  limit  as  a — » 0 , the  noise 

component  in  equation  (2.36)  may  be  neglected.  Thus 

7wk(L)  =|Vwk(L)bk(L)  + Ik(l)|  (2.37) 

The  error  in  this  estimate  is  then 


Awk(L)  = A/wk(L)  — |A/wlc(L)bk(L)  + Ik(l)|  (2.38) 

Equation  (2.37)  is  of  the  general  form 

Vw  - |±Vw  ± l|  (2.39) 

which  is  shown  in  Appendix  A to  be  bounded  by 

— |l|  < Vw  - |±Vw  ± l|  < |l|  (2.40) 

Therefore  the  estimation  error  will  be  bounded  by 

|Awk(L)|<|lk(l)|  (2.41) 


The  worst-case  estimation  error  will  occur  when  each  term 
inside  the  summations  in  equation  (2.35)  has  the  same  sign, 
resulting  in 


K— 1 LK 


|Awk(L)|max  ^|C(L-l)K+l,j^:HL-l)K,V  1)  VW 1^  + (2.42) 


1=1  j=(L-l)K+l+l 


To  proceed,  the  vector  of  maximum  possible  energy  estimation 
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errors  obtained  from  decorrelating  the  uth  window  is  defined 
as 


Aw(u)  = 


I Aw,  (u(L  + 1)  - l)|roax 
|Aw2(u(L  + l)-lLx 


|Awk  (u(L  + 1)  - l)([oax 


(2.43) 


Note  that  the  k^-h  element  of  Aw(l)  will  be  given  by  equation 


(2.42) . 


Now  consider  the  second  window  (u=2),  spanning  bits  L+2 
through  2L+2 . Recall  that  the  (L+l)th  bit  of  each  user  is  a 
pilot  bit  and  is  discarded  by  the  multiuser  receiver.  The 
decision  statistic  vector  for  the  second  window  will  be 

zL=WLbL+yL  + R^p(webe-Webe)  (2.44) 


In  an  AWGN  channel,  the  received  energies  will  be  changing 
slowly  relative  to  a bit  duration.  For  small  window 
lengths,  it  will  be  assumed  that  the  received  energies  are 
constant  over  a span  of  L+2  bits;  in  other  words 


^w.(i  + m)  = ^w.(i),  |m|  < L + 2 


(2.45) 


With  the  above  simplifying  assumption,  equation  (2.25) 


becomes 


Aw,(L) 

Aw2(L) 


Awk(L) 

Aw,(L) 

Aw2(L) 


Awk(L) 


(2.46) 
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Pre -multiplying  equation  (2.46)  by  R^P  and  using  the  worst 
case  estimation  errors  given  by  equation  (2.42)  results  in 
an  expression  for  the  vth  element  of  the  out-of-window  MAI 
term  affecting  the  second  window's  decision  statistic 
[ra]v  =[RiP(w,b,-W.b.)£ 

K-l  LK 

1(l)[Awi(L)]max  XCV.jRj-(L-l)K,l 

1=2  j=l  1=1  j=(L-l)K+l+l 

By  re-grouping  terms,  equation  (2.47)  can  be  written  as 


itAwittOL, 


K 1-1 

=11'./) 


1-1 


K-l 


[l(2)]v  =E[Aw1<L)USeVJR,I(l|+2[AwI(L)]. 


LK 


1=2 


/Cv.jRj-(L-l)K,l^  ^ 


3=1 


1=1 


j=(L— 1)K+1+1 


(2.48) 


Finally,  the  energy  estimates  will  be  based  upon  the  last 
bit  of  each  user  within  the  window  so  that  the  element 

in  the  vector  of  worst  case  energy  estimation  errors 
obtained  from  decorrelating  the  second  window  becomes 


[[Aw(2)Ljk  2 [[Aw(1)Lx  1 XC(L— DK+k.jR  j.i  d) 

+ S[[^W^^]max]  SC(L-l)K+k,jRj-(L-l)K,l(”^) 

1=1  1 j=(L-l)K+l+l 


(2.49) 


A relationship  between  the  vector  of  worst  case  energy 
estimation  errors  obtained  from  consecutive  windows  is  now 
evident  from  equation  (2.49).  The  initial  out-of-window  MAI 
term  corrupting  the  first  window's  decision  statistic  vector 
has  been  reduced  after  the  second  window  has  been 
decorrelated.  This  process  will  repeat  itself  for  each 
additional  window  which  is  processed  by  the  receiver;  in 
other  words,  the  initial  out-of -window  MAI  term  will 
decrease  as  the  number  of  windows  processed  by  the  receiver 
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increases.  The  generic  form  of  this  relationship  will  be 

Aw(u)  = AAw(u  - 1),  u>2  (2.50) 

where  A is  the  appropriate  recursion  matrix.  To  determine 
an  expression  for  the  contents  of  the  K by  K matrix  A, 
equation  (2.50)  can  be  expanded  out  to  obtain 


Awx  (u(L  + 1)  - 1) 
Aw2  (u(L  + 1)  — l) 

_AwK(u(L  + l)-l) 


ai.l  ai,2 
a2,l  a2,2 


a K,1  aK,2 


1,K 


2,K 


lK,K 


Awx  ((u  - 1)(L  + 1)  - 1) 
Aw  2 ((u  — 1)(L  + 1)  — l) 

_AwK((u-l)(L + !)-!)_ 


(2.51) 


It  is  clear  from  equation  (2.51)  that  the  pth  row  of  the 
matrix  A must  correspond  to  the  pt*1  user,  so  the  subscript  k 
in  equation  (2.49)  can  be  replaced  with  p.  Also,  the  qth 
column  of  the  matrix  A must  correspond  to  the  1th  element  in 
the  previous  error  vector,  so  the  subscript  1 in  equation 
(2.49)  can  be  replaced  with  q.  Finally,  since  the  ranges  of 
1 in  the  summations  in  equation  (2.49)  overlap,  each 
summation  will  contribute  a component  to  each  apq  term  and 

thus  their  contributions  must  be  summed.  The  resulting 
coefficient  of  [Awfu-l)^  in  equation  (2.49)  (for  which  u=2) 
must  then  be  apq.  Therefore,  the  element  in  the  pth  row 
and  qth  column  of  the  matrix  A is  given  by 


ap.q  = 


I 

j=(L— l)K+l+q 


q-i 

C(L-l)K+p,j  (L— l)K,q(_  + ^ |C 

j=l 


(L-l)K+p,j  ^j,q^ 


(2.52) 


From  equation  (2.50)  the  vector  of  worst  case  energy 
estimation  errors  obtained  from  decorrelating  the  uth  window 
can  be  related  to  those  obtained  from  decorrelating  the 
first  window  by 
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Aw(u)  = Au  Aw(l),  u > 1 


(2.53) 


Under  the  assumption  of  equation  (2.45),  equation  (2.46)  can 
be  substituted  into  the  expression  for  asymptotic  efficiency 
given  in  equation  (2.31)  to  obtain 

, r r v’vu  ,lAwj((u-1)L)| 

•{  0.  l-EEK^.dlMr  1 

L 1=2  j=i  v^M1) 

|AWj((u-l)L)| 


=max2 


K-l  LK 


'jLl  ^ lCv,j^j-(L-l)K,l^  ^1" 
1=1  j=(L— DK+l+l 


(2. 54) 

Vwk<i)  ^7^77  } 

Finally,  using  the  relationship  between  the  u*-*1  window  and 
the  first  window  given  by  equation  (2.53),  the  windowed 
decorrelator's  worst  case  asymptotic  efficiency  for  the  ith 
bit  of  the  user  is 


K(i)L„ =max2  { °'  -!==[  1-IIK,r,iH 

V^v'v 


1=2  j=l 


[au-2Aw(1)]i 

7wk<i) 


K-l 


_y  yu  R i pi  Aw(1)]i 


(2.55) 


]} 


1=1  j=(L— 1JK+1+1 

where  as  before  the  notation  [x]x  represents  the  1th  element 

in  the  vector  x and  the  window  index  u is  related  to  the  bit 
index  i through  equation  (2.13). 

A very  important  observation  from  equation  (2.55)  is 
that  if 

(2.56) 


lim Au =0 

U— >oo 


then 


limr|k(i)  = - 


(2.57) 


Thus,  if  equation  (2.56)  is  satisfied,  then  the  asymptotic 
efficiency  of  the  windowed  decorrelating  receiver  will 
converge  to  the  same  value  as  for  the  full  length 
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decorrelating  receiver.  In  other  words,  after  a sufficient 
number  of  windows  has  been  processed,  the  windowed 
decorrelator  will  for  all  practical  applications  be  near-far 
resistant  and  hence  immune  to  the  near-far  effect.  It  is 
therefore  important  to  determine  under  what  conditions 
equation  (2.56)  will  hold;  however,  it  is  easier  to  first 
consider  the  conditions  under  which  equation  (2.56)  will  not 
hold. 

It  was  shown  by  Lupas  [38]  that  a stable  realization  of 
the  decorrelating  receiver  will  exist  as  the  message  length 
approaches  infinity  if  and  only  if  the  signature  waveform 
cross-correlations  satisfy 

det[RT(l)e:“  + R(0)+R(l)e~ja>]*0,  Vcoe[0,2rc)  (2.58) 

One  way  to  arrive  at  this  result  is  by  considering  equation 
(2.9),  without  the  term  due  to  the  edge  effects  of  the 
windowing  operation,  in  the  limit  as  the  message  length 
approaches  infinity  and  taking  Z-transforms . Defining 

S(z)  = Rt(1)z  + R(0)+R(1)z-1  (2.59) 

it  is  clear  that  the  transfer  function  of  the  infinite 
length  decorrelating  receiver  must  be 

G(z)  = [S(z)]  1 (2.60) 

A stable  transfer  function  will  therefore  not  exist  unless 
equation  (2.58)  is  satisfied.  Clearly,  if  equation  (2.58) 
is  not  satisfied  then  equation  (2.56)  cannot  be  satisfied; 
that  is,  if  a stable  infinite  length  decorrelator  does  not 
exist,  then  there  will  be  no  receiver  for  the  windowed 
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decorrelator  to  converge  to  as  the  window  index  increases. 
For  a K user  system,  the  appropriate  cross-correlations  will 
need  to  be  checked  according  to  equation  (2.58)  to  determine 
if  an  infinite  length  decorrelator  exists;  however,  for  a 
two  user  system  equation  (2.58)  reduces  to  a rather  simple 
condition.  For  this  special  case  the  cross-correlation 
matrices  may  be  defined  as 


R(0)  = 


_ri2 


P 12 
1 


and  R(l)  = 


0 P 21  * 

0 0 


(2.61) 


so  that  equation  (2.58)  becomes 


det 


Pl2  + P2ie 


+j© 


*0 


l_P  12  + P2ie 

which  upon  taking  the  determinant  yields 


P 12  P21  ^P 12 P21  c°s(CD)  ^ 1 


(2.62) 


(2.63) 


Now, 

P12  P21  2p 12P 21  cos(co)  < 0Pi2 1 |p2i|)  ( 2 . 64 ) 

so  that  if 

|Pi2|  + |p2i|<1  (2.65) 

then 

P12 +P2i +2p12P21cos(C0)<l  (2.66) 

and  equation  (2.58)  is  always  satisfied  and  a stable 
infinite  length  decorrelator  will  exist.  Note  that  equation 
(2.65)  is  a rather  loose  constraint  upon  the  signature 
sequence  cross-correlations;  in  practice  it  is  not  difficult 
to  find  a pair  of  signature  sequences  satisfying  this 
condition. 
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Unfortunately,  the  existence  of  an  infinite  length 
decorrelator  would  not  seem  to  guarantee  that  the  windowed 
decorrelator  will  converge  to  the  same  performance.  It 
therefore  remains  to  determine  under  what  conditions 
equation  (2.56)  will  be  valid.  To  do  this,  some  well  known 
theorems  from  linear  algebra  will  be  used. 

If  the  matrix  A is  diagonalizable  then  it  can  be 
written  as  [60] 

A = BAB_1  (2.67) 


where  A is  the  diagonal  matrix  whose  diagonal  elements 

are  the  eigenvalues  of  A,  and  B is  a matrix  of 
corresponding  eigenvectors.  It  can  be  shown  [60]  that 

Au=BAuB_1  (2.68) 


so  that  equations  (2.56)  and  (2.57)  become 


limrik(i)  = if  limAu=0 


(2.69) 


Since  A is  a diagonal  matrix,  the  above  condition 
simplifies  to 


limri)c(i)  = if  \ <1,  j=l,2,-”,K 


(2.70) 


Thus,  if  the  matrix  A defined  by  equation  (2.52)  is 
diagonalizable  and  its  eigenvalues  all  have  magnitudes  less 
than  one,  then  the  asymptotic  efficiency  of  the  windowed 
decorrelator  will  approach  a nonzero  constant  as  the  window 
index  increases.  Due  to  the  analytical  intractability  of 
determining  the  eigenvalues  of  a K by  K matrix,  the  focus 
here  will  be  on  a two  user  system.  In  addition  to  being 
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mathematically  convenient,  the  two  user  case  may  provide 
valuable  insight  into  the  general  K user  scenario. 

It  is  shown  in  Appendix  B that  the  matrix  A will  always 
be  diagonalizable  for  the  case  of  two  users;  it  therefore 
only  remains  to  compare  the  magnitudes  of  the  eigenvalues  of 
A with  unity.  For  a window  of  length  three  bits,  the 
eigenvalues  of  A are  the  solutions  to  the  equation 

^ — dC  5.6p2X  I 1^- 6.1P2I  l)^1  p21  d<-'5,6<“6.l|  |C6,6C5,l|)  = ^ (2.71) 


Upon  substituting  equations  (B.13)  and  (B.17)  into  equation 
(2.71)  and  using  the  relationships  given  by  equations  (B.4) 
and  (B.5),  the  eigenvalues  of  the  matrix  A are  found  from 
the  quadratic  formula  to  be 


1,2 


_ 1 

P12P21 

+ 

P21N 

+ 

7 

P12P21 

+ 

P21N  ) 

2 

D 

D 

_v 

D 

D J 

f 3 3 

P21N 

P12P21 

P21 

l D 

Jt 

D 

(2.72) 


where 

N = “Pi,  + 2Pi2  + P12P21  ~ P12  (2.73) 

and 

D = PL  - 3pJ2  - P21  + 3p212  + 2p21  - 2p22p21  - 1 (2.74) 

Unfortunately,  equation  (2.72)  does  not  appear  to  lend 
itself  to  a more  compact  form.  In  an  attempt  to  gain 
further  insight  into  the  behavior  of  equation  (2.72)  a 
numerical  approach  was  taken.  The  values  of  the  cross- 
correlations p12  and  p21  were  varied  between  -1  and  +1  in 
increments  of  0.1.  For  each  combination  of  values,  the 
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matrices  R,  C,  and  A were  computed  and  the  eigenvalue  of  A 
with  the  largest  magnitude  was  then  found.  The  results  are 
summarized  in  Table  2.1,  from  which  it  is  readily  apparent 
that  the  condition  for  convergence  of  the  windowed 
decorrelator  given  in  equation  (2.70)  is  equivalent  to 

limT|k(i)  = — if  |p12|  + |p  |<1,  k = 1,2  (2.75) 

i"~  Cv,v 

Interestingly,  the  condition  given  in  equation  (2.75)  is 
equivalent  to  that  given  in  equation  (2.65)  for  the 
existence  of  an  infinite  length  decorrelator.  Thus,  if  a 
stable  infinite  length  decorrelating  receiver  exists  for  a 
two  user  system,  then  the  windowed  decorrelating  receiver 
will  always  become  near-far  resistant  as  the  window  index 
increases . 

Repeating  the  previous  analysis  for  the  more  general  K 
user  system  is  not  feasible.  One  reason  for  the  difficulty 
in  such  an  analysis  is  that  for  a K user  system,  there  will 
be  K(K-l)  cross-correlations  to  deal  with;  finding  a 


Table  2.1  - Dominant  Eigenvalue  of  the  Matrix  A For  a Two 

User  System 


|Pi,l 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

o 

00 

0.9 

1.0 

|P  21 1 

0.1 

.010 

.021 

.036 

.055 

.082 

.123 

.195 

.349 

1.00 

1.18 

0.2 

.021 

.045 

.075 

.116 

.176 

.272 

.457 

1.00 

2.12 

1.45 

0.3 

.033 

.072 

.122 

.192 

.301 

.497 

1.00 

35.5 

1.14 

1.89 

0.4 

.049 

.106 

.183 

.300 

.505 

1.00 

9.01 

1.25 

1.41 

2.74 

IT) 

O 

.069 

.153 

.276 

.489 

1.00 

8.00 

1.26 

1.33 

2.05 

5.03 

0.6 

.098 

.228 

.450 

1.00 

16.4 

1.20 

1.36 

2.00 

3.69 

32.5 

0.7 

.148 

.381 

1.00 

9.84 

1.14 

1.50 

2.22 

3.97 

16.4 

7.02 

CO 

o 

.261 

1.00 

2.07 

1.28 

1.88 

3.00 

6.09 

130 

6.69 

3.05 

0.9 

1.00 

1.13 

1.95 

3.43 

7.56 

975 

8.44 

4.20 

2.70 

1.84 

1.0 

19.9 

9.89 

6.51 

4.79 

3.73 

3.00 

2.44 

2.00 

1.59 

1.00 
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constraint  upon  these  cross-correlations  similar  to  equation 
(2.75)  will  simply  not  be  practical  as  K increases.  Of 
course,  equation  (2.58)  must  be  satisfied  in  order  for  a 
stable  full  length  decorrelator  to  exist  in  the  limit  as  the 
message  length  becomes  infinite,  and  therefore  any  set  of 
signature  sequences  not  satisfying  equation  (2.58)  should  be 
discarded  for  use  with  the  windowed  decorrelator  as  well. 

In  an  attempt  to  gain  some  insight  into  the  K user 
scenario,  a numerical  approach  was  again  taken.  A system 
with  ten  active  users  was  considered  with  signature 
sequences  derived  from  Gold  codes  of  length  31  chips.  For 
convenience,  a three  bit  wide  window  was  used.  Ten  Gold 
codes  were  selected  at  random  from  the  set  of  33  possible 
codes,  and  nine  relative  delays  were  randomly  chosen  as 
integers  between  0 and  30  chips.  This  selection  process 
yielded  one  particular  timing  configuration  for  the  ten  user 
system,  and  the  corresponding  matrix  A was  computed  along 
with  its  eigenvalues.  This  procedure  was  then  repeated 
10,000  times,  thus  yielding  10,000  different  timing 
configurations.  For  each  realization,  the  magnitude  of  the 
dominant  eigenvalue  of  the  matrix  A , , was  computed 

and  compared  with  unity.  If  |A,max|  < 1,  then  that  particular 

timing  configuration  was  labeled  as  "good";  otherwise  it 
received  the  label  "bad."  The  results  of  the  10,000  trials 
are  displayed  in  Table  2.2,  where  it  is  seen  that 
approximately  23%  of  the  tested  timing  configurations  fell 
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Table  2.2  - Numerical  Investigation  of  10,000  Timing 
Configurations  For  a Ten  User  Asynchronous  System 


No . of 
■good" 
timing 
scenarios 

No . of 
■bad" 
timing 
scenarios 

No . of 
■necessary 
not 

satisfied" 

No . of 
"sufficient 
not 

satisfied" 

1^- maxl^ 

max  |mir 

| max | aVg 

7661 

(76.61%) 

2339 

(23.39%) 

0 

2339 

0.99998 

0.27306 

0.75981 

into  the  "bad"  category,  meaning  that  under  those  timing 
configurations  the  windowed  decorrelator  would  not  converge 
to  a near- far  resistant  receiver.  Of  course,  there  may  not 
always  be  a stable  infinite  length  decorrelator  for  the 
windowed  decorrelator  to  converge  to;  therefore,  for  each 
"bad"  timing  configuration  found,  equation  (2.58)  requires 
checking.  Due  to  the  fact  that  the  variable  go  is  defined 
along  a continuum  of  values,  numerically  evaluating  equation 
(2.58)  is  not  possible  and  an  alternative  method  must  be 
used.  It  has  been  shown  by  Lupas  [61]  that  a necessary 
condition  for  equation  (2.58)  is 

det[R(0)  + R(l)  + RT(l)]*0  and  det[R(0)-R(l)-RT(l)]*0  (2.76) 


and  a sufficient  condition  for  equation  (2.58)  is 


Kin  (R(0))  > max{l!MX  (R  + ),  x2min  (R+ )}+  K*x  (R-  ) 

(2.77) 

where 

R+  =Rt(1)  + R(1) 

(2.78) 

and 

R_=j[RT(l)-R(l)] 

(2.79) 

For  each  "bad"  timing  configuration  found 

in  Table  2.2, 

equations  (2.76)  and  (2.77)  were  evaluated;  those  cases  for 
which  the  necessary  condition  of  equation  (2.76)  was  not 
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satisfied  were  labeled  "necessary  not  satisfied",  and  those 
for  which  the  sufficient  condition  of  equation  (2.77)  was 
not  satisfied  were  labeled  "sufficient  not  satisfied."  As 
shown  in  Table  2.2,  it  was  found  that  the  necessary 
condition  for  the  existence  of  the  infinite  length 
decorrelator  was  always  satisfied  for  the  "bad"  timing 
configurations,  whereas  the  sufficient  condition  was  never 
satisfied.  Thus,  the  results  of  the  numerical  investigation 
listed  in  Table  2.2  cannot  be  interpreted  as  conclusive 
proof  that  the  infinite  length  decorrelator  did  not  exist 
whenever  the  windowed  decorrelator  failed  to  converge  to  a 
near-far  resistant  receiver.  The  fact  that  the  sufficient 
condition  was  never  satisfied  when  the  windowed  decorrelator 
failed,  however,  together  with  the  proof  given  for  the  case 
of  two  users,  strongly  implies  this  relationship.  It  is 
therefore  conjectured  that  as  long  as  a stable  version  of 
the  full  length  decorrelator  exists  in  the  limit  as  the 
message  length  becomes  infinite,  then  the  windowed 
decorrelator  will  always  converge  to  a near-far  resistant 
receiver  as  the  window  index  increases.  This,  of  course, 
means  that  equation  (2.58)  represents  the  conditions  imposed 
upon  the  signature  sequences  for  the  windowed  decorrelator 
as  well  as  the  full  length  decorrelator. 

Also  listed  in  Table  2.2  are  the  largest,  smallest,  and 
average  values  found  from  the  10,000  timing 

configurations  tested.  This  information  is  of  interest 
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since  the  value  of  the  dominant  eigenvalue  will  control  the 
rate  at  which  the  asymptotic  efficiency  of  the  windowed 
decorrelator  will  converge  to  a constant  value.  The  smaller 

the  value  of  , the  faster  the  convergence  of  the 

windowed  decorrelator  will  be. 

Since  the  relative  time  delays  between  the  users' 
received  signals  are  defined  along  a continuum  of  values,  an 
infinite  number  of  possible  timing  configurations  exists; 
fortunately,  it  was  reasoned  by  Lupas  [61]  that  the 
probability  of  error  based  upon  integer  chip  delays  will  be 
worse  than  for  continuous  delays.  Indeed,  this  claim  was 
proven  for  the  two  user  case  and  numerical  results  support 
the  extension  to  more  users  [61] . Thus,  restricting  the 
delays  in  the  previous  numerical  analysis  to  be  integer 
multiples  of  a chip  not  only  reduces  the  number  of  total 
timing  configurations  possible,  but  also  represents  a worst 
case  scenario  with  respect  to  the  time  delays. 

The  results  presented  in  Table  2.2  would  suggest  that 
the  decorrelating  solution  will  not  be  applicable 
approximately  23%  of  the  time  for  a system  with  ten  active 
users  employing  Gold  codes  of  length  31  chips  as  their 
signature  sequences.  It  is  important  to  recall,  however, 
that  the  analysis  of  this  chapter  has  been  based  upon  each 
user  employing  baseband  polar  signaling,  since  it  was  shown 
in  section  2.1  that  baseband  transmission  represents  worst 
case  performance  compared  to  RF  transmission  due  to  the  fact 
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that  the  signature  waveform  cross-correlations  are  maximized 
for  baseband  transmission.  Therefore  the  results  presented 
in  Table  2.2  should  be  considered  to  be  worst  case  results 
since  they  represent  a combination  of  baseband  transmission 
and  integer  chip  delays.  To  obtain  results  more 
representative  of  a practical  communication  system,  the 
previous  numerical  analysis  was  repeated  for  a system 
employing  RF  signature  waveforms.  For  each  of  the  10,000 
timing  configurations  tested,  each  user  was  also  assigned  a 
random  received  carrier  phase  uniformly  distributed  over  the 
interval  [0, 27t)  . The  results  of  the  RF  system  investigation 
are  presented  in  Table  2.3,  where  it  is  seen  that  no  timing 
configurations  were  encountered  for  which  the  decorrelating 
strategy  would  not  work.  Furthermore,  the  values  of 

are  now  very  small,  indicating  that  the  windowed 
decorrelator  will  converge  to  a near- far  resistant  receiver 
very  rapidly  in  an  RF  system.  Also,  the  results  in  Table 
2.3  are  still  worst  case  with  respect  to  the  time  delays, 
since  integer  chip  delays  were  again  used;  thus,  the 
preceding  analysis  would  suggest  that  the  windowed 
decorrelator  would  be  very  robust  in  a near- far  RF  scenario. 


Table  2.3  - Numerical  Investigation  of  10,000  Timing 
Configurations  For  a Ten  User  Asynchronous  RF  System 


No . of 
"good" 
timing 
scenarios 

No.  of 
"bad" 
timing 
scenarios 

No.  of 
"necessary 
not 

satisfied" 

No.  of 
'sufficient 
not 

satisfied" 

1^- maxima 

l^max  |mir 

l^maxlavg 

10,000 
(100  %) 

0 

(0  %) 

0 

0 

0.1800 

0.0286 

0.0777 
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Numerical  Examples 

In  this  section  several  numerical  examples  are 
presented  to  highlight  different  features  of  the  windowed 
decorrelator  and  to  compare  its  performance  to  that  of  the 
conventional  DS-CDMA  detector  and  the  full  length 
decorrelator.  All  examples  represent  worst  case  performance 
with  respect  to  carrier  phases  and  time  delays,  as  baseband 
waveforms  and  integer  chip  delays  are  used. 

For  the  first  example,  the  dependence  of  the  asymptotic 
efficiency  of  the  windowed  decorrelator  on  bit  position 
within  the  window  was  investigated  for  a system  with  two 
active  users.  A window  of  length  five  bits  was  used,  and 
the  cross-correlations  were  arbitrarily  chosen  to  be 
Ri2(0)=  0-3  and  R12(l)  = 0.35 . Equation  (2.55)  was  used  to  plot 
the  asymptotic  efficiency  of  user  2 for  the  bits  contained 
within  the  second  window  processed  by  the  receiver.  The 
results  for  all  five  bit  positions  v=2,4,6,8,  and  10  are 
shown  in  Figure  2.8,  from  which  several  interesting  features 
of  the  windowed  decorrelator  are  apparent.  First, 
regardless  of  the  bit  position,  all  the  asymptotic 
efficiency  curves  decay  to  a value  of  zero;  consequently, 
the  near-far  resistance  for  all  the  bits  processed  in  the 
second  window  is  zero.  Define  the  near-far  spread,  Wf  , to 

be  the  largest  difference  in  received  energies  for  which  the 
windowed  decorrelator  has  a nonzero  asymptotic  efficiency. 

As  long  as  the  near-far  spread  is  larger  than  the  worst  case 
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Figure  2.8  - Asymptotic  Efficiency  as  a Function  of  Bit 
Position  for  Two  Users  with  a Five  Bit  Wide  Window 

difference  in  received  energies  encountered  in  practice, 
then  for  all  practical  purposes  the  windowed  decorrelator 
will  be  near-far  resistant.  Secondly,  it  is  seen  in  Figure 
2.8  that  there  can  be  a large  difference  in  the  value  of 
near-far  spread  obtained  for  different  bit  positions  within 
the  same  window.  In  this  example,  the  value  of  near- far 
spread  ranges  from  approximately  24  dB  for  positions  v=2  and 
v=10  to  approximately  53  dB  for  v=6.  Note  that  the  largest 
near-far  spread  occurs  for  the  bit  in  the  middle  of  the 
window,  whereas  the  smallest  near- far  spread  occurs  for  the 
first  and  last  bits  within  the  window  (the  edge  bits) . This 
phenomenon  makes  sense  since  the  bits  at  either  edge  of  the 
window  are  directly  interfered  with  by  the  out-of -window  MAI 
and  therefore  should  feel  the  largest  amount  of  out-of- 
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window  MAI  influence.  Multiplying  the  window  vector  by  the 
inverse  of  the  cross-correlation  matrix  has  the  effect  of 
distributing  the  out-of -window  MAI  across  all  the  bits 
within  the  window;  however,  bits  further  away  from  the 
source  of  the  out-of-window  MAI  will  be  affected  less.  The 
fact  that  there  is  out-of-window  MAI  both  preceding  and 
following  each  window,  with  the  exceptions  of  the  very  first 
and  last  windows,  accounts  for  why  bits  near  the  middle  of 
the  window  will  have  the  largest  near- far  spread  and  why 
there  is  a high  degree  of  symmetry  about  the  center  of  the 
window  (positions  v=2  and  v=10  have  the  same  near-far 
spread,  as  do  positions  v=4  and  v=8) . Thirdly,  another 
interesting  feature  of  the  windowed  decorrelator  is  that  the 
maximum  value  of  asymptotic  efficiency  varies  with  bit 
position  as  well.  From  equation  (2.55),  the  maximum  value 
of  asymptotic  efficiency  will  be  l/cvv  , and  for  this  example 
l/c22  =0.7732  , l/c4i4  =0.7592  , l/c66  = 0.7590  , l/c8  8 = 0.7611  , and 
l/cio.io  =0-8957.  In  summary,  the  performance  of  the  windowed 
decorrelator  will  be  limited  by  its  performance  for  its 
worst  case  bit  position,  which  will  occur  for  either  the 
first  or  the  last  bit  in  the  window,  with  approximately  even 
symmetry  about  the  center  of  the  window. 

The  second  example  illustrates  the  convergence  of  the 
windowed  decorrelator  to  a full  length  decorrelator  as  the 
number  of  windows  processed  increases.  A two  user  system  is 
again  investigated,  and  a window  length  of  three  bits  is 
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used.  The  cross-correlations  are  R12(0)=2/7  and  R12(l)  = l/7, 
corresponding  to  signature  sequences  derived  from  Gold  codes 
of  length  seven  chips  with  a fixed  relative  delay  of  one 
chip.  For  this  example,  the  initial  amplitude  estimation 
error  vector  is 


’|AW!(3Lx" 

|C5,6R2l(-1)|VWl(4) 

0.0455^,(4) 

JAW2(3Lx_ 

0.1559^,(4) 

and  the  matrix  A is 


C5,6R2l(  -0 

C5,lR12(-0 

'0.0455 

0.0003' 

C6,6R2l(-1) 

C6,1R12  C1) 

0.1559 

0.0001_ 

The  eigenvalues  of  A are  easily  found  to  be  A.,  =0.0466  and 
X,2  =-0.0010;  since  equation  (2.70)  is  satisfied,  the 

convergence  of  the  windowed  decorrelator  to  a full  length 
decorrelator  is  guaranteed.  Figure  2.9  displays  the 
asymptotic  efficiency  of  the  windowed  decorrelator  for  bits 
within  the  second,  third,  and  fourth  windows  processed. 

These  curves  represent  the  worst  case  bit  position,  which 
was  found  to  be  v=2  for  this  example.  The  asymptotic 
efficiency  of  the  conventional  detector  is  also  shown,  as  is 
the  asymptotic  efficiency  for  a window  index  of  infinity. 

For  this  example  the  conventional  detector  has  a near- far 
spread  of  only  7.3  dB,  whereas  the  windowed  decorrelator  has 
near-far  spreads  of  43  dB,  70  dB,  and  96  dB  for  the  second, 
third,  and  fourth  windows,  respectively.  Furthermore,  note 
that  even  if  perfect  power  control  was  used  in  conjunction 
with  the  conventional  detector,  the  best  asymptotic 
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Figure  2.9  - Asymptotic  Efficiency  as  a Function  of  Window 
Index  for  Two  Users  with  a Three  Bit  Wide  Window 


efficiency  it  could  achieve  would  be  approximately  0.3.  On 
the  other  hand,  the  convergence  of  the  windowed 
decorrelator's  asymptotic  efficiency  to  a constant  value  of 
0.91  is  readily  apparent,  and  thus  the  windowed  decorrelator 
is  capable  of  obtaining  an  asymptotic  efficiency 
approximately  three  times  larger  than  the  conventional 
detector  with  perfect  power  control.  This  convergence  is 
even  more  readily  apparent  from  Figure  2.10,  which  displays 
the  near-far  spread  versus  window  index.  Indeed,  after  only 
the  fourth  window,  corresponding  to  processing  only  12 
received  data  bits  in  this  example,  the  windowed 
decorrelator  is  able  to  tolerate  differences  in  received 
energy  of  up  to  96  dB. 


60 


In  Figure  2.11,  the  results  of  the  second  example  are 
extended  to  include  the  effects  of  additional  users.  Here 
the  near-far  spread  of  user  2 is  plotted  as  a function  of 
the  number  of  active  users.  The  cross-correlations  between 
all  users  were  kept  equal;  that  is,  Rjlc (0)=  2 / 7,  j ^ k and 
Rjk(l)  = l/7,k> j . Also,  all  users  except  for  the  desired  user 

(user  2)  were  assumed  to  be  at  the  same  received  power 
level,  thus  representing  a worst  case  scenario.  For 
instance,  for  the  case  of  five  active  users,  a near-far 
spread  of  30  dB  means  that  the  four  other  users  each  had  a 
received  power  level  30  dB  greater  than  that  of  the  desired 
user.  As  can  be  seen  from  Figure  2.11,  the  windowed 
decorrelator  maintains  its  ability  to  achieve  any  desired 
level  of  near- far  spread  as  the  number  of  active  users 


Figure  2.10  - Near-far  Spread  versus  Window  Index  for  the 
System  of  Example  Two 


180  T 
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Figure  2.11  - Near-far  Spread  as  a Function  of  the  Number  of 
Active  Users  for  the  System  of  Example  Two 

increases;  however,  the  rate  of  obtaining  this  near- far 
spread  decreases  with  each  additional  user.  For  example,  a 
three  user  system  can  achieve  a near- far  spread  of  90  dB 
after  processing  six  windows  (18  bits),  whereas  a four  user 
system  will  require  nine  windows  (27  bits) . It  should  be 
noted  that  for  this  particular  example  in  which  the  cross- 
correlations were  derived  from  Gold  codes  of  length  seven 
chips,  the  case  of  five  active  users  represents  a very 
heavily  congested  system. 

The  third  example  is  identical  to  the  second  example, 
but  now  the  cross-correlations  are  given  larger  magnitudes. 
Now  R12(0)=0.3  and  R12(l)  = 0.5  are  arbitrarily  chosen  and  the 
matrix  A becomes 

[0.2402  0.0312' 

~ |_0.5  720  0.00  9 3_ 

The  eigenvalues  of  A are  now  found  to  be  =0.3012  and 
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l2  =-0.0517  , so  that  again  equation  (2.70)  is  satisfied  and 

convergence  is  guaranteed.  The  asymptotic  efficiencies  for 
the  second,  third,  and  fourth  windows  are  displayed  in 
Figure  2.12.  Note,  however,  that  the  rate  of  convergence  is 
not  as  rapid  as  in  the  previous  example,  due  to  the  higher 
cross-correlations.  Also,  the  maximum  value  of  asymptotic 
efficiency,  0.6246,  is  significantly  lower  than  in  the 
previous  example.  The  conclusion  is  that  even  when  a 
decorrelating  strategy  is  used,  it  is  still  important  to 
choose  signature  waveforms  with  low  cross-correlations. 


Figure  2.12  - Asymptotic  Efficiency  of  System  in  Figure  2.9 
with  Larger  Cross-correlations 


Since  all  of  the  previous  examples  have  been  based  upon 
the  results  for  asymptotic  efficiency  derived  in  this  work, 
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the  fourth  example  is  a computer  simulation  which  serves  to 
verify  the  analytical  results.  A two  user  asynchronous 
system  was  simulated  with  the  same  cross-correlations  used 
in  the  second  example.  Constant  cross-correlations  and 
received  energies  were  used.  The  simulation  sent  100  bits 
through  a windowed  decorrelator  and  found  the  magnitude  of 
the  final  decision  statistic  zL  given  by  equation  (2.11). 
This  decision  statistic  is  the  sum  of  three  terms:  the 
desired  signal  component,  the  remaining  MAI  component,  and 
the  noise  component.  If  the  receiver  is  operating  as 
expected,  the  MAI  component  in  the  decision  statistic  will 
become  smaller  as  the  window  index  increases,  until  it 
eventually  disappears.  Thus,  |zK(i)|  should  converge  to 

Vwk(i)  plus  zero-mean  noise.  Figure  2.13  shows  the  results 

of  the  simulation  when  no  background  noise  was  present  and 
201og(wl/w2)  = 96  dB.  The  solid  line  indicates  the  value  of 
201og-y/w2  / chosen  to  be  8 here.  Indeed,  after  approximately 

13  received  data  bits  the  receiver  has  converged  for  user  2, 
agreeing  very  strongly  with  the  analytical  results  presented 
in  the  second  example. 

Figure  2.14  shows  the  results  of  the  simulation  when 
user  2 is  operating  at  a signal- to-noise  ratio  of  8 dB  and 
the  near-far  effect  is  again  96  dB.  Here  the  magnitude  of 
the  decision  statistic  converges  to  a noisy  value  of  ^w^  . 

Note  that  Figure  2 . 14  was  plotted  on  a log  scale  due  to  the 
enormous  near-far  effect  of  96  dB;  for  this  reason,  Figure 
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Figure  2.13  - Simulation  Results  for  the  Magnitude  of  the 
Decision  Statistic  vs.  Number  of  Received  Bits  for  Two  Users 
with  a Near-far  Effect  of  96  dB,  and  no  Background  Noise 


Received  Bit  Index  i 

Figure  2.14  - Results  for  the  Simulation  of  the  Decision 
Statistic  Magnitude  vs.  Number  of  Received  Bits  with  User  2 
Having  a Signal-to-noise-ratio  of  8 dB 
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Figure  2.15  - Steady-State  Portion  of  Figure  2.14  on  a 
linear  scale 


2.15  displays  the  steady-state  portion  of  Figure  2.14  on  a 
linear  scale  so  that  it  is  more  readily  apparent  that  |z2(i)| 
has  converged  to  plus  zero-mean  noise. 


CHAPTER  3 

WINDOWED  DECORRELATING  MULTIUSER  RECEIVER  FOR  FREQUENCY- 
NONSELECTIVE  RAYLEIGH  FADING  CHANNELS 

System  Description 

In  the  presentation  of  the  windowed  decorrelator  for 
AWGN  channels  given  in  chapter  2,  it  was  assumed  that  the 
multiuser  receiver  had  knowledge  of  the  received  carrier 
phase  of  each  user,  thus  allowing  for  coherent  demodulation 
of  the  BPSK  modulated  signals.  In  an  AWGN  channel,  the 
received  phases  will  be  changing  slowly  relative  to  a symbol 
interval,  and  tracking  of  the  phases  may  therefore  be 
possible.  It  should  be  noted,  however,  that  the  standard 
single-user  phase-lock-loop  structure  will  not  be  adequate 
to  accomplish  this  task  due  to  the  presence  of  the  MAI; 
thus,  some  other  method  of  obtaining  knowledge  of  the  phases 
will  be  required.  In  fading  channels  the  received  phases 
will  be  changing  much  more  rapidly,  in  many  cases  precluding 
the  possibility  of  successfully  tracking  these  quantities. 
Therefore,  DPSK  will  be  the  assumed  modulation  format  of  all 
users  in  the  fading  channel.  As  mentioned  in  chapter  2,  the 
use  of  DPSK  in  conjunction  with  the  windowed  decorrelator  is 
somewhat  self-defeating,  since  the  received  phases  of  all 
users  will  need  to  be  estimated  between  consecutive  windows 
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in  order  to  correct  for  the  out-of -window  MAI.  DPSK  does, 
however,  allow  for  the  data  bits  contained  within  each 
window  to  be  demodulated  without  knowledge  of  their 
corresponding  phases,  and  a phase  estimation  scheme  is 
proposed  in  this  chapter  which  allows  the  received  phases  of 
bits  between  consecutive  windows  to  be  estimated  while 
requiring  no  additional  computations  beyond  those  required 
for  decorrelating  the  current  window. 

The  system  model  for  the  fading  channel  is  illustrated 
in  Figure  3.1.  Each  user  transmits  a sequence  of 
differentially  encoded  data  bits  through  the  channel.  The 
ifck  differentially  encoded  bit  of  the  kt*1  user,  dk(i)  , is 
related  to  the  ith  data  bit  of  the  kth  user,  bk(i)  , through 
the  relationship 

dk(i)  = dk(i-l)bk(i)  (3.1) 

The  signature  waveforms  are  now  baseband  polar  sequences 
normalized  to  unit  energy  and  zero  outside  the  time  interval 
[ 0 , T ] , where  T is  again  the  common  symbol  duration.  Each 
user  experiences  statistically  independent  flat  Rayleigh 
fading,  represented  as  a time-varying  multiplicative 
distortion  of  each  user's  transmitted  signal.  The 
multiplicative  coefficients  ak(i)  are  modeled  as  independent 
zero-mean,  complex-valued  Gaussian  random  variables.  Since 
the  fading  processes  of  different  users  are  statistically 
independent,  the  following  equality  holds 

E{a*(i)ak(i)}=0,  Vj*k 


(3.2) 
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Figure  3.1  - System  Model  for  Flat  Rayleigh  Fading  Channel 


The  complex  envelope  of  the  received  signal  at  the  base 
station  is  written  as 

R(t)  = s(t,d)+N(t)  (3.3) 

where  N(t)  is  complex  zero-mean  Gaussian  noise  with  power 
spectral  density  a2  and 

S(t,  b)  = 2 Idk(i)ak  (i)Vwk(i)Sk (t  - iT  - Tk  )ejMi)  (3.4) 

i=l  k=l 

where  ^/wk(i)  is  the  received  power  of  the  user  in  the 

it*1  symbol  interval,  sk(t)  is  the  signature  sequence  of  the 
kt^L  user,  and  (J)k(i)  is  the  received  phase  of  the  k^h  user  in 
the  ith  symbol  interval.  The  users  are  again  ordered 
according  to  their  relative  time  delays  0 < x1  < x2  <•••<  tK  < T . 
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The  cross-correlation  matrices  are  defined  by  their 
(m,n)t^1  element  as 

oo 

R ».n  (1)  = j sm  (t  - Tm  )ln  (t  + IT  — xn  )lt  (3.5) 

— oo 

where  it  is  again  clear  from  the  time-limited  nature  of  the 
signature  sequences  that  R(1)  = 0 for  |1|  > 1 and  R(-1)  = RT(1). 

The  front  end  of  the  multiuser  receiver  will  again 
consist  of  a matched  filter  bank,  one  filter  being  matched 
to  each  allowed  signature  waveform.  The  output  of  the  kfck 
user's  matched  filter  sampled  at  the  end  of  the  ith  bit 
interval  is 

(i+l^+i* 

yk(i)=  Jr(t)Ik(t-iT-Tk)dt  (3.6) 

iT+1k 

The  vector  of  matched  filter  outputs  sampled  at  the  end  of 

the  i^h  bit  interval  can  be  written  as 

y(i)  = R(— l)w(i  + l)<D(i  + l)A(i  + l)d(i  + 1)  + R(0)w(i)O(i)A(i)d(i) 

+ R(l)w(i  - l)0(i  - l)A(i  - l)d(i  - 1)  + n(i)  ( 3 ' 7 } 

where  d(i)  = [d , (i),  d2  (i),  • • • , dK  (i)]T  , A(i)  = diag[ai  (i),  a 2 (i),  • • • , a K (i)]  , 

w(i)  = diag^w1(i),  Vw2(i),-- • , ^wK(i)] , O(i)  = diagfe^1  (i) , ej*2(i) , ••  • , ej^(i) ] , 
and  n(i)  = [n1(i),n2(i),”-/nK(i)]T . The  covariance  matrix  of  the 
matched  filter  output  noise  vector  is 

|E{n*(i)nT(j)}  = a2R(i-j)  (3.8) 

For  an  L bit  wide  window  spanning  bits  i through  i+L-1, 
the  entire  sequence  of  matched  filter  outputs  can  be  written 
as 

Yl  = RLKWL0LALdL  + nL  + PWeOeAede 


(3.9) 
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where  the  quantities  WL  =diag[w(i),W(i+l),---,W(i+L-l)] , 

Ol  =diag[0(i),0(i  + l),*“,0(i  + L-l)] , AL  = diag[A(i),A(i+l),---,A(i+L-l)] , 

dL  =[d(i)T<d(i+l)T,---,d(i+L-l)T]T  , nL  =[n(i)T,H(i  + l)T,---,H(i+L-l)T]T  , 

W e =diag[w(i-l),w(i+L)]  f Oe  = diag[0(i-l),"-0(i+L-l)] , 

A.  =diag[A(i-l),A(i  + L-l)],  de  =[d(i-l)T,d(i  + L)T]T  , and  the  LK  by 
LK  and  LK  by  2K  matrices  RLKand  P are  given  by 


rR(0)  R(— 1)  0 •••  0 N 

'R(l) 

0 ' 

R(l)  R(0)  R(— 1)  : 

0 

0 

0 \ *•.  0 

and  P= 

• 

• 

: R(l)  R(0)R(— 1) 

0 

0 

, 0 - 0 R(1)R(0), 

, o 

R(-l), 

Proposed  Multiuser  Receiver  Structure 

The  structure  of  the  windowed  decorrelator  for  a flat 
Rayleigh  fading  channel  is  essentially  the  same  as  for  an 
AWGN  channel,  with  two  exceptions.  First,  the  estimate  of 
the  out-of -window  MAI  affecting  each  window  must  now  include 
an  estimate  of  each  user's  phase  immediately  preceding  and 
following  the  window  (recall  that  the  phases  were  assumed  to 
be  known  in  the  AWGN  channel) . Secondly,  each  output  vector 
of  the  windowed  decorrelator  will  contain  differentially 
encoded  data  bits  of  each  user;  therefore  the  windowed 
decorrelating  algorithm  must  be  followed  by  a bank  of 
standard  DPSK  detectors,  one  for  each  user  as  shown  in 
Figure  3.2.  The  same  technique  can  be  used  as  in  the  AWGN 
channel  to  form  the  required  MAI  estimates,  only  now  the 
intermediate  decision  statistic  vector  resulting  from 
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decorrelating  each  window  will  contain  complex  numbers.  The 
last  K elements  in  the  current  window's  intermediate 
decision  statistic  vector  will  be  used  as  the  estimates  of 
the  out-of-window  MAI  affecting  the  next  window:  the 

magnitude  of  each  complex  number  will  be  the  estimated 
magnitude  of  the  corresponding  user's  MAI  and  the  phase  of 
each  complex  number  will  be  the  estimated  phase  of  the 
corresponding  user's  MAI,  as  illustrated  in  Figure  3.3. 

This  estimation  technique  will,  however,  result  in  a 180° 
phase  ambiguity  due  to  the  fact  that  the  polarity  of  the 
differentially  encoded  data  bits  may  be  either  +1  or  -1  and 
are  unknown  to  the  receiver.  One  way  to  remove  this  phase 
ambiguity  is  to  have  the  last  bit  of  every  user  within  each 
window  be  a pilot  bit.  The  polarities  of  the  pilot  bits  are 
known  so  there  will  be  no  phase  ambiguity  in  the  MAI 
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ML  = MAI  magnitude  estimate  for  user  1 
<t>1  = MAI  phase  estimate  for  user  1 

M^ = MAI  magnitude  estimate  for  user  2 
<J>2  = MAI  phase  estimate  for  user  2 


M^  = MAI  magnitude  estimate  for  user  K 
(j)  = MAI  phase  estimate  for  user  K 

K 


Figure  3.3  - MAI  Estimation  Technique  of  Fading  Channel 
Windowed  Decorrelator 


estimates.  Each  user  must  still  transmit  a pilot  bit 
between  consecutive  data  windows  so  that  the  out-of-window 
MAI  may  be  subtracted  off  properly,  and  thus  each  user  in 
the  fading  channel  must  transmit  two  consecutive  pilot  bits 
after  every  block  of  L-l  data  bits,  as  shown  in  Figure  3.4. 
The  first  pilot  bit  is  grouped  with  the  L-l  differentially 
encoded  data  bits  to  form  a window  L bits  wide;  this  window 
is  then  decorrelated.  The  second  pilot  bit  serves  to 
separate  the  windows  and  will  be  discarded  by  the  receiver, 
just  as  the  lone  pilot  bit  was  in  the  AWGN  analysis. 


L-l  Data 

Pilot 

Pilot 

L-l  Data 

Pilot 

Pilot 

Bits 

Bit 

Bit 

Bits 

Bit 

Bit 

Decorrelate 


Discard 


Decorrelate 


Discard | 


Figure  3.4  - Data  Block  Structure  of  the  Mobile  Units 


The  estimation  technique  described  above  will  be 
referred  to  as  a zeroth  order  scheme  since  it  assumes  each 
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user's  fading  process  will  be  approximately  constant  over  a 
span  of  L+2  consecutive  symbol  intervals.  A smarter 
estimation  technique  might  try  to  predict  future  fading 
coefficients  for  each  user  based  upon  past  values.  In  order 
to  choose  an  appropriate  prediction  algorithm,  some 
characteristics  must  be  known  about  the  fading  processes.  A 
computer  simulation  of  the  magnitude  and  phase  of  a typical 
flat  Rayleigh  fading  process  is  shown  in  Figure  3.5  as  a 
function  of  the  transmitted  bit  index.  To  represent  numbers 
which  might  be  typical  of  a realistic  wireless  communication 
system,  a transmitted  bit  rate  of  20  kbps  and  a carrier 
frequency  of  900  MHz  were  chosen.  In  addition,  the  user  is 
assumed  to  be  traveling  at  a speed  of  65  mph,  thereby 
representing  the  worst  case  fading  which  might  be 
encountered  by  a user  in  a cellular  phone  system.  As  can  be 
seen,  in  between  the  maxima  and  minima  the  curves  are 
approximately  linear;  for  instance,  between  transmitted  bits 
460  and  540  of  Figure  3.5.  Linear  extrapolation  would 
therefore  seem  a logical  choice  for  a prediction  technique. 
Thus,  as  an  alternative  to  using  the  zeroth  order  estimation 
strategy,  linear  extrapolation  could  be  employed.  The 
multiuser  receiver  would  need  to  store  the  last  K elements 
in  the  previous  window's  decorrelated  output  vector.  After 
decorrelating  the  current  window  the  receiver  will  have  two 
points,  each  a complex  number,  for  each  user  from  which  it 
can  then  compute  a slope  and  an  intercept.  The  estimate  of 
each  user's  MAI  contribution  to  the  next  window  can  then  be 
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Figure  3.5  - Magnitude  and  Phase  of  a Flat  Rayleigh  Fading 
Process 

determined  by  plugging  the  appropriate  future  transmitted 
bit  indices  into  the  equation  for  that  user's  line.  Of 
course,  this  estimation  scheme  will  produce  poor  estimates 
near  the  maxima  and  minima  shown  in  Figure  3.5,  but  these 
points  occur  relatively  infrequently,  particularly  for 
smaller  window  lengths.  To  facilitate  the  following 
performance  analysis,  the  zeroth  order  estimation  strategy 
will  be  assumed.  The  linear  extrapolation  method  will  be 
referred  to  as  the  first  order  scheme  later  in  this  chapter. 
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The  overall  block  diagram  of  the  windowed  decorrelating 
receiver  for  flat  Rayleigh  fading  channels  is  shown  in 
Figure  3.6.  The  vector  of  decorrelated  outputs  is  given  by 


Figure  3.6  - Block  Diagram  of  Fading  Channel  Windowed 
Decorrelator 
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(3.10) 


Defining 

yL=R^HL  (3.11) 

the  k^h  user's  decorrelated  intermediate  decision  statistic 
is  given  by 

dk(i)=A/wk(i)ak(i)e^(i)dk(i)+Yk(i)+[R^p(weOeAede-WeOeAed)]v 


(3.12) 
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where  as  before  the  notation  [x]“  represents  the  element 

in  the  u^  window  of  the  vector  x . As  in  the  AWGN 
analysis,  the  window  index  u and  the  window  element  index  v 
are  given  by  equations  (2.13)  and  (2.14),  respectively.  The 
final  bit  decisions  are  determined  as 

(i) = s ign^Re  jdk  (i  - l)d*  (i)||  (3.13) 

Performance  Analysis 


The  decision  variable  for  the  ifch  bit  of  the  user 

given  by  equation  (3.13)  can  be  written  in  the  form 

D = Re{xY*}  (3.14) 

and  therefore  the  probability  of  error  is  given  by  [55, 
Appendix  B] 


Pk(i)=Kl-H) 


where 


M + M 

xy  yx 


I 


M -M  ) +4M  M 

xy  yx ) xx  yy 


and  the  covariances  M„,MW,MW, 

yy  *y 

Mxx=|e{xx*} 
Myy=|E{YY*} 
Mxy  =iE{xY*} 
Myx=|E{YX*} 

As  a starting  point,  it  is 


(3.15) 

(3.16) 

and  are  given  by 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

of  interest  to  determine  the 


best  performance  possible  from  the  proposed  windowed 
decorrelator,  corresponding  to  the  ideal  scenario  in  which 
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the  multiuser  receiver  has  perfect  knowledge  (or  perfect 
estimates)  of  the  out-of -window  MAI.  In  this  scenario 


equation  (3.12)  reduces  to 

dk(i)  = ylwy(i)ak(i)e3*Al)ay(i)  + Yk(i)  (3.21) 

From  equations  (3.13)  and  (3.14),  the  variables  X and  Y are 
then  given  by 

X = Vwk(i-l)a  k (i  - l)e  (i-1}d  k (i  - 1)  + y k (i  - 1)  (3.22) 

Y = Jwjjija  k(i)ej*k(i)dk  (i)+Yk(i)  (3.23) 

Applying  equation  (3.17)  to  find  yields 

Mxx=iwk(i-l)E{ak(i-l)a;(i-l)}  + iE{Yk(i-l)Y;(i-l)}  (3.24) 

where  the  fact  that 

E{ak(i)Yk(3)}=E{a;(i)Yk(j)}  = 0 V i,j  (3.25) 


has  been  used  since  the  fading  processes  are  statistically 
independent  from  the  background  noise. 

The  time  autocorrelation  function  for  the  k^h  user's 
fading  process  is  defined  as 

R«k  (T) = 2 E{ak  (t)ak  (t  - t)}  (3.26) 

and  from  equations  (3.8)  and  (3.11) 

1 E{yk  (i  - l)Yk  (i  ~ 1)}  = cr2c  V.K,  V_K  (3.27) 

where  as  in  the  AWGN  analysis,  the  (ir^n)^  element  of  the 
matrix  R“k  is  denoted  by  cmn.  Using  equations  (3.26)  and 
(3.27),  the  expression  for  reduces  to 

= wk (i  - l)Rak  (0)  + o2cv_K  V_K  (3.28) 

Similarly,  Myy  is  found  to  be 

Myy  =wk(i)Rak(0)+a2v,v 

Applying  equation  (3.19)  to  find  results  in 

xy 


(3.29) 
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(3.30) 


M*v  = !i/w«(i_1)Vwk(i)  e’1*''1  “ *''u|E{ae(i  - l)a;(i)}dK(i  - l)dj(i) 

+ iE{yt(i-l)Y;(i)} 

For  convenience  it  will  be  assumed  that  the  received 
energies  and  phases  are  constant  over  two  consecutive 

signaling  intervals,  so  that  equation  (3.30)  simplifies  to 

=iwkblc(i)Rak(T)  + a2cv.K<v  (3.31) 

Similarly,  is  found  to  be 

= 1 wkbk  (i)  RSk  (T)  + a2c  v_Kv  (3.32) 

which  is  identical  to  the  expression  found  for  in 

equation  (3.31)  and  thus  Mxy=Myx.  Plugging  equations 

(3.28),  (3.29),  (3.31),  and  (3.32)  into  the  expression  for 

|i  given  by  equation  (3.16)  yields 


M, 


P = 


/m  m" 

\ XAXX 


(3.33) 


By  symmetry,  Pr(error)  = Pr(error  | bk(i)  = l)  and  from  equation 

(3.15)  the  probability  of  error  will  be  given  by 

f \ 


Rk  (-*-)min  „ 


1- 


wk(i)Rak  (T)+o2cv_Ki, 


^k(i)Rak(0)+d2cv_K,v_K)(wk(i)Rat  (0)  + a2cv  v) 


(3.34) 


Equation  (3.34)  represents  the  probability  of  error  obtained 
by  an  ideal  windowed  decorrelator  with  perfect  estimates  of 
the  out-of-window  MAI  affecting  each  window.  Define  the  k^h 
user's  average  received  signal-to-noise  ratio  as 

WkRa„  (0) 


SNRk  =• 


(3.35) 


and  the  k^h  user's  fading  correlation  coefficient  as 

Ra.(T) 


Pat(T)  = 


Ra,(0) 


(3.36) 
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The  fading  correlation  coefficient  is  a measure  of  how 
correlated  the  fading  is  between  consecutive  signaling 
intervals.  If  the  fading  was  constant  over  two  consecutive 
intervals,  then  pa  (T)  = l.  On  the  other  hand,  if  the  fading 
was  completely  uncorrelated  from  one  interval  to  the  next, 
then  pa(T)  = 0.  Using  the  definitions  given  in  equations 


(3.35)  and  (3.36),  equation  (3.34)  may  be  re-written  as 

( ._i  'N 


^(■Omin  — „ 


1-- 


Pak  (T)  + (SNR  k )_1  c V_K  v 


(3.37) 


(l  + (SNR  k )_1  c V_K.V_K  )(l  + (SNR  k )_1  cv>v)^ 

In  the  limit  as  the  k^h  user's  average  signal-to-noise 
ratio  becomes  infinite,  the  windowed  decorrelator's  best 
possible  probability  of  error  will  become 


lif?.P«(i).i.=i(1-P..(T)) 


(3.38) 


The  effect  of  the  fading  channel  on  the  k^h  user's 
performance  is  readily  apparent  from  equation  (3.38).  For 
pak(T)*l,  there  will  be  a fading- induced  irreducible  error 


floor,  the  level  of  which  is  dependent  upon  how  highly 
correlated  the  fading  is  between  consecutive  signaling 
intervals.  For  instance,  if  pa(T)  = 0.99,  corresponding  to 

fairly  highly  correlated  fading,  then  the  k^  user  will 
encounter  an  irreducible  error  floor  of  5x1 0-3  . Of  course, 
the  windowed  decorrelator  is  not  attempting  to  combat  the 
effects  of  the  fading;  its  goal  is  simply  to  eliminate  the 
MAI  and  thereby  achieve  the  same  probability  of  error  that  a 
single  user  would  have.  To  reduce  the  level  of  the  error 
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floor  imposed  by  the  fading,  diversity  techniques  may  be 
employed  [28]  . For  the  multiuser  scenario  considered  here, 
each  user  could  repeat  his  or  her  message  in  frequency, 
time,  or  space.  Referring  to  Figure  3.2,  the  appropriate 
re-combining  of  each  user's  replicas  would  then  occur  at  the 
output  of  the  decorrelating  algorithm  prior  to  the  final 
data  demodulation. 

The  probability  of  error  obtained  by  a single  user 
employing  DPSK  modulation  in  a flat  Rayleigh  fading  channel 
is  given  by  [28] 


' e,single  user 


1 — 


Pa(T) 


v-1 


l + (SNReff) 

In  the  limit  as  the  single  user's  signal-to-noise  ratio 
becomes  infinite,  the  probability  of  error  becomes 


(3.39) 


lim 

SNRef  f “ *°° 


‘ e,  single  user 


= f(l-P.(T)) 


(3.40) 


which  is  identical  to  the  probability  of  error  for  the  ideal 
windowed  decorrelating  given  in  equation  (3.38),  and  thus 
the  ideal  windowed  decorrelator  is  able  to  asymptotically 
achieve  the  same  probability  of  error  as  a single  user. 

The  performance  of  the  ideal  windowed  decorrelator  is 
compared  to  that  of  the  single  user  in  Figure  3.7.  A five 
user  system  with  a five  bit  wide  window  was  chosen  for  the 
windowed  decorrelator.  User  two  was  chosen  as  the  desired 
user  and  equal  cross-correlations  of  Rj)c(0)=  2 / 7,  j * k and 
Rjk(l)  = l/7,  j<k  were  assumed.  The  results  displayed  in 


Figure  3.7  are  for  window  bit  position  v=7,  corresponding  to 
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Figure  3.7  - Comparison  of  Ideal  Windowed  Decorrelator  and 
Single  User 

the  first  data  bit  of  the  second  user  within  the  window  (the 
second  differentially  encoded  bit  within  the  window) . As 
can  be  seen,  the  windowed  decorrelator  with  perfect 
knowledge  of  the  out-of-window  MAI  is  able  to  track  the 
probability  of  error  of  a single  user  very  closely. 
Interestingly,  the  results  shown  in  Figure  3.7  were  found  to 
be  almost  identical  for  different  values  of  window  bit 
position  v,  indicating  that  the  ideal  windowed 
decorrelator's  performance  is  for  all  practical  purposes 
independent  of  bit  position.  This  behavior  is  expected, 
since  the  ideal  windowed  decorrelator  is  able  to  completely 
subtract  off  the  out-of-window  MAI.  Figure  3.7  also  displays 
the  error  floor  associated  with  the  fading  for  pa(T)  = 0.99  . 
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To  analyze  the  performance  of  the  windowed  decorrelator 
which  forms  its  MAI  estimates  according  to  the  "dumb" 
estimation  strategy  described  in  section  3.2,  the  last  term 
on  the  right  hand  side  of  equation  (3.12)  must  be  evaluated. 
To  begin,  consider  the  very  first  window  (u=l) . Since  there 
have  been  no  previous  windows,  the  receiver's  estimates  are 
zero  and  the  vfch  element  of  the  MAI  term  becomes 


[R^PWe<DeAede]^  = 


K— 1 LK 

X Xcv.jRj-(L-DK,i(-1)Vwi(L  + l)ejML+1)a1(L  + 1) 


1=1  j=(L— UK+l+l 


(3.41) 

where  it  has  again  been  assumed  for  convenience  that  all 
pilot  bits  have  a polarity  of  +1,  thus  allowing  the  de 
terms  to  be  dropped.  For  compactness,  an  interference 
vector  is  defined  for  the  first  window  as 

i^T)  = [R^PWe<DeAe]u=1  (3.42) 

The  estimate  of  the  user's  contribution  to  the  MAI 
affecting  the  second  window  is  formed  as 

wk  (L  + l)ej$k(L+1)ak  (L  + 1)  = wk  (2L  + 2)ed"4’k(2I,+2)ak  (2L  + 2) 

= w,(L)e»-at(L)+,>(L)+^](wk  <3-«> 

Now  consider  the  second  window  (u=2) . Recall  that  the 
decorrelated  output  vector  is  given  by  equation  (3.12), 
which  is  repeated  below  for  convenience 

(i)  = VMiK (i)ej*k{i)dk (i)  + yk (i)  + [r^ p(weOe Aede  - WeOe Aed)]“  (3.12) 

To  proceed  further,  the  last  term  on  the  right-hand  side  of 
equation  (3.12)  must  be  evaluated.  This  term  represents  the 
MAI  corruption  in  the  intermediate  decision  statistic 
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resulting  from  the  fact  that  the  channel  characteristics 
immediately  before  and  after  the  second  window  are  not  the 
same.  By  multiplying  the  matrices  together,  this  MAI  term 
may  be  expanded  out  as 


(we<PeAe-WeO>eAe)= 


>/w1(L+l)e:,*l(I,+1)a1(L+l)-A/w1(L)e;Wl(L)a1(L)-Y1(L)-  Ix(l) 
7^L^ej^(L+1)a2(  L+l)-V^e^(%(L)^2(L)-k(r) 


. (L-lJK+l 

| 

. (Ii-l)K+2 


^^IjI^e:W'lt(L+1)aK  (LHd)-.^w^(L)ej*K(L)aK  (L)-yK  (L)-^  (l)j 
yjv!1(2h+2)e^'  (2L+2)a1  (2L+2)-^/w1(L)e:<l>1  (l) 

7w2(2L+2)ej02(2L+2,a2  (2L+2)-^(L)eA<L,a2  (L)-y2  (L)- 


Ii(l) 


Ii(l) 


. (L-1)K+1 
. (L-l)K+2 


>K(2L+2)ej0K(2L+2)aK(2L+2)-^^ej0K<L,aK(L)-yK(L)-[i7i)]L 


(3.44) 

To  simplify  the  math  somewhat,  it  is  necessary  to  make  the 
assumption  that  the  fading  characteristics  for  each  user 


will  be  approximately  constant  over  two  consecutive  bit 
intervals.  In  other  words,  it  will  be  assumed  that 


aj  (L  + 1)  — a j (L)  = 0 (3.45) 

Note  that  this  assumption  is  justified  in  the  sense  that 
from  equation  (3.44),  the  receiver's  performance  will 
ultimately  depend  upon  how  small  the  quantity  aj(L  + 2)-aj(L) 

is.  Thus,  if  the  assumption  given  in  equation  (3.45)  is  not 
valid,  then  aj(L  + 2)-aj(L)  will  not  be  very  small  and  the 

receiver  is  doomed  to  fail.  It  is  therefore  only  of 
practical  interest  to  study  the  case  in  which  the  assumption 
is  valid  and  to  evaluate  the  receiver's  performance  as  a 


84 


function  of  how  correlated  the  fading  process  of  each  user 
is  over  a span  of  L+2  consecutive  bit  intervals. 

Under  the  above  assumption,  the  first  K elements  in  the 
vector  of  equation  (3.44)  reduce  to  just  the  sum  of  the 
noise  component  and  the  initial  MAI  component;  thus 
equation  (3.44)  simplifies  to 


(wo  a-wAaJ= 

\ e ee  e eel 


*D 


*D 


(L— 1)K+1 
(L-l)K+2 


-7i(L)- 


A/w1(2L+2)e:)<t,;  (2L+2)a1  (2L+2)-1/w1(L)e:"t>l(L)a1  (l)-y,  (L)— 

V w2  (2L+2)e:i412  (2Lf  2)a2  (2l+2)_^/w2  (L)ejlML)a2  (L)-y2  (L)-  1,(1) 


li(l) 


(L— 1)K+1 
(L-l)K+2 


7wK(2lrf2)ej^(2w2)aK(2L+2)-7^ej^<L)aK(L)-YK(L)-[l1(l)]L 


(3.46) 

To  further  evaluate  the  MAI  term  present  after  decorrelating 
the  second  window,  the  vector  given  in  equation  (3.46)  must 
be  pre-multiplied  by  R~*P  • First,  note  that  the  vfc^ 
element  of  the  vector  [r~*Px]  will  be  of  the  form 

K 1-1  _ K— 1 LK 

[RSPXL=ZXc..lRi.><1>P]1+I  Xcv.JRJ-<..-,„<.I(-l)[x]I  (3.47) 

1=2  j=l  1=1  j=(L-l)K+l+l 

where  X is  an  arbitrary  vector  and  as  before  [x]^  denotes 

the  m^h  element  in  the  vector.  Indeed,  it  is  the  form  of 
equation  (3.47)  which  led  to  the  recursion  matrix  for 
consecutive  windows  used  in  the  AWGN  analysis.  Upon 
substituting  X = (weOeAe  - WeOeAe ) into  equation  (3.47),  the 
following  result  is  obtained 
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K— 1 IK  r -I 

+Z  S cVJ^i^1)K.1(-^{^(2*-J+2)Oe(2Lrf2)Ae(2Lff2)— V'^(L)Oe(l^Ae(L)^y1(L)-[l1(l)]^^^| 

(3.48) 

By  re-grouping  terms,  equation  (3.48)  may  be  re-written  as 

[/  » - ^ ,1U=2  K 1-1  K— 1 LK 

R^p(We<DeAe-WeOeAe)  =-Hcv.jR,i(1>Yi(L)-X  XcvdRj-<L-i»K.i(-1)Yi(L) 

K-1  LK  V 1=2  j=l  1=1  j=(L— 1)K+1+L 

^CV,j^j-(L-l)K,l  (— l){we  (2L+2)Oe  (2L+2)Ae  (2L+2)-We  (L)Oe  (L)Ae  (L)} 

l=lj=(L-l)K+l+L 

K 1-1  K— 1 LK  

1)K-X  X 


1=2  j=l 


1=1  j=(L— 1)K+1+L 


Jl+<L-1)K 


(3.49) 

The  first  two  terms  on  the  right  hand  side  of  equation 

(3.49)  are  noise  components.  The  third  term  is  an 
"acquired"  interference  component  introduced  by  the  second 
window,  and  the  last  two  terms  are  due  to  the  "initial" 
interference  introduced  by  the  first  window.  To  more 
clearly  understand  the  significance  of  each  of  the  two  kinds 
of  interference  terms  present  in  equation  (3.49),  suppose 
that  there  was  no  fading  in  the  channel.  Then  equation 

(3.49)  would  still  contain  the  "initial"  interference  terms 
since  the  first  window  would  still  have  no  previous 
estimates  available  with  which  it  could  attempt  to  subtract 
off  its  out-of -window  MAI.  The  "initial"  interference  terms 
are  thus  inevitable  for  a windowed  decorrelator  operating  in 
any  channel.  Now,  however,  the  "acquired"  interference  term 
would  not  be  present  in  equation  (3.49)  since  in  the  absence 
of  fading  aj(2L  + 2)-aj(L)  = 0 (for  small  L)  . Note  that  the 
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"initial"  interference  term  has  been  reduced  after  the 
second  window  has  been  decorrelated  due  to  the 
multiplication  by  R^P . This  is  almost  precisely  the  same 
situation  which  occurred  in  the  AWGN  channel  and  thus  the 
"initial"  interference  term  will  become  smaller  and  smaller 
as  more  and  more  windows  are  processed  until  it  finally 
becomes  negligible.  The  only  difference  between  the 
behavior  of  the  "initial"  interference  term  in  the  AWGN 
channel  and  the  flat  fading  channel  is  that  in  the  former, 
the  worst  case  estimation  errors  of  the  ut^1  window  were 
determined  as  a function  of  the  worst  case  estimation  errors 
of  the  (u-l)t^1  window,  and  so  on.  As  a consequence  of  using 
worst  case  estimation  errors,  the  sign  preceding  the 
recursion  matrix  A was  always  positive  and  thus  equation 
(2.53)  was  eventually  arrived  at.  For  the  present  fading 
channel  analysis,  worst  case  estimates  are  not  being  used 
since  the  estimates  are  no  longer  obtained  simply  by  taking 
the  absolute  value  of  the  appropriate  element  in  the 
decorrelated  output  vector.  Consequently,  an  equivalent 
bound  on  the  estimation  errors  cannot  be  found  as  in  the 
AWGN  analysis.  Now  the  actual  errors  are  being  used,  and 
decorrelating  the  first  window  yields  an  estimate  which  is 
then  subtracted  from  the  second  window's  matched  filter 
output  vector,  thus  reversing  the  sign  on  the  "initial" 
interference  term.  The  estimate  obtained  from  decorrelating 
the  second  window  is  in  turn  subtracted  from  the  third 


87 


window's  matched  filter  output  vector,  thus  reversing  the 
sign  on  the  "initial"  interference  term  once  again.  This 
process  repeats  itself  and  as  a result  the  fading  channel 
counterpart  to  equation  (2.53)  is  therefore 

i;(u)=(-du-1au-1i;(i)  (3.50) 

where  li(u)  is  a K element  vector  containing  the  last  K 
elements  of  the  interference  vector  I1(u),  and  I^l)  is 

defined  in  equation  (3.42).  The  recursion  matrix  A is 
defined  by  its  (P/q)*-*1  element  as 

LK  q-1 

ap,q  — SCp,j^j-(L-l)K,q(_ ^)+  5rfCp,j^j,q(-*-)  (3.51) 

j=(L-l)K+l+q  j=l 

The  matrix  A in  equation  (3.50)  should  not  be  confused  with 
the  matrix  of  out-of -window  fading  coefficients  Ae  . 

The  real  difference  between  the  AWGN  channel  and  the 
fading  channel  is  the  presence  of  the  "acquired" 
interference  term  in  equation  (3.49).  This  term  is  a result 
of  the  fact  that  the  fading  characteristics  of  each  user 
immediately  before  the  second  window  are  not  the  same  as 
those  immediately  after  the  second  window.  For  compactness 
define  an  interference  vector  due  to  the  second  window  as 

Ll2^Jv=S  Xcv.jRj-<L-1,K,i(-1)[We(2L+2)<^e(2L+2)Ae(2L+2)-We(L)<De(L)Ae(L)] 

V 1=1  j=(L-l)K+l+l 

(3.52) 

Note  that  this  "acquired"  interference  term  has  the  same 
form  as  the  "initial"  interference  term  did  after  the  first 
window  was  decorrelated.  Consequently  the  second  window's 
"acquired"  interference  term  will  obey  the  same  relationship 
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at  the  first  window's  "initial"  interference  term,  except 
now  the  sign  on  this  term  will  be  positive  after  the  second 
window,  negative  after  the  third  window,  positive  after  the 
fourth  window,  and  so  on.  Thus 

I2  (u)  = (— l)u~2  Au-2 1^(2),  u > 2 (3.53) 

The  structure  of  equation  (3.49)  will  be  the  same  for 
each  additional  window  processed  by  the  receiver;  that  is, 
after  the  nfch  window  is  decorrelated,  the  MAI  term  present 
in  the  decorrelated  output  vector  will  contain  an  "acquired" 
interference  component  introduced  by  the  n^  window  as  well 
as  reduced  versions  of  the  "acquired"  interference  terms 
introduced  by  each  of  the  preceding  (n-2)  windows.  There 
will  be  also  be  a reduced  version  of  the  "initial" 
interference  term  introduced  by  the  first  window.  From 
equations  (3.50)  and  (3.53),  the  total  MAI  contained  in  the 
nth  window's  intermediate  decision  statistic  vector  will  be 


[r^  p(weoeAe  - weoeAe  )]^"  = (-lr1  An_i  1;  (1) + (-i)n-2  An-2 1;  (2) 


(3.54) 


+ (-D 


n-3  , n-3 


i;(3)+  •••  +I„(n) 


Equation  (3.54)  can  be  written  more  compactly  as 
[R^p(we<DeAe  -WeOeAe)]^"  =!>)  + £(— l)n-jAn-jI^) 


(3.55) 


j-l 


where  the  "acquired"  interference  term  introduced  by  the  m^h 
window  is  given  by 

r j^i  lk 

— ^ ^Cv,j^j-(L-1)K,1  (-l)[we(m(L+l)>De(m(L+l))Ae(m(L+l)) 


l=lj=(L-l)K+l+l 

-We((m-l)(L+l)>De((m-l)(L+l))Ae((m-l)(L+l))] 
and  [l»  (m)j  contains  the  last  K elements  in  ^Im(rn)j. 


(3.56) 
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All  that  remains  to  be  considered  are  the  noise 
components  in  equation  (3.49).  Since  these  terms  are 
identical  to  the  "initial"  interference  terms  in  equation 


(3.49)  with  [i^l)]^  replaced  by  [y(l)J  , they  will  behave  in 

exactly  the  same  manner.  The  MAI  term  in  each  additional 
window  will  contain  a reduced  noise  component  from  each  of 
the  preceding  windows.  Therefore  the  total  noise  component 
contained  in  the  MAI  term  of  the  n^h  window's  intermediate 
decision  statistic  vector  will  be 


[*« PfW-°eA«  ~ WeOA  )jU“n  = X (-l)n-j  An-j Nd (j)  (3.57) 

L N 'JNoise 

where  the  noise  vector  N^j)  is  defined  as 

ML 

K-l  LK  

^ v,j^  j— (I*— 1)K,1  L 

1=1  j=(L— 1JK+1+1  L 

Finally,  by  combining  the  MAI  component  and  the  noise 
component,  the  vt^1  element  of  the  MAI  vector  in  the  n^ 
window's  intermediate  decision  statistic  vector  is  then 
given  by 


>1 


(3.58) 


l(L-l)K+l 


-W.®.A.)£  ’ =[l„(n)]v 


.j»l 


n-1 


X(-irjAn-jNj(j) 


(3.59) 


U=i  Jv 

Equations  (3.14)  through  (3.20)  can  now  be  applied  to 
obtain  an  expression  for  the  probability  of  error  resulting 
after  decorrelating  the  n^h  window.  Here  the  variable  X in 
equation  (3.14)  will  be  given  by 
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X = -v/wl(1-l)e-',-‘1':iat(i-l)d1(i-l)  + Yk(i-l) 

j=i 

and  the  variable  Y will  be  given  by 
If  = t/w*  (i)e’*‘<‘)a,  (i)d,  (i)  + yk  (i) 

^)+I(-irJA"-j[Jie)+s^)] 

j-1 


(3.60) 


(3.61) 


It  is  clear  from  the  preceding  equations  and  discussion 
that  unlike  its  AWGN  counterpart,  the  fading  channel 
windowed  decorrelator  will  encounter  an  irreducible  error 
floor  due  to  the  presence  of  the  "acquired"  interference 
introduced  by  the  current  window.  In  other  words,  even  if  a 
"steady-state"  condition  could  be  reached  in  which  the 
influence  of  the  out-of -window  MAI  encountered  from  all 
previous  windows  was  negligible,  the  "acquired"  MAI  in  the 
current  window  will  introduce  an  irreducible  error  floor, 
thereby  limiting  the  receiver's  performance  and  resulting  in 
a near-far  resistance  of  zero.  It  therefore  becomes 
important  to  determine  what  the  near- far  spread  will  be  for 
the  fading  channel  windowed  decorrelator.  The  best  possible 
performance  for  the  fading  channel  windowed  decorrelator 
with  imperfect  MAI  estimates  will  correspond  to  the  scenario 
in  which  all  "acquired"  MAI  from  previous  windows  is 
negligible,  and  thus  only  the  MAI  introduced  by  the  current 
window  will  be  considered.  In  reality,  the  effects  of  the 
MAI  from  the  several  most  recent  windows  would  probably  not 
be  negligible,  but  this  analysis  will  provide  an  upper  bound 
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on  the  receiver's  performance.  Ignoring  the  MAI  introduced 
by  previous  windows,  equation  (3.59)  becomes 

[R^p(weOeAe  -W>eAe)£'"  =[ljn)]v  (3.62) 

For  notational  compactness,  define 

hk(i)  = Vw  k(i)ej<Ml)ak(i)  (3.63) 

and 

S =[R^P(weOeAe  -W>.A.)]  (3.64) 

Then  from  equations  (3.56),  (3.63),  and  (3.64),  equation 

(3.62)  will  become 

_ K— 1 LK 

M7  = S Z c v.i : R (-1)[hi  (n(L  + 1))-hi  ((n  - 1)(L  + 1))]  (3.65) 


1=1  j=(L-l)K+l+l 

The  decision  variables  X and  Y are  then  given  by 

X = hk  (i  - l)dk  (i  - 1)  + yk  (i  - 1)  + [s]"_r  (3.66) 

Y = hk  (i)dk  (i)  + Yk  (i)  + [s]"  (3.67) 

Finally,  the  covariances  and  Mxy  are  then 

=2wk(i-l)Rak(0)+a2cv_2>v_2  +E{[s]nv2[s*TL}  (3.68) 

Mw  = 2wk(i)RSk  (0)+  a2cv  v + E{[s];[s*Tfv  } (3.69) 

= 2Vwk(i-l)Vwk(i)Rak  (T)  + a2cv_2  v + E{[S];_K[S*T];}  (3.70) 


where  the  notation  X*T  represents  the  conjugate  transpose 
operation.  By  symmetry,  . Next,  the  quantity 

E{(XIv-k[S*T]v-k}  f°un<3  to 

E{[CKp'T.»}=E{(lj=(L!titic.-K,JRJHw*.i<-1){hiWL+l))-h1((n-lXL+l))]^ 

(K— 1 LK  'N  'J 

Z Zcv-K,qRq-a-i,K.p(-l)[hp(n(L+l))-hp((n-lXL+l))] 
p=lq=(L-l)K+l+p  ) J 


(3.71) 
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which  can  be  re-written  as 

, - K— 1 LK  K— 1 LK 

EPLk[S  ]v_kJ  S X X ^^-K.jRjHL-lIKa^l^v-K.qRq-U-lJK.pC-!)' 

l=lj=(L-l)K+l+lp=lq=(L-l)K+l+p 

Ejfh,  (n(L+l))-h1  ((n-l)(L+l))][hp  (n(L+l))-hp  ((n-l)(L+l))]} 

(3.72) 

Since  the  fading  processes  of  different  users  are 
statistically  independent, 

E-jhi (m)h*  (n)}  = 2^  (m)  ^ (n)RSi  (m - n)8lp  (3.73) 


where  the  delta  function  is  defined  as 

s f a = b 

ab  [0,  otherwise 

Using  equation  (3.73),  equation  (3.72)  simplifies  to 


(3.74) 


e{PLP‘,L  }= ‘I'.sK  R.,  «L + 2W] 

1=1 


LK 


-|2 


5rfCv-K,jRj-(Ii-l)K,l(  1) 

j=(L-l)K+l+l 


(3.75) 


where  it  has  been  assumed  for  simplicity  that  the  received 
energies  are  constant  over  the  span  of  L+2  bits;  that  is, 

wk((n-  1)(L  + 1))  = wk(n(L  + 1))  = wk  (i)  (3.76) 

This  assumption  will  be  valid  for  small  window  lengths  since 
the  wk  represent  received  energies  in  the  absence  of 
fading.  Similarly, 


E{[E]"[E'T]"}=4Sw1(i)[R.,  (°)-R.,  (fr+2W] 

1=1 


-\2 


rfCv,jRj-(L-l)K 


,x(-D 


(3.77) 


j=(L— 1)K+1+1 


and 


r . R-i  lk 

E{[sL[^-T]:}=4SwI(iK(0)-R.1((L  + 2)r)]  I;  cv.K,RjH„MJ(-l). 


1=1 

LK 


j=(L-l)K+l+l 


ZCv.qRq-l 


(3.78) 


(L-DK+l+l 


(-D 


q=(L— DK+l+l 
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The  final  expression  for  the  k^h  user's  probability  of  error 


is  then  given  by  equation  (3.15),  with  equations  (3.75), 
(3.77),  and  (3.78)  appropriately  substituted  into  equations 
(3.68),  (3.69),  and  (3.70).  Note  that  the  resulting 

probability  of  error  is  a function  of  the  bit  position 
within  the  window.  Averaging  the  above  mentioned 
expression  over  all  bit  positions  results  in  an  average 
probability  of  error 


Recall  that  the  average  probability  of  error  given  in 
equation  (3.79)  is  a lower  bound  due  to  the  fact  that  only 
the  out-of-window  MAI  introduced  by  the  current  window  is 
considered.  In  order  to  determine  how  tight  this  bound  is, 
a computer  simulation  of  an  asynchronous  DS-CDMA  system 
employing  DPSK  modulation  through  a flat  Rayleigh  fading 
channel  was  created.  In  the  simulation,  each  user 
transmitted  100,000  bits  through  the  channel  and  used  the 
windowed  decorrelating  receiver  to  decode  the  received  bits. 
The  total  number  of  errors  was  counted  and  then  divided  by 
100,000  to  determine  the  average  probability  of  error.  To 
verify  that  the  simulation  was  working  properly,  the  number 
of  users  was  set  equal  to  one  and  the  resulting  probability 
of  error  performance  was  compared  with  equation  (3.39).  The 
results  are  listed  in  Table  3.1  and  plotted  in  Figure  3.8, 
from  which  the  validity  of  the  simulation  is  confirmed. 


v-v  + K 


(3.79) 
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Table  3.1  - Comparison  of  Computer  Simulation  Results  with 

True  Single  User  Error  Rates 


Average  SNR 
(dB) 

Analytical  Prob. 
of  Error 

Simulation  Prob. 
of  Error 

0 

0.25 

0.2512 

5 

0.120126 

0.12233 

10 

0.04545 

0.046217 

15 

0.015327 

0.014079 

20 

0.0049505 

0.0050397 

25 

0.0015762 

0.0014799 

30 

0.0004995 

0.00043997 

35 

0.00015806 

0.00015999 

Figure  3.8  - Comparison  of  Computer  Simulation  Results  with 
True  Single  User  Error  Rates 


Figure  3.9  compares  the  simulation  with  the  lower  bound  for 
the  case  of  two  active  users  with  cross-correlations 
Ri2(0)=%  and  R12(l)=>4  . A nineteen  bit  wide  window  was  used, 
resulting  in  a 10%  overhead  for  each  user.  To  model  the 
fading,  both  users  were  assumed  to  be  traveling  at  65  mph, 
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transmitting  a bit  rate  of  20  kbps,  and  using  a carrier 
frequency  of  900  MHz.  User  two  was  designated  as  the 
desired  user  and  given  an  average  received  signal- to-noise 
ratio  of  30  dB.  As  can  be  seen  from  Figure  3.9,  the 


Figure  3.9  - Comparison  of  Analytical  Lower  Bound  with 
Computer  Simulation 


simulation  indicates  a larger  probability  of  error  than  the 
lower  bound.  This  is  expected,  of  course,  since  the 
simulation  will  contain  the  out-of-window  MAI  influence  from 
all  previous  windows.  Both  the  simulation  and  the  lower 
bound  indicate  that  the  windowed  decorrelator  can  track  the 
single  user  error  probability  for  near-far  effects  of 
approximately  20  dB  or  less.  For  larger  near-far  effects, 
the  windowed  decorrelator's  error  probability  diverges  from 
that  of  the  single  user  and  eventually  approaches  a value  of 
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one-half.  For  comparison  purposes,  the  performance  of  the 
conventional  DS-CDMA  detector  is  also  displayed  in  Figure 
3.9. 

Figure  3.10  displays  the  probability  of  error  for  user 
two  versus  the  average  signal-to-noise  ratio  of  user  two. 
Here  the  average  near-far  effect  is  held  at  30  dB.  Note 
that  the  lower  bound  is  not  too  far  off  from  the  computer 
simulation.  The  simulation  reveals  an  error  floor  of 


Figure  3.10  - Error  Floors  of  Lower  Bound  and  Simulation  for 
a 30  dB  Near-Far  Ratio 

. -4 

approximately  5x10  whereas  the  lower  bound  predicts  an 
error  floor  of  approximately  2xl0~4.  Figure  3.11  displays 
the  same  scenario  as  in  Figure  3.10,  except  now  the  average 
near-far  effect  is  increased  to  40  dB.  The  result  is  that 
the  lower  bound  is  not  as  tight  as  before:  the  simulation 
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_3 

shows  an  error  floor  of  approximately  5x10  whereas  the 
lower  bound  predicts  an  error  floor  of  approximately  1x10  3. 
This  increase  is  expected,  however,  since  the  MAI  acquired 
by  past  windows  will  become  less  negligible  as  the  severity 
of  the  near-far  effect  increases. 


Figure  3.11  - Error  Floors  of  Lower  Bound  and  Simulation  for 
a 40  dB  Near-Far  Ratio 


The  lower  bound  derived  in  the  preceding  analysis  may 
be  tightened  by  extending  the  previous  analysis  to  include 
the  effects  of  the  out-of-window  MAI  acquired  from  the 
window  immediately  preceding  the  current  window,  since  this 
window's  MAI  influence  will  be  the  next  largest  after  that 
of  the  current  window.  Equation  (3.59)  will  then  become 

[s]"=[R^p(we<DeAe-WeOeAe)]“  " =p^)]v-[Ai;_1(n-l)]v+[Nn_1(n-l)]v  (3.80) 
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Using  equation  (3.58)  to  expand  out  the  noise  term  on  the 
right  hand  side  of  equation  (3.80),  the  expression  for  the 
noise  component  in  the  MAI  term  contained  in  the  nth 
window's  intermediate  decision  statistic  vector  becomes 


[n^Oei-d]  = -XXc^R^CDY.CCn-lKL  + D-l) 

v 1=2  j=l 

K— 1 LK 


(3.81) 


- X Xc  vjR  i-(L-i,K.i  (-l)Yi  ((n  - 1)(L  + 1)  - 1) 

1=1  j=(L— 1JK+1+1 

The  middle  term  on  the  right  hand  side  of  equation  (3.80) 
represents  the  MAI  acquired  from  the  (n-l)st  window; 
expanding  out  this  term  results  in 


[1  IV  K — 1 LIV 

^in-l(n  — 1)1  — Xav,qX  XC<I'-l>K+q,j^j-<L-l)K,l  ’ 

V q=l  1=1  j=(L— DK+l+l 

[h,  ((n  - 1)(L  + 1))  - hx  ((n  - 2XL  + 1)  - 1)] 


(3.82) 


The  first  term  on  the  right  hand  side  of  equation  (3.80) 
represents  the  out-of-window  MAI  acquired  by  the  n1-*1  window 
and  as  before  it  is  given  by  equation  (3.65) . 

Recall  that  the  average  probability  of  error  for  the 
fading  channel  windowed  decorrelator  is  given  by  equation 
(3.79).  The  decision  variables  X and  Y are  again  given  by 
equations  (3.66)  and  (3.67)  and  the  covariances  M^,MW/  and 

M^  are  again  given  by  equations  (3.68)  through  (3.70).  The 
only  difference  in  the  analysis  of  the  second  lower  bound  is 


that  the  expression  for  [s]  given  in  equation  (3.80) 

contains  more  terms  than  it  did  for  the  first  lower  bound, 
thereby  resulting  in  more  complicated  expressions  for  the 
covariances.  Repeating  the  same  procedure  used  to  derive 


99 


equations  (3.71)  through  (3.78)  yields 

3{pl rLHi". Mk «»- r., «l + 2wf 

J 1=1  [jNl-dk+i+i 


h2 


— K , j R j— ( L— 1 )K  ,1  ( 1) 


K— 1 J”  K LK 

+ 4Xwi(i)[R.l(°)-R.1((L  + 2)r)1  £av_K,q  £ 

1=1  1q=l  j=(  L— 1)K 

K— 1 


n2 


C(L-l)K+q,jRj-<L-l)K,l( 


,-llK+l+l 


- wi(i)[2Rai  ((L  + 1)T)-  Rai  (T)  - Rai  ((2L  + l)r)]- 


1=1 


K LK 

LK 

£av-K,q  £C(L-l)K+q,jRj-(L-l)K,l  (~ ^) 

^Cv-K,s^s-(L-1)K,1  (”^-) 

_q=l  j=(L— 1)K+1+1 

_s=(L— 1)K+1+1 

K K 


1-1  r-1 


££C(L-l)K+l,(L-l)K+r  £ICv-K.jCv-K,Rj.l(1)Rs.r(l) 

1=2  r=2  j=l  s=l 

K-l  K-l  LK  LK 

££C(L-l)K+l,(L-l)K+r  £ XCv,jCv,sRj-|L-llK,l(_^)Rs-(L-l|K,r(_^) 

1=1  r=l  j=(L-l)K+l+l  s=(L— l)K+l+r 

K K-l  1-1  LK 

££C<L-l)K+l,(L-l)K+r  £ £Cv,jCv,sRj,l(l)Rs-(L-l)K,r(—  1) 

1=2  r=l  j=l  s=(L-l)K+l+r 

(3.83) 

r — r_  in  1 K_1  LK 

EM-K  KTU=4I  W1  (ifr-x  (°)  - : R.»  «L  + 2>r)]  £c  v— k,r  wl-uk.1  (-1)  • 

C J 1=1  j=(L— 1)K+1+1 

LK  K-l 

£cv,sRs_a_1)K1(-l)  + 4£w1(i)[Rai(0)-Rai((L  + 2)T)]- 


s=(L— 1)K+1+1 
K LK 


LK 


^av-K,q  ^C(L-l)K+q,j^j-(L-l)K,l(  ^)^av,r  ^Cr 


v— K,q 

q=l  j=(L-l)K+l+l 

K-l 


s s— (L-l)K 


a(“D 


r=l  s=(L-l)K+l+l 


“ 2X  wx  (i)[2Rai  ((L  + 1)t) - Rai  (T)  - Rai  ((2L  + 1)t)]- 

i=i 


K 

LK 

~ 

LK 

2^a  v-K,q 

£C(L-l)K+q,jRj-(L- 

XCv,s^s-(L-l)K,l  (“”1) 

_q=1 

j=(L— DK+l+l 

_s=(L— 1)K+1+1 

K K 

1-1  r-1 

l L-l)K+l,(L-l)K+r  ^ / v 
1=2  r=2  j=l  s=l 

K-l K-l  LK 


+ a2X£c(L. 

1=1  r=l 
K K-l 
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(3.84) 

The  expression  for  e|[s][s  T]  | will  be  identical  to  equati 
(3.83),  except  the  subscript  v-K  will  be  replaced  with  v. 
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The  resulting  lower  bound  on  the  average  probability  of 
error  is  then  found  by  using  equations  (3.83)  and  (3.84)  in 
the  expressions  for  M^,  Myy,  and  M^  given  by  equations 

(3.68)  through  (3.70),  which  in  turn  are  substituted  into 
equation  (3.79)  . 

The  second  lower  bound  is  compared  to  the  first  lower 
bound  and  the  simulation  results  in  Figures  3.12,  3.13,  and 
3.14.  In  Figure  3.12,  the  second  lower  bound  is  seen  to 
represent  a substantial  improvement  over  the  first  lower 
bound:  approximately  half  of  the  total  out-of -window  MAI 

influence  due  to  all  previous  windows  is  accounted  for  by 
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Figure  3.12  - Comparison  of  Second  Lower  Bound  with  First 
Lower  Bound  and  Simulation 


the  MAI  incurred  from  the  most  previous  window.  In  Figure 
3.13,  the  second  lower  bound  predicts  an  error  floor  of 
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approximately  3.3x10  which  compares  very  favorably  with 

-4 

the  error  floor  of  5x10  indicated  by  the  simulation. 

Similarly,  Figure  3.14  also  shows  a substantial  improvement 

in  the  second  lower  bound  over  the  first:  an  error  floor  of 

. -3 

approximately  2 . 3x10  is  now  predicted  which  is  not  too  far 

-3 

off  from  the  simulation's  value  of  5x10 

The  second  lower  bound  could  be  tightened  by  extending 
the  previous  analysis  to  include  the  effects  of  out-of- 
window MAI  from  even  more  past  windows;  however,  such  an 
analysis  would  soon  become  analytically  intractable. 
Therefore,  the  second  lower  bound  will  be  used  as  an 
approximation  for  the  average  probability  of  error  of  the 


Figure  3.13  - Error  Floors  of  Second  Lower  Bound,  First 
Lower  Bound,  and  Simulation  for  an  Average  Near-Far  Ratio  of 
30  dB 
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Figure  3.14  - Error  Floors  of  Second  Lower  Bound,  First 
Lower  Bound,  and  Simulation  for  an  Average  Near-Far  Ratio  of 
40  dB 


windowed  decorrelator.  Of  course,  another  way  to  view  the 
results  presented  in  this  section  is  that  the  lower  bounds 
serve  to  further  confirm  the  validity  of  the  computer 
simulation;  therefore,  the  simulation  will  be  used  as  well 
in  obtaining  performance  results  for  the  windowed 
decorrelator  in  a flat  Rayleigh  fading  channel. 

Unfortunately,  equations  (3.83)  and  (3.84)  are  so 
unwieldy  that  they  preclude  the  possibility  of  using 
equation  (2.10)  to  solve  for  the  fading  channel  windowed 
decorrelator's  asymptotic  efficiency.  The  results  shown  in 
Figures  3.9  and  3.12,  however,  are  sufficient  to  allow  the 
performance  of  the  windowed  decorrelator  to  be  completely 
characterized;  moreover,  they  suggest  that  the  overall 
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system  structure  proposed  previously  will  require 
modification  in  order  for  the  fading  channel  windowed 
decorrelator  to  be  near-far  resistant. 

Modified  Multiuser  Receiver  Structure 

Figures  3 . 9 and  3 . 12  show  that  the  windowed 
decorrelator  is  capable  of  tracking  the  bit  error  rate  of  a 
single  user  up  to  a certain  level  of  near-far  effect  (about 
20  dB  for  the  system  considered  in  Figure  3.9);  for  greater 
near- far  effects,  the  windowed  decorrelator  loses  this 
ability.  This  observation  leads  to  the  conclusion  that  the 
multiuser  receiver  structure  presented  earlier  requires  some 
modification  in  order  to  become  near-far  resistant;  namely, 
an  average  power  control  scheme  is  required  which  would 
function  to  keep  the  average  received  power  levels  of  all 
users  to  within  a certain  range  of  one  another  (20  dB  in 
Figure  3.9).  At  first,  the  idea  of  using  a power  control 
scheme  in  conjunction  with  the  windowed  decorrelator  may 
seem  self-defeating,  since  the  original  goal  of  this  work 
was  to  overcome  the  near- far  problem  without  the  use  of 
conventional  power  control;  however,  this  combination  is  not 
self-defeating  for  three  important  reasons.  First,  the 
fading  channel  windowed  decorrelator  requires  an  average 
power  control  scheme,  not  a bit-by-bit  power  control  scheme 
as  is  required  by  the  conventional  DS-CDMA  receiver.  Thus 
the  rate  at  which  the  power  control  updates  must  be  made 


104 


will  be  much  less  for  the  windowed  decorrelator,  thereby- 

overcoming  the  necessity  of  updating  powers  during  fast 

fading.  Secondly,  the  conventional  DS-CDMA  receiver 

requires  a difference  of  0 dB  between  all  received  signals 

for  optimum  performance;  as  the  difference  in  received 

powers  increases,  the  performance  of  the  conventional 

detector  degrades.  On  the  contrary,  as  demonstrated  by 

Figures  3.9  and  3.12,  the  windowed  decorrelator's 

performance  remains  unaffected  by  changes  in  the  received 

powers  as  long  as  they  are  all  within  the  required  range  (20 

dB  in  this  case) . Of  course,  the  exact  power  control  range 

required  will  depend  upon  the  number  of  active  users  and  the 

magnitudes  of  their  signature  waveform  cross-correlations. 

Finally,  even  if  perfect  bit-by-bit  power  control  was 

achievable  for  the  conventional  detector,  there  would  still 

be  an  irreducible  error  floor  due  to  the  presence  of  the  MAI 

component  in  each  user's  decision  statistic.  As  seen  in 

Figure  3.9,  for  a near-far  effect  of  0 dB  the  conventional 

-2 

detector  can  only  achieve  a bit  error  rate  of  4.5x10  . On 

the  contrary,  as  long  as  the  average  received  energies  are 
kept  within  20  dB  of  one  another,  the  windowed  decorrelator 
is  able  to  remove  the  MAI  and  approach  the  single  user  error 
rate . 

To  summarize,  the  windowed  decorrelating  receiver  for 
flat  Rayleigh  fading  channels  must  be  used  in  conjunction 
with  an  average  received  power  control  scheme.  The 
estimates  of  each  user's  received  power  level  may  be 
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obtained  by  the  MAI  estimation  procedure  illustrated  in 
Figure  3.3  and  thus  no  additional  computational  complexity 
is  required  to  track  the  received  power  levels.  The 
multiuser  receiver  can  simply  compare  the  magnitudes 
obtained  from  the  estimation  procedure  and  then  send  out 
appropriate  control  signals  to  the  users,  instructing  them 
to  adjust  their  transmitter  output  powers  accordingly.  The 
block  diagram  of  the  fading  channel  windowed  decorrelator 
incorporating  power  control  is  shown  in  Figure  3.15. 


Figure  3.15  - Block  Diagram  of  Fading  Channel  Windowed 
Decorrelator  Incorporating  Power  Control 
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Numerical  Examples 

Numerical  examples  are  presented  in  this  section  to 
highlight  the  performance  of  the  windowed  decorrelator  in 
flat  Rayleigh  fading  channels.  The  results  for  two  users 
displayed  in  Figures  3.16  through  3.20  were  obtained  from 
the  computer  simulation  discussed  previously.  The  results 
for  more  than  two  users  presented  in  Figures  3.21  through 
3.25  were  obtained  from  the  second  lower  bound  derived 
earlier  in  this  chapter. 

The  performance  of  the  windowed  decorrelator  as  a 
function  of  the  window  length  was  investigated;  the  results 
are  displayed  in  Figure  3.16  for  window  lengths  of  19  bits 
(10%  overhead)  and  5 bits  (33%  overhead) . Here  zeroth 


Figure  3.16  - Performance  of  the  Windowed  Decorrelator  as  a 
Function  of  the  Window  Length  with  Zeroth  Order  Estimation 
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order  estimation  was  employed.  For  error  rates  of  1x10  or 
less,  both  window  sizes  yield  almost  identical  performance. 
For  the  two  user  system  considered  here,  the  power  control 
scheme  would  need  to  keep  the  average  received  power  levels 
to  within  about  20  dB  of  one  another  in  order  to  achieve  an 
error  rate  approaching  that  of  a single  user. 

Figure  3.17  compares  the  zeroth  order  and  first  order 
MAI  estimation  strategies  when  a window  length  of  nineteen 
bits  is  used.  Interestingly,  for  error  rates  of  1x10  3 or 
less  there  is  no  noticeable  difference  in  the  receiver's 
performance;  thus,  the  zeroth  order  technique  would  be 
preferred  since  it  results  in  computational  savings  over  the 
use  of  linear  extrapolation. 


Figure  3.17  - Comparison  of  Zeroth  Order  and  First  Order 
Estimation  Techniques  for  a Window  Length  of  Nineteen  Bits 
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Figure  3.18  compares  the  zeroth  order  and  first  order 
MAI  estimation  strategies  when  a window  length  of  five  bits 
is  used.  Now  there  is  a significant  difference  in 
performance.  Note  that  when  the  zeroth  order  strategy  is 
used,  the  power  scheme  must  keep  the  averaged  received  power 
levels  to  within  about  32  dB  of  one  another  to  obtain  a 
1x10  3 error  rate.  The  first  order  strategy  only  requires 
the  average  received  power  levels  to  be  within  about  40  dB 
to  achieve  the  same  level  of  performance,  and  so  for  this 
scenario  the  first  order  technique  would  be  preferred. 

Figure  3.19  compares  window  lengths  of  five  bits  and 
nineteen  bits  when  the  first  order  estimation  strategy  is 
employed.  For  a window  length  of  nineteen  bits,  a power 

_3 

control  range  of  about  32  dB  is  required  to  achieve  a 1x10 


Figure  3.18  - Comparison  of  Zeroth  Order  and  First  Order 
Estimation  Techniques  for  a Window  Length  of  Five  Bits 
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Figure  3.19  - Performance  of  the  Windowed  Decorrelator  as  a 
Function  of  Window  Length  When  First  Order  Estimation  is 
Used 


error  rate,  whereas  for  a window  length  of  five  bits  this 
range  becomes  about  40  dB.  The  conclusion  from  Figures  3.15 
through  3.18  is  that  for  smaller  window  lengths,  the  first 
order  MAI  estimation  strategy  (linear  extrapolation)  can 
significantly  relax  the  level  of  average  power  control 
required  by  the  system.  As  the  window  length  increases,  the 
advantage  offered  by  the  first  order  strategy  becomes 
smaller . 

The  irreducible  error  floors  encountered  by  the  fading 
channel  windowed  decorrelator  for  various  average  near- far 
effects  are  displayed  in  Figure  3.20.  As  can  be  seen,  an 

-3 

error  rate  of  1x10  or  less  may  be  achieved  as  long  as  the 
average  near- far  ratio  is  kept  to  within  roughly  30  dB  or 


less . 
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Figure  3.20  - Error  Floors  of  the  Windowed  Decorrelator  for 
Various  Average  Near-Far  Ratios 


The  performance  of  the  windowed  decorrelator  as  a 

function  of  the  number  of  active  users  is  examined  in 

Figures  3.21  through  3.25.  In  Figure  3.21,  Gold  codes  of 

length  seven  chips  were  used  for  the  signature  sequences. 

The  results  indicate  that  average  near- far  ratios  of  less 

than  35  dB,  25  dB,  and  17.5  dB  are  required  in  order  for  the 

-3 

receiver  to  obtain  a 1x10  error  rate  for  two,  three,  and 

four  active  users,  respectively.  In  Figure  3.22,  Gold  codes 

of  length  31  chips  were  used  for  the  signature  sequences, 

resulting  in  lower  cross-correlation  magnitudes.  Here  it  is 

seen  that  average  near-far  ratios  of  less  than  57.5  dB,  42 

dB,  and  22  dB  are  required  in  order  for  the  receiver  to 

-3 

maintain  a 1x10  error  rate  for  two,  three,  and  four  active 
users,  respectively.  The  conclusion  from  Figures  3.21  and 
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3.22,  therefore,  is  that  even  when  a decorrelating  strategy 
is  used,  it  is  still  important  to  choose  signature  waveforms 
with  low  cross-correlations;  specifically,  lower  cross- 
correlation magnitudes  will  result  in  a relaxation  of  the 
power  control  range  required  by  the  communication  system. 

Figures  3.21  and  3.22  were  obtained  by  assuming  that 
each  user  was  moving  at  a speed  of  65  mph,  corresponding  to 
a cellular  telephone  scenario.  For  many  applications, 
however,  such  as  personal  communication  systems  or  indoor 
wireless  networks,  the  mobile  units  will  be  traveling  at 
much  lower  speeds.  Since  the  rapidity  of  the  fading 
processes  is  a function  of  the  speed  of  the  mobile  units 
[28] , the  fading  will  occur  at  a much  slower  rate  in  such 
communication  systems.  As  a result  of  the  slower  fading, 
the  estimates  of  the  future  MAI  should  become  better,  and 
the  performance  of  the  windowed  decorrelator  would  be 
expected  to  improve.  Figure  3.23  compares  a computer 
simulation  of  the  time  autocorrelation  function  of  a typical 
flat  Rayleigh  fading  process  for  user  speeds  of  65  mph  and 
15  mph.  Indeed,  for  the  user  traveling  at  the  slower  speed, 
the  fading  process  remains  highly  correlated  over  a much 
wider  time  span.  Figure  3.24  repeats  the  system  of  Figure 
3.19  (Gold  codes  of  length  seven  chips)  but  with  the  users 
moving  at  15  mph  instead  of  65  mph.  Likewise,  Figure  3.25 
repeats  the  system  of  Figure  3.20  (Gold  codes  of  length 
thirty-one  chips)  but  with  user  speeds  of  15  mph.  Clearly, 
there  is  a noticeable  improvement  in  performance  as  the 
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Figure  3.23  - Comparison  of  the  Autocorrelation  Functions  of 
Flat  Rayleigh  Fading  Processes  for  User  Speeds  of  65  mph  and 
15  mph 


Figure  3.24  - Error  Rate  of  the  Windowed  Decorrelator  vs. 
the  Number  of  Active  Users  for  Gold  Codes  of  Length  Seven 
Chips  and  User  Speeds  of  15  mph 
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Figure  3.25  - Error  Rate  of  the  Windowed  Decorrelator  vs. 
the  Number  of  Active  Users  for  Gold  Codes  of  Length  Thirty- 
one  Chips  and  User  Speeds  of  15  mph 


speed  of  the  mobile  units  decreases;  specifically,  the 
required  power  control  range  is  relaxed.  It  is  reasonable, 
therefore,  to  expect  that  the  windowed  decorrelator  would 
tend  to  be  more  robust  in  a personal  communication  system  or 
indoor  wireless  network  than  in  a cellular  telephone  system. 


CHAPTER  4 

IMPLEMENTATIONAL  CONSIDERATIONS 

In  this  chapter,  two  important  practical  considerations 
related  to  the  implementation  of  a windowed  decorrelating 
receiver  are  examined.  First,  the  sensitivity  of  the 
receiver  to  errors  in  the  timing  information  is 
investigated.  Secondly,  the  computational  requirements  for 
computing  the  inverse  of  the  cross-correlation  matrix  are 
examined. 

Sensitivity  to  Propagation  Delay  Estimation  Errors 

In  the  previous  chapters,  the  performance  analyses  of 
the  windowed  decorrelator  for  AWGN  and  frequency 
nonselective  Rayleigh  fading  channels  assumed  that  the 
multiuser  receiver  had  perfect  knowledge  of  the  propagation 
delay  from  each  user  to  the  base  station.  In  practice  it 
will  not  be  possible  for  the  receiver  to  obtain  perfect 
estimates;  even  if  the  MAI  corruption  could  be  removed,  the 
background  noise  would  still  result  in  estimation  errors. 

It  is  therefore  of  interest  to  determine  how  the  windowed 
decorrelator's  performance  will  be  affected  by  imperfect 
timing  information. 

Imperfect  timing  information  will  affect  the  windowed 
decorrelator  in  two  ways:  the  vector  of  matched  filter 
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sampled  outputs  will  be  different  than  before,  and  the 
computed  inverse  cross-correlation  matrix  will  no  longer 
represent  the  inverse  of  the  true  cross-correlation  matrix. 
It  should  be  noted  that  the  operation  of  windowing  the 
received  sequences  does  not  in  any  way  affect  the 
decorrelator's  sensitivity  to  timing  errors;  in  other  words, 
the  windowed  decorrelator  should  have  the  same  sensitivity 
to  timing  errors  as  does  the  full-length  decorrelator.  Also, 
since  the  decorrelating  operation  itself  is  the  same  for 
both  the  AWGN  and  flat  Rayleigh  fading  channels,  the 
following  analysis  will  be  for  a full-length  decorrelator 
operating  in  an  AWGN  channel.  A similar  analysis  has  been 
presented  by  Strom  [59],  although  in  that  work  Monte  Carlo 
simulations  were  used  to  obtain  results.  The  goal  here  will 
be  to  derive  an  analytical  expression  for  the  asymptotic 
efficiency  of  a decorrelating  receiver  as  a function  of  the 
severity  of  the  timing  estimation  errors. 

For  convenience,  the  notation  introduced  by  Strom  [59] 
will  be  adopted.  The  system  is  essentially  described  by  the 
same  notation  as  in  chapter  2,  but  now  only  one  window  of 
length  M bits  is  transmitted  and  decorrelated . Of  course, 
pilot  bits  are  no  longer  employed  since  no  windowing  is 
being  done.  The  received  signal  at  the  base  station  is 


where  n(t)  is  white  Gaussian  noise  with  power  spectral 
density  c2  and 


(4.1) 
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S(t'b)=  ZXVwk(i)bk(i)sk(t-iT-'0  <4. 2) 

i=l  k=l 

where  as  before  wk(i)  is  the  received  energy  of  the  kth  user 
during  the  ith  interval,  bk(i)e{±l}  is  the  polarity  of  the 
kth  user's  ith  bit,  and  sk(t)  is  the  kth  user's  RF  signature 
waveform  normalized  to  unit  energy  and  zero  outside  the  time 
interval  [0,T]. 

Let  the  time  delay  analyzer  shown  in  Figure  (2.7) 

th 

produce  an  estimate  of  the  k user's  propagation  delay 
denoted  by  Tk  . It  will  be  assumed  that  xke[0,T).  The  ith 

fch 

sampled  output  of  the  k user's  matched  filter  now  becomes 


(i+l)T+Tk 

yk(i)  = Jr(t)2k(t  — iT  — Xk)dt 

iT+'xk 

i+l)T+tk 

= J s(t,  b)ik(t  - iT  - Xk)dt  + nk(i) 

iT+Tk 


(4.3) 


where 


(i+l)T+Tk 

nk(i)  = j*n(t)sk(t  - iT  - xk)dt 

iT+Tk 


is  a zero-mean  Gaussian  random  variable  with  cross- 
correlation given  by 

oo 

E[nk(i)n1(j)]  = a2  Jsk(t-iT-xk)s1(t- jT-xJdt 


(4.4) 


(4.5) 


Next,  the  following  two  cross-correlation  matrices  are 
defined  by  their  (k,j)th  elements 


/V  oo 

Ki  (!)=  |sk(t-tk)l.(t  + lT-Xj)dt  (4.6) 

oo 

Rkj(1)=  |sk(t-xk)i.(t  + lT-x.)dt 


(4.7) 
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so  that  equation  (4.5)  may  be  rewritten  as 

E[nk(i)n1(j)]  = a2Rkl(i- j) 


(4.8) 


Due  to  the  time-limited  nature  of  the  signature  waveforms, 
it  is  again  clear  that  Rkj(l)  = Rkj(l)  = 0 for  |1|  > 1 . Upon 


substituting  for  sft,b),  equation  (4.3)  becomes 

K 

y*  (i)  = nk  (i)  + X bi  (i  - R, cl  (1) 

1=1 

+Zbi(i)Vwi(i)Rki(0)  + Xb1(i  + 1)Vwi(i  + 1)^ki(-1) 

1=1  1=1 

where  the  fact  that  bk(0)  = bk(M  + 1)  = 0 has  been  used. 


(4.9) 


Following  the  same  steps  as  in  chapter  2,  the  complete 
vector  of  matched  filter  sampled  outputs  can  be  written  as 

Y = RWb  + n (4.10) 


where  Y = [yT(l),yT(2 ),•••, yT(M)]T  , y(i)  = [yx(i),y2(i)  ,--,yK(i)]T  , 
n = [nT(l),nT(2 ),•••, nT(M)]T  , n(i)  = [n1(i),n2(i),-»,nK(i)]T  , 


W = diag[w(l),W(2 ),•••, W(M)] , W(i)  = diag[A/w1(i),^/w2(i),---,iywK(i)j  , 


R = 


R(0)  R(l)  0 0 

R(l)  R(0)  R(— 1)  i 

0 R(l)  R(0)  0 

1 ’*■  > ^ 

0 0 R(l)  R(0) 


and 


The  noise  vector  n is  a Gaussian  zero-mean  vector  with 
a covariance  matrix  given  by 

E[nnT]  = o2R  (4.11) 

where  the  matrix  R is  defined  as  R with  the  appropriate 
substitutions . 

Recall  that  the  goal  of  the  multiuser  receiver  is  to 
determine  the  vector  of  transmitted  polarities  b from  the 
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vector  of  matched  filter  sampled  outputs  Y given  by 

equation  (4.10).  The  matrix  R is  unknown  to  the  receiver, 

however,  since  its  elements  are  functions  of  the  true  time 

delays  through  equation  (4.7);  thus,  the  receiver  cannot 
~ - 

compute  R . Of  course,  the  multiuser  receiver  does  not 
know  that  its  timing  information  is  in  error,  and  thus  it 
will  form  its  decisions  according  to 

(4.12) 


b = sign  R 1RWb  + R 


It  is  immediately  clear  from  equation  (4.12)  that  since 
R-1R*I,  the  multiuser  receiver  will  not  be  able  to 
completely  decorrelate  the  received  signals  and  will 
consequently  have  a near-far  resistance  of  zero.  It  is  of 
interest,  then,  to  determine  the  near- far  spread  obtainable 
by  the  receiver.  Define  the  matrix 


_i 

E = R XR  = 


and  the  vector 


Then  equation  (4.12)  becomes 


el.l 

ei,2  ei,MK 

e2,l 

e2,2  '*’  e2,MK 

eMK,l 

P • • • P 

^MK,2  ^MK,MK 

II 

R_1n 

b = sign[EWb  + y] 


th  th 

and  the  estimate  of  the  i bit  of  the  k user  is 


bk(i)  = sign  Xev.i[wbi  + [y]v 


Li=i 


= sign 


ev.vVwk(i)bk(i)  + Xev.i[wb]1  +[Y]V 


(4.13) 


(4.14) 


(4.15) 


(4.16) 
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where  as  before  the  vector  element  v is  related  to  the  bit 
index  i through  equation  (2.14)  and  the  notation  [x]x 
denotes  the  1th  element  in  the  vector  x . 

The  conditional  probability  of  bit  error  will  then  be 
given  by 


pk(i|b)  = Pr{bk(i)  > °|bk(i)  = -!'b} 

MK  

[y]v  > ev,vVMi)  - Xev.i[wbl 


= Pr- 


1=1 

l*v 


(4.17) 


where  the  notation 


r»i 

L -Jm,n 


denotes  the  element  in  the  mth  row 


th 


and  n column  of  the  matrix  R 1 . The  average  probability 
of  error  is  found  by  averaging  equation  (5.17)  over  all 
possible  vectors  b with  bk(i)  = -1  . Doing  so  yields 


pk(ilb)  = ^=r_X  Q 


be{±l}“ 

bfc(i)=-l 


(4.18) 


v 7 

Since  computing  equation  (4.18)  will  not  be  practical 
for  large  MK,  the  performance  measure  of  choice  becomes 
asymptotic  efficiency.  In  the  limit  as  o —> 0,  the  summation 
in  equation  (4.18)  will  be  dominated  by  the  Q-function  with 
the  smallest  argument,  thus 
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1 im  P„  ~ Q 

CT — >0  e'avg 


V 

where  the  vector  VL  has  been  defined  as 


ev.v 

MK 

Vw*(i)_Z 

hv^KlJ 

[«-] 

v,v 

(4.19) 


H,  = [7wx (1) , 7 W2  (1) , • • •, ^/wK  (1) , (2) , ^/w2  (2) , • • •, 7wk (2)/"  •/ VW1  (M) » Vw2  (M) / 1 ' ' */ Vwk (M)f  • 


From  equation  (2.10)  the  asymptotic  efficiency  for  the 
windowed  decorrelator  with  imperfect  timing  estimates  is 


T|k(i)  = max2 


0,-s 


R 


-l 


f 

r — i 

MK 

e , WvJ 

|ev,v|-2. 

cv,lLv,bJi 

Vwic(i) 

1=1 

) 

V 1*V 

(4.20) 


Of  course,  if  the  receiver's  timing  estimates  were  perfect 
then  E=I,  R-1  = cvv,  and  equation  (4.20)  would  reduce  to 

L -*v,v 

the  expression  found  by  Lupas  [38]  for  the  full-length 
decorrelator . 

To  obtain  some  insight  into  the  sensitivity  of  the 
decorrelator  to  timing  errors,  equation  (4.20)  was  evaluated 
for  a three  user  system  employing  Gold  codes  of  length  31 
chips.  Each  user  transmitted  a total  of  five  bits.  Three 
of  the  possible  33  Gold  codes  were  chosen  at  random,  and 
three  true  propagation  delays  were  randomly  chosen  to  be 
X1  =17.1921  chips,  x2  = 23.1971  chips,  and  X3  = 29.6162  chips. 

The  asymptotic  efficiency  of  user  two  was  computed  for  the 
third  transmitted  bit  (v=8);  this  was  an  arbitrary  choice 
since  there  is  no  reason  why  the  timing  sensitivity  should 
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be  dependent  upon  the  bit  index,  even  though  the  value  of 
asymptotic  efficiency  will  be.  The  results  are  displayed  in 
Figure  4.1  for  three  different  timing  error  scenarios: 

Tk  = Tk  ± 0.1  chips,  Tk=Tk±0.01  chips,  and  Tk=Tk±0.001  chips. 


Figure  4.1  - Sensitivity  of  Decorrelator  to  Timing  Errors 
for  One  Particular  Timing  Configuration 


As  can  be  seen  from  Figure  4.1,  when  the  timing  estimates 
are  perfect,  the  decorrelator  has  an  asymptotic  efficiency 
of  0.967  and  a near-far  spread  of  infinity.  When  the  timing 
estimates  of  all  users  are  off  by  one- thousandth  of  a chip, 
the  decorrelator's  asymptotic  efficiency  drops  slightly  to 
0.95  and  its  near-far  spread  becomes  approximately  67  dB. 
When  the  timing  estimates  of  all  users  are  off  by  one- 
hundredth  of  a chip,  the  decorrelator's  asymptotic 
efficiency  drops  to  0.925  and  its  near-far  spread  becomes 
approximately  45  dB.  Finally,  when  the  timing  estimates  of 
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all  users  are  off  by  one-tenth  of  a chip,  the  decorrelator's 
asymptotic  efficiency  drops  to  0.7  and  its  near-far  spread 
becomes  approximately  24  dB. 

The  results  displayed  in  Figure  4 . 1 are  for  one 
particular  timing  configuration;  that  is,  the  three  chosen 
Gold  codes  and  the  three  chosen  true  propagation  delays.  Of 
course,  different  Gold  codes  and  different  propagation 
delays  will  result  in  different  cross-correlations  which  in 
turn  will  result  in  different  values  of  asymptotic 
efficiency.  It  is  therefore  of  interest  to  average  the 
decorrelator's  timing  sensitivity  over  all  possible  timing 
configurations;  however,  since  an  infinite  number  of  timing 
configurations  exist,  the  asymptotic  efficiency  was  averaged 
over  1000  different  timing  configurations.  For  each 
configuration,  three  of  the  33  possible  Gold  codes  were 
chosen  at  random  along  with  three  true  propagation  delays. 
The  resulting  average  asymptotic  efficiency  for  the  third 
bit  of  user  two  is  displayed  in  Figure  4.2  for  the  same 
timing  error  scenarios  as  in  Figure  4.1.  Interestingly,  the 
results  shown  in  Figure  4 . 1 are  fairly  representative  of  the 
average  timing  sensitivity  of  the  decorrelator.  Note  from 
Figures  4.1  and  4.2  that  for  this  three  user  system,  as  long 
as  the  received  energies  are  kept  within  about  25  dB  of  one 
another,  then  a nonzero  asymptotic  efficiency  is  achievable 
when  each  user's  propagation  delay  estimate  is  off  by  as 
much  as  a tenth  of  a chip.  The  conclusion,  therefore,  is 
that  the  use  of  a transmitter  power  control  scheme  in 
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Figure  4.2  - Sensitivity  of  Decorrelator  to  Timing  Errors 
Averaged  Over  1000  Timing  Configurations 


conjunction  with  the  windowed  decorrelator  will  relax  the 
accuracy  requirements  of  the  timing  information.  Thus,  even 
in  an  AWGN  channel,  a practical  decorrelating  receiver  will 
need  to  include  an  average  power  control  scheme. 

The  fact  that  the  decorrelating  strategy  is  sensitive 
to  errors  in  the  timing  information  should  not  be 
surprising,  as  a decorrelating  receiver  is  essentially 
performing  the  operation  of  inverse  filtering.  The  cross- 
correlation matrix  RLK  may  be  thought  of  as  the  transfer 

matrix  of  the  multi-input/multi-output  channel.  In  order  to 
successfully  compute  the  inverse  of  the  channel  transfer 
matrix,  perfect  timing  information  will  be  necessary.  An 
important  conclusion,  therefore,  is  that  for  channels  in 
which  it  is  not  feasible  to  obtain  reliable  timing 
information,  a method  other  than  a decorrelating  strategy 
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should  be  considered  for  demodulating  the  received  signals. 
Of  course,  any  demodulation  method  which  includes  a matched 
filter  bank  inherently  requires  knowledge  of  the  system 
timing  configuration. 

Inversion  of  the  Cross-correlation  Matrix 

The  windowed  decorrelator  needs  to  compute  the  inverse 
of  the  symmetric  block  tri-diagonal  Toeplitz  matrix  RLK , 
which  is  of  dimension  LK  by  LK  where  L is  the  window  length 
in  symbols  and  K is  the  number  of  active  users.  As  noted  by 
Lupas  [61],  the  straightforward  inversion  of  RLK  by  Cholesky 
decomposition  or  Gaussian  elimination  will  require  on  the 
order  of  2(LK)3  operations.  To  obtain  some  numerical 
results  on  the  computational  complexity  of  performing  the 
matrix  inversion,  the  computer  mathematical  modeling  package 
MATLAB  was  utilized.  A useful  feature  of  MATLAB  is  its 
ability  to  count  the  number  of  floating  point  operations 
(flops)  it  performs.  Additions  and  subtractions  are  each 
one  flop  if  real  and  two  if  complex.  Multiplications  and 
divisions  count  one  flop  each  if  the  result  is  real  and  six 
flops  if  it  is  complex. 

A two  user  asynchronous  system  was  constructed  and  the 
matrix  RLK  was  inverted  for  various  window  lengths.  The 
results  are  summarized  in  Table  4.1.  MATLAB  was  then  used 
on  the  data  in  Table  4.1  to  determine  a polynomial  curve  fit 
of  order  three.  The  result  is 
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Table  4.1  - Number  of  Flops  Required  to  Invert  RLK  for  Two 
Users  with  Various  Window  Lengths 


L 

LK 

Flops 

2 

4 

179 

3 

6 

495 

4 

8 

1047 

5 

10 

1899 

6 

12 

3115 

7 

14 

4759 

8 

16 

6895 

No.  of  flops  = 1.3333(LK)3  + 5.5(LK)2  - L6  6 6 7(LK)  - 1 (4.21) 

Note  that  equation  (4.21)  is  an  exact  curve  fit,  not  an 
approximation.  In  a practical  multiple  access  communication 
system,  many  users  may  be  simultaneously  active.  For 
instance,  if  K = 35  users  and  L = 8 bits  then  from  equation 
(4.21)  the  computation  of  the  inverse  of  R280  will  require 
29,699,334  flops. 

One  simple  way  to  reduce  the  number  of  flops  required 
to  compute  R“^  is  to  choose  a small  window  length  . If  a 
window  length  of  L = 4 bits  is  used  in  the  previous  example, 
then  the  number  of  flops  required  to  invert  R140  is 
3,766,141,  or  approximately  one-eighth  of  the  number  of 
flops  required  when  the  window  length  was  twice  as  wide. 
Keeping  the  window  length  small  will  also  result  in  better 
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estimates  of  future  MAI  since  the  received  energies  and 

phases  will  have  had  less  time  to  change  since  being 

estimated.  Of  course,  the  major  disadvantage  in  choosing  a 

small  window  length  is  that  the  percentage  of  overhead  will 

increase  due  to  the  presence  of  the  pilot  bits. 

The  formula  presented  in  equation  (4.21)  is  based  upon 

straightforward  computation  of  the  matrix  inverse.  An 

iterative  algorithm  has  been  found  which  can  reduce  the 

3 2 

number  of  computations  from  approximately  2 (LK)  to  n (LK)  , 
where  n is  a number  which  depends  upon  the  eigenvalue  spread 
of  the  matrix  RLK  [61].  Other  methods  also  exist  for 

reducing  the  computational  complexity  of  inverting  Toeplitz 
matrices  [62]. 

Finally,  it  should  be  recalled  that  the  matrix  RLK  only 

requires  inverting  when  the  multiuser  receiver  detects  a 
change  in  either  the  number  of  active  users  or  their 
relative  time  delays.  For  indoor  wireless  communication 
networks  and  personal  communication  systems,  these 
quantities  may  remain  static  for  much  of  the  duration  of  a 


user's  message. 


CHAPTER  5 
CONCLUSIONS 

Summary  of  Time-Windowed  Decorrelating  Receivers 

In  this  work,  a unique  method  of  implementing  a 
multiuser  decorrelating  receiver  for  asynchronous  DS-CDMA 
channels  was  presented.  Unlike  the  original  full-length 
decorrelator,  the  time-windowed  decorrelator  proposed  here 
has  a small  decoding  delay  (on  the  order  of  L bits,  where  L 
may  be  chosen  as  desired) , only  needs  to  store  LK  elements 
in  memory  (where  K is  the  number  of  active  users) , and  only 
needs  to  invert  a matrix  of  dimension  LK  by  LK. 

In  chapter  two,  a windowed  decorrelating  multiuser 
receiver  was  proposed  and  analyzed  for  AWGN  channels.  The 
front-end  of  the  receiver  is  a matched  filter  bank.  The 
receiver  stores  L consecutive  vectors  of  matched  filter 
sampled  outputs,  thus  forming  a "window".  Estimates  of  the 
out-of-window  multiple-access-interference  (MAI)  affecting 
the  current  window  are  then  subtracted  from  the  current 
window  vector.  The  resulting  modified  window  vector  is  then 
decorrelated  (i.e.,  multiplied  by  the  inverse  of  the  system 
cross-correlation  matrix) . Sign  decisions  are  then  made 
upon  the  resulting  decorrelated  output  vector  to  determine 
the  transmitted  polarities.  To  form  the  required  out-of- 
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window  MAI  estimates,  each  mobile  transmitter  inserts  a 
pilot  bit  into  its  transmitted  sequence  after  every  L data 
bits.  Knowledge  of  the  polarities  of  the  pilot  bits  thus 
provides  the  multiuser  receiver  with  knowledge  of  the 
polarities  of  the  out-of-window  MAI  components.  The 
magnitudes  of  the  last  bit  of  each  user  contained  within  the 
previous  window,  obtained  from  the  previous  window's 
decorrelating  solution,  are  used  as  the  estimates  of  the 
out-of-window  MAI  magnitudes  affecting  the  current  window. 
The  out-of -window  MAI  phases  are  known  since  it  was  assumed 
that  the  carrier  phases  of  the  users  were  being  tracked  for 
the  purpose  of  coherent  BPSK  demodulation. 

The  asymptotic  efficiency  of  the  proposed  window 
decorrelator  was  then  found.  It  was  observed  that  under 
certain  conditions  imposed  upon  the  signature  waveform 
cross-correlations,  the  asymptotic  efficiency  will  converge 
to  a constant  value.  For  the  special  case  of  two  users,  the 
condition  necessary  for  convergence  was  shown  to  be 
identical  to  the  condition  necessary  for  the  existence  of  a 
full-length  decorrelator  in  the  limit  as  the  message  length 
approaches  infinity.  For  the  more  general  K user  scenario, 
the  condition  required  for  convergence  could  not  be 
analytically  derived;  however,  the  results  of  a numerical 
analysis  suggest  that  the  same  condition  as  found  for  the 
two  user  case  also  applies  to  the  K user  case.  It  is 
therefore  conjectured  that  the  windowed  decorrelator  will  be 
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applicable  for  any  communication  system  for  which  the  full- 
length  decorrelator  is  applicable. 

The  asymptotic  efficiency  of  the  windowed  decorrelator 
was  then  evaluated  for  several  numerical  examples  and  seen 
to  provide  a near-far  resistant  solution  up  to  a certain 
level  of  near-far  effect  when  a finite  number  of  windows  was 
processed.  A new  performance  measure,  the  receiver's  near- 
far  spread,  was  defined  as  the  largest  difference  in 
received  energies  for  which  the  receiver  could  operate  at  a 
nonzero  asymptotic  efficiency.  It  was  demonstrated  that  as 
the  number  of  active  users  increases,  the  windowed 
decorrelator  maintains  its  ability  to  achieve  any  desired 
level  of  near- far  spread,  but  the  number  of  windows  which 
must  be  processed  first  increases  accordingly. 

In  chapter  three,  a windowed  decorrelator  was  proposed 
and  analyzed  for  a communication  channel  in  which  each  user 
experiences  statistically  independent  frequency  nonselective 
Rayleigh  fading.  Differential  phase  shift  keying  (DPSK)  was 
chosen  as  the  modulation  technique  so  that  tracking  of  each 
user's  carrier  phase  is  not  required.  The  overall  receiver 
strategy  is  similar  to  that  for  the  AWGN  channel,  but  now 
each  mobile  transmitter  must  insert  two  pilot  bits  between 
every  block  of  L-l  data  bits.  The  first  pilot  bit  is 
grouped  with  the  L-l  data  bits  to  form  a window  L bits  wide. 
Estimates  of  the  out-of-window  MAI  are  subtracted  from  the 
matched  filter  sampled  output  vector,  and  the  resulting 
modified  vector  is  decorrelated.  Standard  DPSK  detection  is 
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then  done  on  the  resulting  decorrelated  output  vector  to 
determine  the  polarities  of  the  transmitted  bit  sequences. 
The  out-of -window  MAI  estimates  are  complex  numbers;  that 
is,  they  must  include  a magnitude  and  phase.  These 
estimates  are  obtained  as  in  the  AWGN  scenario,  only  now  the 
last  K elements  in  the  previous  decorrelated  output  vector, 
corresponding  to  the  first  of  the  two  consecutive  pilot 
bits,  are  used  as  the  MAI  estimates  rather  than  their 
magnitudes . 

Asymptotic  efficiency  was  found  to  be  analytically 
intractable,  so  the  performance  measure  became  probability 
of  bit  error.  It  was  shown  through  numerical  examples  that 
the  windowed  decorrelator  could  approach  the  single  user  bit 
error  rate  up  to  a certain  level  of  average  near- far  effect. 
This  observation  led  to  the  incorporation  of  an  average 
power  control  scheme  into  the  overall  communication  system. 

In  chapter  four,  the  sensitivity  of  the  decorrelating 
strategy  to  errors  in  the  timing  information  was 
investigated.  The  asymptotic  efficiency  of  a full-length 
decorrelator  operating  in  an  AWGN  channel  was  derived  as  a 
function  of  the  error  matrix  resulting  from  incorrectly 
inverting  the  true  cross-correlation  matrix.  It  was  found 
that  the  decorrelating  strategy  is  sensitive  to  timing 
errors;  the  degree  of  sensitivity  increases  as  the  near-far 
effect  increases  and  as  the  magnitude  of  the  timing  errors 
increases.  The  use  of  power  control  will  therefore  also 
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help  to  relax  the  level  of  precision  needed  for  the 
propagation  delay  estimates. 

A brief  discussion  of  the  computational  complexity 
involved  in  computing  the  inverse  of  the  cross-correlation 
matrix  was  also  addressed  in  chapter  four.  It  was 
demonstrated  that  straightforward  inversion  of  the  LK  by  LK 
symmetric  tri-block-diagonal  Toeplitz  matrix  requires  on  the 
order  of  2(LK)3  floating  point  operations.  References  were 
given  which  discuss  more  computationally  efficient 
algorithms,  although  an  in-depth  treatment  of  such 
algorithms  is  beyond  the  scope  of  this  work. 

Areas  For  Future  Work 

In  this  section  several  areas  for  further  work  related 
to  the  windowed  decorrelator  are  suggested.  A fundamental 
requirement  of  any  decorrelating  receiver  is  the  knowledge 
of  the  propagation  delay  from  each  user  (i.e.,  the  system 
timing  configuration) . It  was  assumed  in  this  work  that  an 
ideal  time  delay  analyzer  was  available  to  provide  the 
windowed  decorrelator  with  perfect  timing  information.  In 
practice,  of  course,  such  a device  must  be  built  and  thus 
methods  for  obtaining  knowledge  of  the  timing  configuration 
need  to  be  developed. 

Another  area  for  future  work  is  the  modification  and 
performance  analysis  of  the  windowed  decorrelator  for 
frequency  selective  fading  channels.  RAKE  receivers  are 
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often  suggested  for  use  in  frequency  selective  fading 
channels  [25],  [59].  In  the  multiuser  context,  the  matched 

filter  bank  would  be  replaced  with  a bank  of  RAKE  filters, 
one  for  each  user  [50],  [51].  The  decorrelating  algorithm 

would  then  decorrelate  all  the  received  replicas  from  each 
user,  and  appropriate  recombining  of  the  decorrelated 
replicas  would  then  occur  before  final  decisions  are  made. 

Of  course,  a RAKE  receiver  requires  a lot  of  side 
information  which  may  not  be  readily  available.  Acquisition 
of  any  required  side  information  should  therefore  be 
addressed.  Also,  intersymbol  interference  is  a common 
result  of  digital  communication  systems  operating  in 
frequency  selective  fading  channels  and  would  therefore  need 
to  be  addressed  as  well. 

In  many  communication  networks,  data  and  video 
transmissions  are  present  in  addition  to  voice 
transmissions.  Such  systems  would  typically  allow  the  users 
to  have  different  data  rates,  depending  upon  the  type  of 
information  they  were  transmitting.  Varying  data  rates 
among  the  users  changes  some  of  the  basic  system  assumptions 
made  by  a decorrelating  receiver;  thus,  modifying  the 
decorrelator  to  be  capable  of  operating  in  an  environment 
with  varying  data  rates  should  also  be  investigated. 


APPENDIX  A 

PROOF  OF  EQUATION  (2.40) 


In  this  appendix,  the  upper  and  lower  bounds  given  in 
equation  (2.40)  are  proven.  Although  the  proof  is 
straightforward,  it  is  included  here  for  completeness. 

Claim:  — |l|  < Vw  — |±Vw  ± l|  < |l| 

Proof : The  proof  consists  of  considering  the  four  possible 

combinations  of  the  signs  inside  the  absolute  value. 

Case  I:  +Vw,+I 

The  quantity  inside  the  absolute  value  is  Vw  + I so 
that |Vw  + 1|  = Vw  + 1 and  the  result  of  the  subtraction  is 
then  Vw  - (Vw  + 1)  = -I  . 

Case  II:  — Vw,— I 

The  quantity  inside  the  absolute  value  is  -Vw-I  so 
that  |— Vw  - 1|  = Vw  + |l|  and  the  result  of  the  subtraction 
is  then  Vw  - (Vw  + |l|)  = -|l|  . 

Case  III : +Vw ,-I 

The  quantity  inside  the  absolute  value  is  Vw-I  . 

If  Vw  > |l|  then  | Vw  — 1|  = Vw  - |l|  and  the  result  of  the 
subtraction  is  then  Vw  - (Vw  - |l|)  = +|l|  . 

If  Vw  < |l|  then  |Vw  - 1|  = |l|  — Vw  and  the  result  of  the 
subtraction  is  then  Vw  - (|l|  - Vw ) = 2 Vw  - |l|  < 2|l|  - |l|  = |l|  . 
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Ca.se  XV : — Vw,+I 

The  quantity  inside  the  absolute  value  is  I-Vw  . 
If  I>Vw  then  |l  - Vwj  = |l|  - Vw  and  the  result  of  the 
subtraction  is  then  Vw  — (|l|  — Vw)  = 2 Vw  — |l|  < 2|l|  — |l|  = |l|  . 
If  I<Vw  then  |l  - Vw|  = Vw  - |l|  and  the  result  of  the 
subtraction  is  then  Vw  - (Vw  - |l|)  = |l|  . 


APPENDIX  B 

DIAGONAL I Z AT ION  OF  THE  RECURSION  MATRIX  FOR  TWO  ACTIVE  USERS 

In  this  appendix  it  is  shown  that  the  matrix  A,  defined 
in  equation  (2.52),  will  always  be  diagonalizable  for  the 
case  of  two  active  users.  For  convenience  a window  length 
of  three  bits  is  chosen.  The  matrix  A will  be 
diagonalizable  if  its  determinant  is  nonzero  [60];  therefore 
det[A]  — |c5  6p21  ||c6  ^21 1 ~ I1—  5,i P 21 11*“ 6,6  P21 1 ^ 0 ( B . 1 ) 

which  produces  the  condition 

|c5,6|c6,iHcJK6|  (B  • 2 ) 

where  it  has  been  assumed  that  |p21|*0  since  this  would 

result  in  the  matrix  A being  identically  zero  and  thus  the 
windowed  decorrelator  would  already  be  near- far  resistant 
from  equation  (2.55).  Expanding  out  the  equation  RR_1=I, 
where  I is  the  identity  matrix,  yields 
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It  is  evident  from  equation  (B.3)  that 

Pl2C5,l  +C6.1  =0  (B . 4 ) 

and 

Pl2^5,6  *-6,6  = 1 (B  . 5 ) 
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Solving  equations  (B.4)  and  (B.5)  for  c61  and  c66, 

respectively,  and  substituting  the  results  into  equation 
(B.2),  the  following  condition  is  obtained 

|C5,6Pl2C5,l|^|C5,l  — C5,iPi2C5,g|  (B.  6) 

which  is  satisfied  as  long  as 

c5,i*°  (B . 7 ) 

To  proceed  further,  it  will  be  necessary  to  determine  c51  as 
a function  of  the  cross-correlations  p12  and  p21  . From 
equation  (B.3)  it  is  seen  that 

Cl,l  + Pl2C2.1  =1  (B . 8 ) 


and 


Pi2ci,i  +c2,i  +P21C3.1  =°  (B . 9 ) 

Solving  equations  (B.8)  and  (B.9)  for  cla  and  c21, 

respectively,  and  substituting  the  result  of  equation  (B.8) 
into  equation  (B.9)  produces 


c 


2,1 


Pl2  P 21^*  3,1 

pL-i 


(B.10) 


Next,  from  equation  (B.3)  it  is  clear  that 

P21^"2,l  ^3,1  Pl2<-4,1  — 0 (B.ll) 


Substituting  equation  (B.10)  into  equation  (B.ll)  and 
solving  for  c31  yields 


C3.1  = 


P 12  (p  12  l) 
I-P12-P21 


Cu  + 


P12P 


21 


I-PL-P21 


(B. 12) 


Repeating  this  procedure  by  solving  for  c41  in  terms  of  c51 
and  then  using  the  relationship  given  by  equation  (B.4)  to 
eliminate  c6  3 , the  following  result  is  obtained 
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C 


P 12  3Pl2  P21  ~^Pl2  ^Pl2P21  ^P21  1 


(B . 13 ) 


Applying  the  condition  that  c51*0  given  by  equation  (B.7), 

it  is  obvious  that  the  matrix  A will  be  diagonalizable  if 
and  only  if 


Since  the  cross-correlations  are  functions  of  the 
relative  delay  between  the  received  signals  of  user  1 and 
user  2,  and  since  the  delay  is  defined  along  a continuum  of 
values,  the  converse  of  equation  (B.15)  will  occur  with 
probability  zero.  In  the  event  that  p12p21=0,  one  of  three 
possible  situations  has  occurred:  p12  = 0 , p21  = 0 , or 

Pi2=p2i=0-  If  either  of  the  latter  two  situations  occurs, 
then  the  windowed  decorrelator  will  be  near-far  resistant 
for  the  duration  of  the  occurrence  since  the  matrix  A will 
be  identically  zero.  If  the  first  situation  should  arise, 
then  from  equations  (B.13)  and  (B.4)  c51=0  and  c61=0  so 
that  the  matrix  A will  then  become 


C5.6P21 

O 

1 

C6,6P21 

O 

1 

To  proceed,  it  will  be  necessary  to  solve  for  c56  in  terms 
of  p12  and  p21  . This  can  be  accomplished  by  repeating  the 


(B . 14 ) 


which  is  equivalent  to  the  condition 

P12P21  * 0 


(B.15) 


procedure  used  to  obtain  equation  (B.13)  but  using  the 
identity  R_1R  = I.  Doing  so  results  in 
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c Pl2  ^Pl2  p2lPl2~t~Pl2 . 17. 

5,6  P12  “ 3pJ2  — p21  ~2pi2Ph  +3pi2  + 2p21  — 1 

From  equation  (B.17),  it  is  clear  that  c56  will  equal 
zero  for  p12  = 0 ; thus,  from  equation  (B.16)  the  matrix  A 
will  identically  be  zero  and  again  the  windowed  decorrelator 
will  be  near-far  resistant.  Therefore,  for  the  case  of  two 
users,  the  matrix  A will  always  be  diagonalizable  except  for 
the  trivial  case  when  p12p21  = 0 , for  which  the  matrix  A is 

identically  zero  and  the  windowed  decorrelator  is  near- far 
resistant . 


LIST  OF  REFERENCES 


1.  W.  C.  Lee,  "Overview  of  Cellular  CDMA, " IEEE  Trans,  on 
Vehicular  Technology,  Vol . 40,  No.  2,  pp.  291-302,  May 
1991. 

2.  K.  S.  Gilhousen,  I.  M.  Jacobs,  R.  Padovani,  A.  J. 
Viterbi,  L.  A.  Weaver,  Jr.,  and  C.  E.  Wheatley  III,  "On 
the  Capacity  of  a Cellular  CDMA  System,"  IEEE  Trans,  on 
Vehicular  Technology,  Vol.  40,  No.  2,  pp.  303-312,  May 
1991. 

3.  R.  L.  Pickholtz,  L.  B.  Milstein,  and  D.  L.  Schilling, 
"Spread  Spectrum  for  Mobile  Communications,"  IEEE  Trans, 
on  Vehicular  Technology,  Vol.  40,  No.  2,  pp.  313-322, 

May  1991. 

4.  L.  Milstein,  T.  Rappaport,  R.  Barghouti,  "Performance 
Evaluation  for  Cellular  CDMA, " IEEE  Journal  on  Selected 
Areas  in  Communications,  Vol.  10,  No.  4,  pp.  680-689, 

May  1992. 

5.  M.  Kavehrad,  B.  Ramamurthi,  "Direct-Sequence  Spread 
Spectrum  with  DPSK  Modulation  and  Diversity  for  Indoor 
Wireless  Communications,"  IEEE  Trans,  on  Communications, 
Vol.  COM-35,  No.  2,  pp.  224-236,  February  1987. 

6.  W.  Lam,  R.  Steele,  "Performance  of  Direct  Sequence 
Spread  Spectrum  Multiple-Access  Systems  in  Mobile 
Radio,"  IEE  Proc.-I,  Vol.  138,  No.  1,  pp.  1-14,  February 
1991. 

7.  K.  Johannsen,  "Code  Division  Multiple  Access  Versus 
Frequency  Division  Multiple  Access  Channel  Capacity  in 
Mobile  Satellite  Communication,"  IEEE  Trans,  on 
Vehicular  Technology,  Vol.  39,  No.  1,  pp.  17-26, 

February  1990. 

8.  K.  Gilhouse,  I.  Jacobs,  R.  Padovani,  and  L.  Weaver, 
"Increased  Capacity  Using  CDMA  for  Mobile  Satellite 
Communication, " IEEE  Journal  on  Selected  Areas  in 
Communcations,  Vol.  8,  No.  4,  pp.  503-513. 


140 


141 


9.  M.  Soroushnejad,  E.  Geraniotis,  "Performance  Comparison 
of  Different  Spread- Spectrum  Signaling  Schemes  for 
Cellular  Mobile  Radio  Networks,"  IEEE  Trans,  on 
Communications,  Vol.  40,  No.  5,  pp.  947-956,  May  1992. 

10.  B.  Vojcic,  R.  Pickholtz,  I.  Stojanovic,  "A  Comparison 
of  TDMA  and  CDMA  in  Microcellular  Radio  Channels," 
Proceedings  ICC  '91,  pp.  866-870. 

11.  R.  De  Gaudenzi,  C.  Elia,  R.  Viola,  "Bandlimited  Quasi- 
Synchronous  CDMA:  A Novel  Satellite  Access  Technique 
for  Mobile  and  Personal  Communication  Systems,"  IEEE 
Journal  on  Selected  Areas  in  Communications,  Vol.  10, 
No.  2,  pp.  328-343,  February  1992. 

12.  A.  J.  Viterbi,  "When  Not  to  Spread  Spectrum  --  A 
Sequel,"  IEEE  Communications  Magazine,  Vol.  23,  pp.  12- 
17,  April  1985. 

13.  A.  J.  Viterbi,  "Spread  Spectrum  Communications  --  Myths 
and  Realities,"  IEEE  Communications  Magazine,  pp.  Il- 
ls, May  1979. 

14.  G.  R.  Cooper,  R.  W.  Nettleton,  "A  Spread- Spectrum 
Technique  for  High  Capacity  Mobile  Communications," 

IEEE  Trans,  on  Vehicular  Technology,  Vol.  VT-27,  pp. 
264-275,  November  1978. 

15.  D.  Goodman,  "Second  Generation  Wireless  Information 
Networks,"  IEEE  Trans,  on  Vehicular  Technology,  Vol. 

40,  No.  2,  pp.  366-374,  May  1991. 

16.  M.  Kavehrad,  P.  J.  Mclane,  "Performance  of  low- 
complexity  channel  coding  and  diversity  for  spread- 
spectrum  in  indoor,  wireless  communication, " AT&T 
Technical  Journal,  Vol.  64,  pp.  1927-1965,  October 
1985. 

17.  L.  B.  Milstein,  "Interference  Rejection  Techniques  in 
Spread  Spectrum  Communications,"  Proc . IEEE,  Vol.  76, 
pp.  657-671,  June  1988. 

18.  R.  litis,  L.  B.  Milstein,  "Performance  Analysis  of 
Narrowband  Interference  Rejection  Techniques  in  DS 
Spread  Spectrum  Systems,"  IEEE  Trans,  on 
Communications,  Vol.  COM-32,  pp.  1169-1177,  Nov.  1984. 

R.  Vijayan,  H.  V.  Poor,  "Nonlinear  Techniques  for 
Interference  Suppression  in  Spread  Spectrum  Systems," 


19. 


142 


IEEE  Trans,  on  Communications,  Vol . COM-38,  pp.  1060- 
1065,  July  1991. 

20.  L.  B.  Milstein,  "On  the  Feasibility  of  a CDMA  Overlay 
for  Personal  Communication  Networks,"  IEEE  Journal  on 
Selected  Areas  in  Communications,  Vol. 10,  pp.  655-668, 
May  1992. 

21.  D.  Schiling,  E.  Kanterakis,  "Broadband  CDMA  Overlay  of 
FM  or  TDMA  in  the  Cellular  System, " Conference  Record 
of  the  Communication  Theory  Mini -Conference  (in 
conjunction  with  GLOBECOM  '92),  Orlando,  Florida, 
December  1992,  pp.  61-65. 

22.  L.  A.  Rusch,  H.  V.  Poor,  "Narrowband  Interference 
Suppression  in  CDMA  Spread  Spectrum  Communications, " 
IEEE  Trans,  on  Communications,  Vol.  42,  No.  2,  pp. 
1969-1979,  April  1994. 

23.  J.  Ketchum,  J.  G.  Proakis,  "Adaptive  Algorithms  for 
Estimating  and  Suppressing  Narrowband  Interference  in 
PN  Spread  Spectrum  Systems,"  IEEE  Trans,  on 
Communications,  Vol.  COM-30,  pp.  913-924,  May  1982. 

24.  G.  J.  Saulnier,  "Suppression  of  Narrowband  Jammers  in  a 
Spread  Spectrum  Receiver  Using  Transform-Domain 
Adaptive  Filtering, " IEEE  Journal  on  Selected  Areas  in 
Communications,  Vol.  10,  No.  4,  pp.  742-749,  1992. 

25.  G.  Turin,  "Introduction  to  Spread- Spectrum 
Antimultipath  Techniques  and  Their  Application  to  Urban 
Digital  Radio,"  Proceedings  of  the  IEEE,  Vol.  68,  pp. 
328-353,  March  1980. 

26.  K.  Simon,  J.  Omura,  R.  Scholtz,  B.  Levitt,  Spread 
Spectrum  Communications  Vol  III,  Rockville,  MD: 

Computer  Science  Press  1985. 

27.  J.  Lehnert,  M.  Pursley,  "Multipath  Diversity  Reception 
of  Spread  Spectrum  Multiple-Access  Communications, " 

IEEE  Trans,  on  Communications,  Vol.  COM-35,  pp.  1189- 
1198,  November  1987. 

28.  J.  Proakis,  Digital  Communications,  New  York:  McGraw 
Hill,  1989,  second  edition. 

29.  S.  Verdu,  "Minimum  Probability  of  Error  for 
Asynchronous  Gaussian  Mulitple-Access  Channels, " IEEE 
Trans,  on  Information  Theory,  Vol.  32,  No.  1,  pp.  85- 
96,  January  1986. 


143 


30.  E.  A.  Geraniotis,  M.  B.  Pursley,  "Error  Probability  for 
Direct-Sequence  Spread-Spectrum  Multiple-Access 
Communications --Part  II:  Approximations,"  IEEE  Trans, 
on  Communications,  Vol . COM-30,  pp.  985-995,  May  1982. 

31.  M.  B.  Pursley,  D.  V.  Sawarte,  W.  E.  Stark,  "Error 
Probability  for  Direct-Sequence  Spread- Spec trum 
Multiple-Access  Communications--Part  I:  Upper  and  Lower 
Bounds,"  IEEE  Trans,  on  Communications,  Vol.  COM-30, 
pp.  975-984,  May  1982. 

32.  J.  S.  Lehnert,  M.  B.  Pursley,  "Error  Probabilities  for 
Binary  Direct-Sequence  Spread- Spec trum  Communications 
with  Random  Signature  Sequences,"  IEEE  Trans,  on 
Communications,  Vol.  COM-35,  pp.  87-98,  January  1987. 

33.  P.  K.  Enge,  D.  V.  Sawarte,  "Spread- Spectrum  Multiple- 
Access  Performance  of  Orthogonal  Codes:  Linear 
Receivers,"  IEEE  Trans,  on  Communications,  Vol.  COM-35, 
December  1987. 

34.  L.  F.  Chang,  S.  Ariyavisitakul , "Performance  of  Power 
Control  Method  for  CDMA  Radio  Communication  System, * 

IEE  Electronics  Letters  23rd,  Vol.  27,  No.  11,  May 
1991. 

35.  R.  Diaz,  R.  Agusti,  "Analysis  of  a Fast  CDMA  Power 
Control  Scheme  in  an  Indoor  Environment, * Proceedings 
VTC  '92,  pp.  67-70. 

36.  S.  Soliman,  C.  Wheatley,  R.  Padovani,  "CDMA  Reverse 
Link  Open  Loop  Power  Control,"  Proceedings  GLOBECOM'93, 
pp.  69-73. 

37.  R.  Lupas  and  S.  Verdu,  "Linear  Multiuser  Detectors  for 
Synchronous  Code-Division  Multiple  Access  Channels, " 
IEEE  Trans,  on  Information  Theory,  Vol.  35,  No.  1,  pp. 
123-136,  January  1989. 

38.  R.  Lupas  and  S.  Verdu,  "Near- far  Resistance  of 
Multiuser  Detectors  in  Asynchronous  Channels, " IEEE 
Trans,  on  Communications,  Vol.  38,  No.  4,  pp.  496-508, 
April  1990. 

39.  M.  K.  Varanasi  and  B.  Aazhang,  "Multistage  Detection  in 
Asynchronous  Code-Division  Multiple  Access 
Communications,"  IEEE  Trans,  on  Communications,  Vol. 

38,  No.  4,  pp.  509-519,  April  1990. 


144 


40.  Y.  C.  Yoon,  R.  Kohno,  H.  Imai,  "A  Spread- Spectrum 
Multiaccess  System  with  Cochannel  Interference 
Cancellation  for  Multipath  Fading  Channels,"  IEEE 
Journal  on  Selected  Areas  in  Communications,  Vol.  11, 
No.  7,  pp.  1067-1075,  September  1993. 

41.  S.  Kubota,  S.  Kato,  K.  Feher,  "Inter-Channel 
Interference  Cancellation  Technique  for  CDMA 
Mobile/Personal  Communication  Systems,"  Proceedings  of 
the  IEEE,  pp.  112-117,  1992. 

42.  A.  J.  Viterbi,  "Very  Low  Rate  Convolutional  Codes  for 
Maximum  Theoretical  Performance  of  Spread- Spectrum 
Multiple-Access  Channels,"  IEEE  Journal  on  Selected 
Areas  in  Communications,  Vol.  8,  pp.  641-649,  1990. 

43.  A.  Duel-Hallen,  "Decorrelating  Decision  Feedback 
Multiuser  Detector  for  Synchronous  Code-Division 
Multiple  Access  Channel,"  IEEE  Trans,  on 
Communications,  Vol.  COM-41,  No.  2,  pp.  285-290, 
February  1993. 

44.  P.  Dent,  B.  Gudmundson,  M.  Ewerbring,  "CDMA-IC:  A Novel 
Code  Division  Multiple  Access  Scheme  Based  on 
Interference  Cancellation,"  Proceedings  of  the  IEEE, 
pp.  98-102,  1992. 

45.  Z.  Xie , R.  Short,  and  C.  K.  Rushforth,  "A  Family  of 
Suboptimum  Detectors  for  Coherent  Multiuser 
Communications,"  IEEE  Journal  on  Selected  Areas  in 
Communications,  Vol.  8,  No.  4,  pp.  683-690,  May  1990. 

46.  S.  Miller, "An  Adaptive  Direct-Sequence  Code-Division 
Multiple  Access  Receiver  for  Multiuser  Interference 
Rejection,"  to  appear  in  IEEE  Trans,  on  Communications. 

47.  R.  Kohno,  H.  Imai,  M.  Hatori,  S.  Pasupathy,  "An 
Adaptive  Canceller  of  Cochannel  Interference  for 
Spread- Spectrum  Multiple-Access  Communication  Networks 
in  a Power  Line, " IEEE  Journal  on  Selected  Areas  in 
Communications,  Vol.  8,  No.  4,  pp.  691-699,  May  1990. 

48.  M.  Varanasi,  B.  Aazhang,  "Optimally  Near-Far  Resistant 
Multiuser  Detection  in  Differentially  Coherent 
Synchronous  Channels,"  IEEE  Trans,  on  Information 
Theory,  Vol.  37,  No.  4,  pp.  1006-1018,  July  1991 


145 


49.  M.  Varanasi,  "Noncoherent  Detection  in  Asynchronous 
Multiuser  Channels,"  IEEE  Trans,  on  Information  Theory, 
Vol . 39,  No.  1,  pp.  157-176,  January  1993. 

50.  U.  Fawer,  B.  Aazhang,  "A  Coherent  Diversity  Receiver 
for  Code  Division  Multiple  Access  Communications  over 
Multipath  Fading  Channels,"  Conference  Record  of  The 
Communication  Theory  Mini-Conference  (in  conjunction 
with  GLOBECOM  ’92),  Orlando,  Florida,  December  1992, 
pp.  71-76. 

51.  S.  Vasudevan,  M.  Varanasi,  "Multiuser  Detectors  for 
Asynchronous  CDMA  Communications  over  Rician  Fading 
Channels, " Conference  Record  of  The  Communication 
Theory  Mini-Conference  (in  conjunction  with 
GLOBECOM ' 92 ) , Orlando,  Florida,  December  1992,  pp.  77- 
81. 

52.  Z.  Zvonar,  D.  Brady,  "Suboptimum  Multiuser  Detector  for 
Synchronous  CDMA  Frequency-Selective  Rayleigh  Fading 
Channels, " Conference  Record  of  The  Communication 
Theory  Mini-Conference  (in  conjunction  with 
GLOBECOM ' 92 ) , Orlando,  Florida,  December  1992,  pp.  82- 
86. 

53.  R.  Prasad,  H.  S.  Misser,  A.  Kegel,  "Performance 
Analysis  of  Direct  Sequence  Spread- Spectrum  Multiple 
Access  Communication  in  an  indoor  Rician-Fading  Channel 
with  DPSK  Modulation,"  Electronics  Letters,  Vol.  26, 
pp.  1366-1367,  August  1990. 

54.  S.  Wi jayasuriya,  J.  McGeehan,  and  G.  H.  Norton,  "RAKE 
Decorrelating  Receiver  for  DS-CDMA  Mobile  Radio 
Networks,"  Electronics  Letters,  Vol.  29,  No.  4,  pp. 
395-396,  February  1993. 

55.  Z.  Zvonar,  "Multiuser  Detection  For  Rayleigh  Fading 
Channels, " Ph.D.  Dissertation,  Northeastern  University, 
Boston,  Massachusetts,  September  1993. 

56.  S.  S.  H.  Wi j ayasuriya , J.  P.  McGeehan,  and  G.  H. 

Norton,  "Sliding  Window  Decorrelating  Algorithm  for  DS- 
CDMA  Receivers,"  Electronics  Letters,  Vol.  28,  No.  17, 
pp.  1596-1597,  August  1992. 

57.  R.  Smith,  and  S.  Miller,  "Code  Timing  Estimation  in  a 
Near-Far  Environment  for  Direct  Sequence  Code-Division 
Multiple  Access,"  to  appear  in  Proceedings  of 

MILCOM ' 94 . 


146 


58.  R.  litis, "Joint  Estimation  of  PN  Code  Delay  and 
Multipath  Using  the  Extended  Kalman  Filter,"  IEEE 
Trans,  on  Communications,  Vol . 38,  No.  10,  October 
1990,  pp.  1677-1685. 

59.  E.  Strom,  S.  Parkvall,  and  B.  Ottersten,  "Near-Far 
Resistant  Propagation  Delay  Estimators  for  Asynchronous 
Direct-Sequence  Code  Division  Multiple  Access  Systems," 
Proceedings  of  the  1994  International  Zurich  Seminar  on 
Digital  Communications,  Zurich,  Switzerland,  March 
1994. 

60.  A.  Howard,  Elementary  Linear  Algebra,  Second  edition. 
New  York:  John  Wiley  and  Sons,  Inc.,  1977. 

61.  R.  Lupas,  "Near-Far  Resistant  Linear  Multiuser 
Detection,"  Ph.D.  Dissertation,  Princeton  University, 
Princeton,  New  Jersey,  January  1989. 

62.  T.  H.  Wallace,  "Improving  the  Computational  Efficiency 
of  Subspace  Algorithms  for  Frequency  Estimation  of 
Sinusoidal  Signals,"  Ph.D.  Dissertation,  University  of 
Florida,  Gainesville,  Florida  1994. 

63.  R.  Price,  P.  E.  Green,  Jr.,  "A  Communication  Technique 
for  Multipath  Channels,"  Proceedings  of  IRE,  Vol.  46, 
pp.  555-570,  March  1958. 


BIOGRAPHICAL  SKETCH 


Michael  Feeney  received  the  Bachelor  of  Science  degree 
in  electrical  engineering  from  the  University  of  Missouri- 
Rolla  in  1988  and  the  Master  of  Science  degree  in  electrical 
engineering  from  the  University  of  Florida  in  1989.  From 
1990  to  1991  he  worked  for  Motorola  Inc.  in  the  Personal 
Communication  Systems  Subscriber  Division  developing 
hardware  for  a digital  wireless  PCS  telephone.  Since  1991 
he  has  been  a participant  in  the  Motorola  Distinguished 
Employee-Student  Program  while  working  on  his  Ph.D.  in 
electrical  engineering. 


147 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

^ . 

Scott  L.  Miller,  Chair 
Associate  Professor  of 
Electrical  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Dqcty6r  of  philosophy. 


Leon  W.  Couch,  II,  Cochair 
Professor  of  Electrical 
Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Donald  G.  Childers 
Professor  of  Electrical 
Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

— K yj 

Haniph  A.  Latchman 
Associate  Professor  of 
Electrical  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Richard  E.  Newman-Wolfe 
Assistant  Professor  of 

Computer  and  Information 
Sciences 


This  dissertation  was  submitted  to  the  Graduate  Faculty 
of  the  College  of  Engineering  and  to  the  Graduate  School  and 
was  accepted  as  partial  fulfillment  of  the  requirements  for 
the  degree  of  Doctor  of  Philosophy. 


December,  1994 


Winfred  M.  Phillips 

Dean,  College  of  Engineering 


Karen  A.  Holbrook 
Dean,  Graduate  School 


